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I.  INTRODUCTION 


This  is  the  final  report  on  a  program  to  develop  and  evaluate  a  model  for  an  ionic 
mechanism  of  soot  formation.  The  nucleation  steps  in  soot  formation  remain  unknown,  while  the 
subsequent  steps  are  reasonably  well  understood.  The  nucleation  steps  refers  to  the  formation  of 
an  incipient  soot  particle  from  a  gas  phase  precursor.  In  the  ionic  mechanism  the  gas  phase 
precursor  is  assumed  to  be  the  chemiion,  HCO*;  in  the  free  radical  mechanism  it  is  benzene.  The 
mechanism  of  forming  HCO""  is  well  established;  the  mechanism  of  forming  benzene  is  still  subject 
to  study.  In  the  ionic  mechanism,  a  series  of  ion  molecule  reactions  starting  with  the  precursor 
produce  larger  and  larger  ions  on  the  way  to  soot.  Many  large  soot-like  ions  have  been  observed 
in  sooting  flames.  Our  immediate  concern  is  to  account  for  these  large  ions  by  using  a  reaction 
mechanism  involving  a  series  of  likely  ion-molecule  reactions.  A  major  task  has  been  the 
development  of  adequate  data  bases  for:  (1)  ion  thermochemistry;  (2)  ion-molecule  reaction  rate 
coefficients;  (3)  cation-electron  recombination  coefficients;  and  (4)  diffusion  coefficients.  While 
thermochemical  data  were  available  on  some  large  ions  with  even  numbers  of  carbon  atoms,  there 
were  no  such  data  for  large,  odd  carbon  number  ions.  The  thermochemical  data  base  proved  to  be 
the  most  difficult  to  develop  and  it  also  proved  to  be  the  most  important.  The  rate  coefficients  for 
the  reverse  reactions  were  calculated  from  the  thermochemistry  of  the  species  involved,  and  hence 
depends  on  the  development  of  the  thermochemical  data  base.  For  the  forward  (no  energy  barrier) 
reactions  the  Langevin  theory  of  ion-molecule  reaction  rates  is  consistent  with  experimental  data  at 
ambient  temperatures  and  for  small  ions  but  inadequate  at  flame  temperatures  and  for  large  ions. 
We  thus  extended  this  theory  to  cover  large  ions  and  found  that  a  slightly  positive  temperature 
coefficient  was  introduced  into  the  calculated  rate  coefficients.  However,  use  of  the  extended  theory 
made  very  little  difference  in  the  final  results.  The  other  two  data  bases  are  less  critical  and  were 
easier  to  prepare.  The  data  bases  are  continually  being  updated  as  new  data  appear  in  the  literature 
or  as  our  modeling  effort  requires  data  which  are  unavailable  and  we  estimate  it.  The  model 
predictions  were  compared  with  experimental  data  to  confirm  its  validity. 

In  another  task  the  Sandia  Flame  Code  was  modified  to  handle  ions.  Modifications  were 
made  in  the  code  to  incorporate:  (1)  ion  chemistry;  (2)  ambipolar  diffusion;  (3)  the  ability  to  input 
neutral  and  ion  profiles;  and  (4)  the  ability  to  employ  a  modified  Arrhenius  form  for  ion-molecule 
reaction  (necessary  for  handling  equilibrium  reverse  reactions  without  exceeding  the  Langevin 
reaction  rate). 

The  ionic  mechanism  of  soot  formation  was  applied  to  experimental  acetylene/oxygen  flame 
data  and  compared  with  a  free  radical  mechanism  used  to  model  the  same  data.  The  resultant 
models  were  compared  by  calculating  the  times  to  add  a  carbon  atom  to  the  growing  species  [1,2] 
(the  ionic  mechanism  appears  to  be  faster);  and  by  determining  the  thermodynamic  driving  force 
between  the  two  mechanisms  (the  driving  force  for  the  ionic  mechanism  is  much  greater).  The 
experimental  measurements  made  on  benzene/oxygen  flames  under  a  previous  contract  [3]  have  been 
organized  [4]  to  form  the  experimental  data  base  for  comparing  with  a  computer  model  of  the 
benzene/oxygen  flame.  This  flame  is  much  different  from  the  acetylene/oxygen  flame,  contains 
different  large  ions,  and  will  require  additional  thermochemical  and  reaction  kinetic  estimates.  We 
initiated  an  effort  to  understand  the  relationship  of  observed  soot  formation  and  the 
thermochemistry  of  soot  formation  in  flames  involving  a  variety  of  fuels.  Experimental  data  in  the 
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literature  on  soot  threshold  equivalence  ratios,  soot  yields,  and  acetylene  concentrations  are 
compared  with  thermochemical  equilibrium  calculations  of  these  quantities.  (Appendix  B) 

Two  papers  summarizing  the  work  of  this  program  were  submitted  for  consideration  for 
presentation  at  the  Twenty-Fifth  Combustion  Symposium;  both  were  rejected.  The  reasons  for  the 
rejections  are  interesting,  and  require  a  response,  so  we  use  this  report  as  a  medium  to  record  and 
respond  to  them.  The  two  papers  with  reviewers'  comments  and  our  responses  are  thus  included  as 
Appendices  A  and  B. 


II.  STATEMENT  OF  WORK 


The  ionic  mechanism  of  soot  formation  in  flames  will  be  further  evaluated  and  compared 
with  the  neutral  mechanism  by  pursuing  the  following  phases: 

A.  PHASE  I  -  EXTEND  THE  IONIC  MECHANISM  TO  BENZENE-OXYGEN  FLAMES 

1.  Extend  the  thermodynamic  and  reaction  rate  coefficient  data  base  to  include  those 
species  found  in  benzene/oxygen  flames  which  are  not  present  in  acetylene/oxygen  flames. 

2.  Organize  the  experimental  data  available  on  the  benzene  flame  including  both  neutral 
and  ionic  species. 

3.  Develop  a  detailed  ionic  mechanism  for  the  formation  of  soot  in  benzene  flames  and 
submit  this  to  others  to  run  on  a  large  computer  to  compare  its  agreement  with  experimental  data 
and  to  compare  its  simulation  of  soot  formation  with  the  neutral  mechanism. 

4.  Compare  the  computer  modeling  results  with  experimental  data  and  interpret  the 
results  in  terms  of  the  major  chemical  pathway  to  soot  in  the  benzene  flame. 


B.  PHASE  II  -  MODEL  COAGULATION  AND  AGGLOMERATION 

1.  Extend  the  detailed  mechanism  of  soot  nucleation  to  include  coagulation  of  large  ions 
and  neutrals,  and  charged  with  charged  and  neutral  incipient  particles. 

2.  Program  a  desktop  computer  to  test  the  model  developed  in  Task  I  (in  a  limited 
fashion)  against  experimental  data.  If  warranted,  submit  this  model  to  others  to  test  on  a  mainframe 
computer. 

3.  Interpret  the  results  from  Task  2  in  terms  of  a  simplified  mechanism. 
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C.  PHASE  III  -  DEVELOP  A  THEORY  FOR  LARGE  ION-MOLECULE  REACTIONS 

1.  Extend  the  Langevin  theory  of  ion-molecule  reactions  to  include  large  ions  by 
removing  the  restriction  of  a  point  charge  on  the  ion. 


D.  PHASE  IV  -  COMPARE  THERMODYNAMIC  PREDICTIONS  OF  SOOT  FORMATION 

WITH  EXPERIMENTAL  OBSERVATIONS 

1.  Collect  and  organize  the  literature  data  on  soot  yield  and  acetylene  concentrations 
as  a  function  of  equivalence  ratio. 

2.  Calculate  the  thermodynamic  equilibrium  concentration  of  acetylene  and  soot  as  a 
function  of  equivalence  ratio  for  the  literature  systems  identified  in  Task  1. 

3.  Compare  the  experimental  and  calculated  values  and  interpret  them  in  terms  of 
generalizations  relevant  to  the  mechanisms  of  soot  formation. 


HI.  THERMOCHEMISTRY 


The  thermochemistry  of  cations  continues  to  be  a  major  source  of  uncertainty  in  the  kinetics 
model.  The  thermochemistries  of  several  cations  have  been  revised  according  to  the  most  recent 
literature,  and  we  have  added  a  few  larger  cations  to  the  original  data  base.  As  before,  we  prefer 
evaluated  literature  values  for  enthalpies  of  formation  [5,6];  however,  when  they  are  not  available, 
we  estimate  the  enthalpies  of  formation  from  proton  affinities,  ionization  potentials,  or  theoretical 
calculations.  Heat  capacities  and  entropies  are  calculated  by  either  rigid-rotor-harmonic  oscillator 
statistical  mechanics  (SM)  formulae  [7],  hindered  rotor  harmonic  oscillator  SM  formulae  [8]  or 
Benson's  group  additivity  scheme  [9]  using  the  most  recent  group  values  [10].  SM  formulae  are 
preferred  over  Benson’s  method  when  reliable,  experimental  or  high  level  QM,  structural  and 
vibrational  frequency  information  are  available  for  an  ion.  This  data  is  being  assembled  for  a  major 
publication  which  will  be  submitted  when  the  data  base  is  complete. 


IV.  REACTION  MECHANISM  AND  REACTION  COEFFICIENTS 


The  reaction  mechanism  has  remained  essentially  the  same  as  developed  in  the  previous 
program,  with  some  changes  in  isomers  used  and  with  the  deletion  of  others.  The  only  changes  in 
ion-molecule  reaction  rate  coefficients  resulted  from  use  of  the  extended  Langevin  theory  or  from 
the  effect  of  new  thermochemistry  on  the  reverse  ratio  coefficient. 

To  apply  the  extended  Langevin  theory  to  our  set  of  reactions,  we  used  the  conducting  sphere 
model  of  the  extended  Langevin  theory  [11].  For  acetylene  and  diacetylene,  the  two  major  neutral 
species  reacting  with  ions,  we  calculated  rate  coefficients  as  a  function  of  temperature  for  a  series 
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of  ion  diameters.  The  polarizabilities  of  acetylene  and  diacetylene  (important  parameters  in  the 
calculations)  are:  3.33  x  10“^“*  and  5.99  x  cm^,  respectively.  To  incorporate  the  calculated  rate 
coefficients  into  the  kinetic  model,  we  cast  the  temperature  dependence  in  terras  of  the  classical 

Arrhenius  equation,  (i.e.,  k{T)  =  AT'  e~^).  We  assumed  that  is  small  ( ~  0)  for  exothermic 
ion-molecule  reactions.  The  calculated  rate  coefficients  for  each  ion  diameter  were  then  fit  by 
regression  analysis  to  the  equation: 

logfc  =  logA  +  nlogT  (^) 

This  gave  a  table  of  A's  and  n's  as  a  function  of  ion  diameter.  These  were  then  each  fit  to  simple 
equations  giving  the  dependence  of  A  and  n  on  ion  diameter  (d)  for  both  acetylene  and  diacetylene. 
The  resultant  coefficients  are  given  below. 

Acetylene: 


A  =  (3.37  -I-  12  d)  X  10”,  cmVmol/s 

n  =  0.292  -h  0.0397  x  d,  d  in  nm 

Diacetylene: 

A  =  (2.96  +  9.93  x  d)  x  10”,  cm^mol/s 

n  =  0.298  +  0.033  x  d,  d  in  nm 

For  all  of  the  fits,  the  correlation  coefficients  (r^)  were  0.98  or  greater,  for  d  >  0.7  nm;  there  were 
large  changes  in  values  below  this  ion  diameter  which  is  where  the  extended  Langevin  theory  breaks 
down. 


The  ion  diameters  were  deduced  from  experimental  ion  mobility  (k)  data  for  PAH  molecules 
in  nitrogen  gas  [12].  From  the  mobility  data  and  the  hard  sphere  diameter  of  the  nitrogen  molecule 
(=  0.375  nm.  Ref.  13),  using  Langevin's  ion  mobility  equation  [14], 


K 


A 


( 

1 

\ 


(2) 


we  back  calculated  the  equivalent  ion  diameters  of  the  PAH  ions  in  nitrogen  gas.  Here,  M  is  the 
MW  of  nitrogen  and  m  is  the  MW  of  the  ion,  k  is  the  dielectric  constant  of  nitrogen  and  Pj,  is  the 
gas  density.  ^  is  a  complicated  function  of  a  parameter  (X  -  not  mean  free  path),  gas  pressure,  p, 
/c,  and  the  sum  of  the  neutral  and  ion  collision  radii  (D^).  The  X  parameter  is  given  by: 

^2  ^  ^TipPu  (3) 

(K  -  \)e^ 
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where  e  is  the  electron  charge.  Fitting  these  data  to  various  empirical  functions,  we  found  the 
following  relationship  between  the  deduced  ion  size  and  the  PAH  ion  mass: 

\n{d,A)  =  -1.375  +  2.154  In  (In  (m))  W 

The  correlation  coefficient  (r^)  for  this  fitting  equation  was  0.999-almost  a  perfect  fit  to  the  ion 
diameters. 

For  comparison,  we  plot  ion  diameters  deduced  from  ion  mobility  data  against  ion  diameters 
calculated  for  homogeneous  spheres  having  densities,  p,  of  1.5  or  2.0  g/cm^  -  the  density  range  for 
soots.  The  results  are  plotted  in  Fig.  1  and  they  show  close  agreement.  This  agreement  indicates 
that  the  effective  ion  diameters  derived  from  mobility  data  are  fairly  reliable.  The  fitting  function, 
Eq.  (4),  is  routinely  used  to  calculate  PAH  ion  diameters  for  all  ions  in  our  modeling  work. 

Table  I  includes  a  comparison  of  the  rate  coefficients  by  the  two  theories  for  this  present 
mechanism.  There  is  little  difference  in  the  calculated  ion  concentration  profiles  when  the  rate 
coefficients  from  the  extended  Langevin  theory  are  substituted  for  the  original  Langevin  coefficients. 

In  modeling  experimental  flame  results  using  the  kinetic  scheme  presented  in  Table  I,  we 
input  experimental  profiles  for  and  HjCg"^  (phenyl  cation)  because  we  do  not  have  good 

experimental  profiles  for  the  HCC  ion,  the  initial  ion  produced  by  chemiionization,  from  which  all 
other  ions  are  derived.  This  ion  is  not  observed  in  fuel  rich  flames  because  it  reacts  extremely  fast, 
transferring  a  proton  to  many  neutral  species  including  water.  When  we  calculated  the  concentration 
of  HCO^,  using  the  basic  neutral  reaction  scheme  of  Frenklach,  et.al.  [15]  for  soot  formation  with 
the  addition  of  chemiionization  reactions  to  form  HCO'^  and  experimental  electron  recombination 
rates,  the  calculated  value  was  too  high  compared  to  the  maximum  possible  value  based  on 
experimental  measurements. 

The  question  arises  as  how  to  handle  QHj*.  The  ion  is  identified  as  a  flame  ion  by  mass 
spectrometry  which,  unfortunately,  provides  no  direct  information  on  the  ion  structure.  Two 
isomeric  structures  for  QHj*  have  been  considered  as  potentially  important  [16].  (These  are  the 
linear  propargylium  cation,  denoted  as  C3H3(1)'^,  and  the  cyclic  cyclopropenium  cation,  C3H3(c)'"). 
In  the  mechanism  we  have  chosen  to  identify  the  C3H3(1)'*'  ion  with  the  experimentally  observed  ion. 
Recent  calculations  [17]  show  that  the  cyclic  cation  is  more  stable  than  the  linear  cation  by  116 
kJ/mol;  the  experimental  difference  is  105  kJ/mol  [18].  The  effect  of  this  difference  in  stability 
decreases  with  increasing  temperature;  on  the  basis  of  the  derived  thermochemistry,  at  300  K  the 
equilibrium  ratio  of  linear  to  cyclic  is  5  x  10“^®  but  at  2000  K  the  ratio  is  0.03,  Fig.  2. 

It  has  been  determined  experimentally  that  at  low  temperatures  the  cyclic  isomer  is  not 
reactive  with  acetylene  or  diacetylene,  but  that  the  linear  isomer  is  reactive  with  acetylene  and 
diacetylene  [19].  There  is  no  evidence  that  the  cyclopropenium  cation  isomerizes  to  the 
propargylium  cation  in  the  absence  of  encounter  complexes  with  neutral  molecules,  e.g.,  acetylene 
[17],  Yet  Cameron,  et.al.  [17]  suggest  that  in  sooting  flames  the  cyclic  form  may  not  be  the  most 
abundant.  This  assumes  that  the  linear  ion  is  formed  and  does  not  rapidly  isomerize  to  the  more 
stable  cyclic  ion.  In  the  flame,  the  C3H3*  cation  is  assumed  to  be  formed  in  a  number  of  elementary 
reactions  involving  HCO^[20]. 
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V.  COMPUTER  CODE 


In  our  earliest  work,  the  model  was  run  on  an  IBM  3090  computer  by  M.  Frenklach  and  Hai 
Wang  at  Penn  State  [21].  There  was  a  long  turnaround  time  and  their  code  was  not  adaptable  to 
ions  or  some  of  the  computer  tests  which  we  wanted  to  run.  We  subcontracted  for  one  summer  to 
Prof.  Robert  Brown  at  Iowa  State  University,  who  had  a  code  in  which  ion/electron  diffusion  was 
treated  with  Poissons  equation  [22],  which  treats  ion  transport  more  exactly  than  ambipolar  diffusion, 
an  approximation  we  eventually  used.  There  was  again  a  problem  with  long  turnaround  time,  and 
the  program  not  being  as  flexible  as  we  required.  We  thus  hired  Fokian  Egolfopoulos  at  Princeton 
University  (now  at  University  of  Southern  California)  as  a  consultant  to  modify  the  flame  code 
package,  based  on  Sandia's  CHEMKIN  and  Flame  Codes.  Modifications  were  made  in  the  code  to 
incorporate:  (1)  ion  chemistry;  (2)  ambipolar  diffusion;  (3)  the  ability  to  input  neutral  and  ion 
profiles;  and  (4)  the  ability  to  employ  a  modified  Arrhenius  form  for  handling  rate  coefficients 
calculated  by  equilibrium  constants. 

These  changes  and  the  means  of  using  the  modified  program  are  described  in  the  Symposium 
paper,  Appendbc  A.  To  decrease  run  times,  AeroChem  purchased  a  Silicon  Graphics  IRIS 
workstation.  Run  times  have  been  relatively  short,  from  20  minutes  to  a  few  hours.  Part  of  the 
ability  to  run  on  a  work  station  in  reasonable  times,  rather  than  requiring  a  computer  with  the  power 
of  an  IBM  3090,  is  due  to  using  selected  experimental  mole  fraction  profiles  as  inputs.  Thus  we  were 
able  to  eliminate  a  large  number  of  neutral  species  reactions.  We  note  that  the  neutral  part  of  the 
mechanism  gave  results  which  did  not  agree  with  experimental  profiles.  Clearly  when  we  obtain 
satisfactory  agreement  with  experimental  ion  concentrations  we  will  have  to  consider  a  modified 
neutral  mechanism  in  order  to  model  the  chemiionization  steps. 


VI.  COMPARISON  OF  CALCULATED  AND  EXPERIMENTAL 
ION  CONCENTRATION  PROFILES 


The  availability  of  the  modified  Sandia  Flame  Code  to  handle  ions  and  the  availability  of  a 
workstation  have  greatly  accelerated  the  rate  of  progress  on  this  program.  Thermodynamic 
equilibrium  has  a  much  greater  effect  on  the  ion-molecule  association  reactions  than  expected  (the 
reverse  reaction  rate  constants  are  calculated  from  the  forward  rate  and  the  equilibrium  constants); 
this  shows  up  as  a  stronger  temperature  dependence  than  anticipated.  At  higher  temperatures  the 
ion-molecule  reverse  reactions  become  more  important  as  entropy  effects  become  more  important: 

-  TAfS  (5) 


This  same  effect,  of  course,  controls  neutral  reactions  and  may  very  well  explain  why  the  rate 
of  soot  formation  peaks  in  the  temperature  range  of  1400  to  1900  K  (see,  e.g..  Ref.  23).  The  usual 
explanation  involves  a  competition  between  the  rate  of  soot  formation  and  the  rate  of  soot  oxidation. 
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The  strong  temperature  effect  on  ion  growth  reactions  has  shown  up  on  several  occasions. 
In  some  of  our  first  computer  runs  with  Frenklach  and  Wang,  a  table  with  the  experimental  flame 
temperature  vs.  distance,  in  which  one  point  was  in  error,  was  inadvertently  used  as  input.  This 
showed  up  in  the  calculation  of  ion  profiles  as  a  dip  in  the  concentration  profiles.  In  the  early  ion 
profile  calculations  made  by  Frenklach  and  Wang,  the  ions  always  peaked  much  closer  to  the  burner 
than  the  experimental  data,  occurring  close  to  the  burner  where  the  flame  temperature  is  low. 

Our  findings,  submitted  to  the  25th  Combustion  Symposium  and  presented  in  Appendbc  A, 
can  be  summarized  briefly.  First,  we  find  that  our  calculated  ion  concentrations  are  very  sensitive 
to  the  ion  Gibb's  free  energies  of  formation  used.  Second,  the  calculated  ion  concentrations  are 
almost  always  low  compared  with  experimental  concentrations.  Third,  electron-ion  recombination 
rates  can  be  lowered  by  a  factor  of  100  with  little  effect  on  the  profiles.  This  means  that  the  model 
predictions  would  be  relatively  unaffected  by  the  inclusion  of  other  negatively  charged  species,  which 
have  been  reported  in  sooting  flames,  and  could  result  in  higher  positive  ion  concentrations. 

We  ask  why  the  thermodynamics,  which  controls  the  reverse  reaction  rates,  plays  such  an 
important  role.  We  also  ask  why  the  calculated  concentrations  are  so  much  smaller  than  the 
measured  concentrations.  It  could,  of  course,  be  that  ion  growth  is  not  the  source  of  larger  ions  in 
flames,  but  we  have  eliminated  many  other  possible  mechanisms  [1,24].  Earlier  Homann  was  critical 
of  an  ion  growth  mechanism.  Homann  and  associates  [25]  now  seem  to  believe  that  the  source  of 
large  ions  is  ion  growth,  as  we  proposed  years  ago!  [26]  Another  possibility  is  that  the  experimental 
measurements  are  off  It  is  difficult  to  see  how  they  could  be  off  by  this  much.  Equivalent 
experimental  results  have  been  obtained  by  Homann  and  associates  in  Germany  [27,28]  and  by 
Delfau  and  associates  in  France  [29].  An  attempt  to  bring  the  model  into  agreement  with 
experimental  through  incorporation  of  the  modified  Langevin  theory  was  ineffective. 

Several  rather  "brutal"  modifications  were  made  in  an  attempt  to  determine  what  has  to  be 
done  to  obtain  agreement.  The  best  agreement  was  found  by  reducing  the  flame  temperature  profile 
by  10  %  at  the  maximum  temperature,  and  linearly  decreasing  the  extent  of  reduction  to  zero  at  the 
burner  surface.  The  results  for  modeling  the  profiles  of  ions  up  to  the  perinaphthenyl  ion,  C^H^HS^ 
in  Table  I,  looked  acceptable.  The  data  and  discussion  are  in  one  of  the  rejected  Symposium  papers. 
Appendix  A. 

Subsequent  to  submitting  the  Symposium  paper,  we  extended  the  model  to  the  protonated 
coronene  ion,  C24Hi3''’.  The  temperature  effect  became  even  more  evident.  If  a  reduced  temperature 
maximum  of  1700  K  is  assumed  and  the  temperature  profile  closer  to  the  burner  is  modified  to 
extrapolate  to  the  experimental  burner  surface  temperature,  then  the  ratios  of  the  maximum 
experimental  to  the  maximum  calculated  ion  concentrations  are  obtained  as  shown  in  Fig.  3. 
Negative  log  values  (left  of  zero)  imply  the  observed  concentrations  are  low  compared  with  the 
model  and  are  acceptable  because  reaction  rate  coefficients  involved  in  the  production  of  such 
species  can  be  decreased.  Log  ratios  to  the  right  of  zero  are  of  greater  concern.  In  addition  to 
reducing  the  maximum  temperature,  the  AHf(CiiH8H^)  was  reduced  by  about  3%.  The  apparent 
requirements  to  adjust  the  experimental  temperature  will  be  our  major  consideration  in  continuing 
this  work.  It  would  be  premature  to  abandon  the  concept  of  ion  growth  and  the  ionic  mechanism 
of  soot  formation  because  there  is  so  much  evidence  in  favor  of  both  concepts  [see  e.g,  references 
1,2,30  and  below].  It  is  also  difficult  to  explain  the  rather  good  agreement  between  experimental 
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and  calculated  ion  concentrations  just  by  scaling  the  temperature  profile,  even  by  the  large  amount 
required. 


VII.  COMPARISON  OF  THE  RATE  OF  CARBON  SPECIES  GROWTH 
FOR  THE  FREE  RADICAL  AND  IONIC  MECHANISMS 


In  the  final  report  on  the  last  contract  [31],  we  compared  the  rates  of  increase  in  carbon 
content  by  a  free  radical  and  the  ionic  mechanism,  using  experimentally  measured  neutral  and  ion 
species  and  the  appropriate  rate  coefficients.  That  eomparison  has  been  further  refined  and  has 
been  accepted  for  publication  [1].  The  more  relevant  experimental  neutral  species  data  of  Bittner 
and  Howard  [32]  are  used  rather  than  the  data  of  Vovelle  [33]  (our  burner  and  conditions  were 
designed  to  match  those  of  Bittner  and  Howard  so  such  comparisons  could  be  made).  The 
experimental  data  on  which  this  analysis  was  made  are  presented  in  Fig.  4.  To  the  extent  that  the 
Bittner  and  Howard  data  can  be  smoothly  extrapolated,  beyond  250  u  the  ion  and  neutral 
concentrations  become  equal  and  the  difference  in  growth  rate  will  be  determined  by:  the  reaction 
rate  coefficients  in  the  neutral  and  ion  mechanisms;  the  concentration  of  the  other  reactants 
(hydrogen  atoms  and  acetylene  for  the  free  radical  mechanism,  and  acetylene  for  the  ionic 
mechanism);  and  the  number  of  reaction  steps  to  account  for  the  same  number  of  earbon  atoms 
added,  ratio  2:1  for  the  free  radical  mechanism  over  the  ionic  mechanism.  On  this  basis,  the  ionic 
mechanism  is  faster;  compare  6.7  /rs  for  the  ionic  mechanism  with  25  /ts  for  the  free  radical 
mechanism.  This  is  still  not  a  conclusive  comparison,  recognizing  the  inaccuracies  in  the  data  and 
in  the  rate  coefficients,  but  it  certainly  suggests  that  the  commonly  promoted  radical  mechanism  is 
not  necessarily  the  dominant  path  to  soot  formation. 


VIII.  COMPARISON  OF  THE  THERMOCHEMICAL  DRIVING  FORCE 
FOR  THE  IONIC  AND  FREE  RADICAL  MECHANISMS 


The  ultimate  driving  force  for  any  sequence  of  chemical  reactions,  such  as  must  occur  for 
soot  formation,  is  thermochemistry.  As  water  flows  downhill,  reactions  proceed  in  the  direction  of 
a  negative  gradient  in  free  energy.  To  compare  this  driving  force  for  ions,  free  radicals,  and  neutral 
molecules,  we  eompare  the  free  energy  of  formation  of  these  three  components  of  the  ionic  and  free 
radical  mechanism  of  soot  formation  in  Figs.  5  and  6  for  the  respective  temperatures,  1000  and  2000 
K. 


Several  features  are  obvious  on  inspection.  First,  the  driving  foree  (the  negative  slope  of  the 
curve)  for  ions  is  distinctly  greater  than  for  free  radicals  or  neutral  molecules  in  the  initial  stages  of 
growth.  Second,  the  free  energy  per  carbon  atom  of  any  of  the  species  (free  radical,  neutral 
molecule,  or  ion)  rapidly  levels  off  as  the  molecules  become  larger.  Third,  the  ion  and  neutral 
species  free  energies  of  formation  converge  when  the  species  contain  about  30  carbon  atoms,  i.e., 
a  molecular  weight  of  about  350  to  400  u. 
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The  conclusion  to  be  drawn  is  dear.  There  is  very  little  thermochemical  driving  force  for  the 
growth  of  free  radicals  (or  neutral  molecules)  compared  to  the  corresponding  driving  force  for  ions. 


IX.  IONS  IN  BENZENE/OXYGEN  FLAMES 


The  above  kinetic  scheme  has  been  applied  to  modeling  acetylene/oxygen  flames.  In  a 
previous  program  [3],  we  measured  ion  spectra  of  benzene/oxygen  flames  in  the  same  flames  and  on 
the  same  burner  as  used  by  Bittner  and  Howard  [32]  to  measure  neutral  species  profiles.  This  was 
done  to  have  a  data  base  available  for  modeling  a  different  flame.  We  are  currently  analyzing  that 
data  in  preparation  for  modeling  the  benzene/oxygen  flame. 

The  ion  structures  of  these  flames  are  quite  different.  As  in  acetylene  flames,  is 

assumed  to  be  the  primary  ion,  yet  its  concentration  profile  is  different  in  benzene  flames.  In  both 
nonsooting  {<j>  =  1.8)  and  sooting  {<j)  =  2.0)  benzene/oxygen  flames,  this  ion  exhibits  two  peaks,  a 
small  one  near  the  burner  and  a  relatively  large  second  peak  downstream.  It  slowly  decays  further 
from  the  burner.  Another  unique  feature  of  the  benzene  flame  is  that  a  large  number  of  oxygenated 
ions,  e.g.,  are  observed  which  are  not  observed  in  acetylene/oxygen  flames.  Additionally, 

of  the  ions  observed  in  acetylene/oxygen  flames  are  also  observed  in  the  benzene/oxygen  flame. 

The  apparent  stability  of  the  CjHj^  ion  in  the  benzene  flame  contrasts  with  its  behavior  in 
the  CjHj/Oj  flame.  In  the  benzene  flame  its  concentration  is  an  order  of  magnitude  greater  than 
the  concentration  of  other  ion  species;  in  the  sooting  acetylene  flame  its  concentration  is  similar  to 
those  of  other  ions.  In  the  benzene  flames  its  concentration  decays  with  distance  more  slowly  than 
in  the  acetylene  flame.  The  difference  can  probably  be  explained  by  the  equilibrium  between  the 
linear  reactive  ion  and  the  stable  cyclic  nonreactive  ion  favoring  the  nonreactive  form  at  the  lower 
temperature  1900  K  (maximum)  in  the  benzene  flame  compared  to  2020  K  (maximum)  in  the 
acetylene  flame.  Some  of  the  differences  may  be  due  to  the  fact  that  the  sooting  acetylene  flame 
is  more  fuel  rich,  relative  to  the  critical  equivalence  ratio  for  soot  formation,  than  the  benzene 
flames.  There  are  other  great  differences  (see  Section  X  and  Appendix  B)  between  the  two  flames. 
Soot  formation  in  the  benzene  flame  is  greatly  in  excess  of  equilibrium  and  occurs  at  an  equivalence 
ratio  less  that  the  equilibrium  equivalence  ratio  for  soot  formation.  In  contrast,  soot  formation  in 
the  acetylene  flame  is  less  than  equilibrium  and  occurs  at  an  equivalence  ratio  greater  than  the 
equilibrium  equivalence  ratio.  Further,  since  the  total  ion  concentrations  in  the  benzene  flames  are 
greater  than  in  the  acetylene  flame,  compare  about  10“  with  7  x  10’  ions  cm the  greater  stability 
of  QHj'"  (cyclic)  compared  to  CaHj""  (linear)  apparently  does  not  limit  the  formation  of  large  ions, 
but  may  spread  their  formation  over  a  greater  distance.  These  differences  between  the  two  flames 
will  be  explained  by  successful  models. 

Recently  Loffler  and  Homann  [34,35]  reported  ion  profile  measurements  in  sooting 
benzene/oxygen  flames.  The  conditions  in  their  experiments  and  in  ours  are  compared  in  Table  II. 
In  Table  III  we  compare  the  ions  observed  in  acetylene/oxygen  and  benzene/oxygen  flames  at  both 
AeroChem  and  at  Darmstadt.  In  their  paper,  Loffler  and  Homann  concentrate  [34]  on  oxygenated 
ions,  but  observed  ion  mass  spectra  are  reported  in  Loffler's  thesis  [35].  We  have  therefore  used 
these  spectra  to  determine  other  ions  observed  in  their  flames. 
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In  any  such  comparison  of  mass  spectrometer  data,  one  should  be  aware  that  the  reporting 
of  the  appearance  of  an  ion  involves  a  somewhat  subjective  factor  in  the  choice  of  which  peaks  are 
identified  as  important,  and  an  objective  factor  that  is  dependent  upon  the  sensitivity  and  resolution 
of  the  mass  spectrometer.  The  list  used  here  for  acetylene/oxygen  flames  at  Darmstadt  was  kindly 
furnished  by  Homann  [36]  to  supplement  those  given  in  the  paper  [27].  Many  other  ions  were 
identified  by  Homann;  we  used  only  those  identified  as  "prominent  peaks".  Some  of  these  ions,  e.g., 
71,  83,  and  85  u  are  identified  by  Ldffler  and  Homann  in  the  acetylene  flame  as  an  oxygen  molecule 
attached  to  an  ion.  This  indicates  reactions  in  the  sampling  cone  of  their  instrument;  these  ions 
would  not  be  expected  to  be  stable  in  a  flame.  Our  experience  has  been  that  the  attachment  of 
water,  sometimes  more  than  one  water  molecule,  to  an  ion  indicates  sampling  problems.  Such  ion 
signals  can  be  eliminated  by  changing  the  sampling  cone  geometries,  temperature,  bias  voltage,  or 
the  pressure  behind  the  cone.  Some  of  the  oxygenated  species  observed  by  Gerhardt  and  Homann, 
may  in  fact  be  water  clustered  to  a  hydrocarbon  ion,  e.g.  in  Table  III,  for  ion  masses:  47,  67,  69,  73, 
81,  85,  95,  105,  109,  and  119.  This  could  also  be  the  case  for  some  of  the  oxygenated  species 
observed  in  both  laboratories  in  benzene  flames. 

The  above  consideration  raises  again  [37]  the  problem  of  sampling  ions  which  have  very  high 
reaction  rate  coefficients.  Are  the  identified  ions  the  true  flame  ions  or  ions  that  have  been 
produced  in  the  sampling  system,  by  e.g.,  reacting  with  a  species  in  large  concentration,  e.g.,  H2O, 
as  discussed  above,  Hj,  H*,  or  C2H2?  This  would  muddle  an  already  complicated  problem.  It  would 
mean  that  some  of  the  observed  ions  are  not  flame  ions,  but  are  produced  in  the  sampling  cone. 
To  model  the  flame  chemistry,  one  would  first  have  to  deduce,  from  the  observed  ions,  the  actual 
ions  in  the  flame. 

In  general,  for  the  benzene  flame,  Table  II,  the  agreement  in  observed  ions  between  Loffler 
and  Homann  and  AeroChem  are  consistent,  except  for  CijH,"^.  We  observe  this  ion  in  benzene 
flames  but  they  do  not.  This  is  a  very  stable  ion  and  we  both  observe  it  in  large  concentrations  in 
acetylene  flames.  Instead  of  this  ion,  they  observe  Ci2H90^,  4  u  greater  than  C^H,"",  an  ion  we 
observe  only  in  very  small  concentrations.  Because  we  had  previously  reported  QjH,""  in  benzene 
flames  [38],  Loffler  and  Homann  were  very  careful  to  confirm  the  fact  that  it  was  not  observed  [39]. 
Their  instrumentation  is  better  than  ours,  but  it  is  difficult  to  see  how  we  could  be  off  in  our 
calibration  by  4  u.  If  we  shift  our  mass  scale  by  4  u  to  force  a  match  at  169  u,  then  our  agreement 
with  identifying  larger  ions  falls  from  80  ions  to  54  ions.  Further,  we  would  identify  40  ions  which 
they  do  not  observe  and  are  not  observed  in  the  acetylene  flame. 

One  possibility  is  that  the  two  ions,  and  Ci2H,0'^,  have  similar  thermodynamic 

stabilities,  so  that  the  differences  in  our  two  experiments.  Table  II,  are  sufficient  to  switch  the 
stability  from  one  to  the  other  ion.  To  determine  the  feasibility  of  this  would  require  estimating  the 
thermodynamics  for  other  ions  of  mass  169  u,  which  we  will  do.  This  ion  may  derive  from  an  ion 
of  mass  152  u  with  H2O  replacing  a  H  atom  or  151  u  with  H2O  added.  The  question  remains,  why 
does  Darmstadt  not  observe  mass  165?  The  differences  in  the  two  flames,  Table  II,  do  not  explain 
this  observation  difference.  Hopefully,  the  difference  will  become  clear  when  we  model  this  flame. 

The  next  step  in  analyzing  the  benzene/oxygen  flame  will  be  to  calibrate  our  mass 
spectrometer  comparing  the  mass  spectral  data  with  our  Langmuir  probe  data,  as  we  did  for  the 
acetylene  flame  [40]. 
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X.  THERMODYNAMIC  CONSIDERATIONS  OF  SOOT  FORMATION 


Data  in  the  literature  on  soot  and  acetylene  concentration  variation  through  the  threshold 
equivalence  ratio  for  soot  formation  were  compared  with  thermodynamic  equilibrium  concentrations. 
One  of  the  reasons  for  exploring  this  is  that  at  soot  threshold,  the  only  other  observed  property  of 
the  flame  which  changes  dramatically  with  equivalence  ratio  is  the  ion  concentration  and  size  of  the 
ions,  Fig.  7. 

It  was  found  that  in  most  flames  soot  is  observed  at  C/0  ratios  very  much  lower  than  where 
soot  is  predicted  to  form  by  equilibrium,  and  the  concentrations  of  soot,  acetylene  and  other 
hydrocarbons  are  in  great  excess  of  equilibrium.  However  in  the  well  studied  acetylene/oxygen  flame, 
the  "standard  flame",  the  opposite  is  true.  For  the  acetylene/oxygen  flame,  acetylene  is  in  excess  of 
equilibrium  concentrations,  and  soot  concentrations  are  less  than  equilibrium.  To  our  knowledge, 
this  is  the  first  time  this  has  been  noted.  This  work  is  presented  in  one  of  the  papers  rejected  for 
the  Combustion  Symposium,  see  Appendix  B.  We  will  expand  the  manuscript,  explaining  the 
implications  for  the  ionic  mechanism  of  soot  formation,  and  submit  it  for  publication  in  Combustion 
and  Flame. 
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Calcote,  H.F.,  "Modeling  Study  to  Evaluate  the  Ionic  Mechanism  of  Soot  Formation,"  presented  at 
AFOSR  Contractors  Meeting  on  Propulsion  Research,  La  Jolla,  CA,  15-19  June  1992. 

Calcote,  H.F.  and  Gill,  R.J.,  "A  Detailed  Computer  Model  of  Ion  Growth  in  a  Fuel  Rich 
Acetylene/Oxygen  Flame,"  presented  at  Joint  Meeting  of  the  British  and  German  Sections  of  the 
Combustion  Institute,  Queens  College,  Cambridge,  England,  29  March  to  2  April  1993. 

Calcote,  H.F.  and  Gill,  R.J.,  "A  Model  of  Ion  Growth  in  a  Fuel  Rich  Acetylene/Oxygen  Flame," 
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Appendix  C. 
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Mechanisms  and  Models  of  Soot  Formation,  Ruprecht-Karls-Universitat,  Heidelberg,  Germany,  29 
September  -  2  October  1991. 

Calcote,  H.F.  and  Gill,  R.G.,  "Comparison  of  the  Ionic  Mechanism  of  Soot  Formation  with  a  Free 
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TABLE  I 

REACTIONS  AND  RATE  COEFFICIENTS  (LANGEVIN  AND  EXTENDED  LANGEVIN)  FOR 

IONIC  MECHANISM  OF  SOOT  FORMATION 

kCT)  =  A  T"  exp  (-E^^t/RT) 
units  J,  cm^  mole,  s 


Langevin  Rate  Extended  Langevin 

Parameters  Rate  Parameters 


REACTIONS 

A 

n 

Eact 

A 

n 

^act 

+ 

C2H2 

=  C5H2H+ 

-1-  H2 

6.50e+14 

0 

0 

8.60e+13 

0.3 

0 

+ 

C2H2 

=  CsHs^ 

6.50e+14 

0 

0 

8.60e+13 

0.3 

0 

CaHa"' 

+ 

C4H2 

=  C5H2H+ 

+  C2H2 

7.40e+14 

0 

0 

7.28e+13 

0.3 

0 

CaHa"' 

+ 

C4H2 

=  C7H4H+ 

7.40e+14 

0 

0 

7.28e+13 

0.3 

0 

C4H5'" 

+ 

C2H2 

=  HsCg-" 

+  H2 

6.20e+14 

0 

0 

9,17e+13 

0.3 

0 

C5H2H+ 

+ 

H2 

=  CsHs-" 

9.00e+14 

0 

0 

8.97e+14 

0.0 

0 

C5H2H+ 

+ 

C2H2 

=  C7H4H+ 

6.00e+14 

0 

0 

9.51e+13 

0.3 

0 

CsHsH-" 

+ 

CaH4 

=  H5C6+ 

-P  C2H2 

6.50e+14 

0 

0 

9.51e+13 

0.3 

0 

+ 

C2H2 

=  C7H4H+ 

•+  H2 

6.00e+14 

0 

0 

9.58e  +  13 

0.3 

0 

CsHs^ 

+ 

C2H2 

=  C7H5+ 

4-  H2 

2.40e+14 

0 

-319 

9.58e+13 

0.3 

0 

CsHs^ 

C4H2 

=  C9H7+ 

6.50e+14 

0 

0 

9.58e+13 

0.3 

0 

CsHs-^ 

+ 

C4H2 

=  H7C9+ 

6.50e+14 

0 

0 

8.09e+13 

0,3 

0 

CsHs^ 

+ 

C2H2 

=  C7H7+ 

6.50e+14 

0 

0 

8.09e4-13 

0.3 

0 

+ 

C2H2 

=  C3H7+ 

5.80e+14 

0 

0 

9.94e+13 

0.3 

0 

C7H4H+ 

+ 

H2 

=  C7H7+ 

8.90e+14 

0 

0 

8.93e+14 

0.0 

0 

C7H4H+ 

+ 

H2 

=  H7C7^ 

8.90e+14 

0 

0 

8.93e+14 

0,0 

0 

CyHs" 

+ 

C2H2 

“  *^9^7 

5.70e+14 

0 

0 

1.03e+14 

0.3 

0 

+ 

C2H2 

=  H7C9+ 

5.70e  +  14 

0 

0 

1.03e+14 

0.3 

0 

+ 

M 

=  C7H4H  + 

-P  M 

4.60e  +  14 

0 

0 

1.03e+14 

0.3 

0 

C7H5+ 

+ 

CaH4 

=  C3H7+ 

+  C2H2 

6.l5e+14 

0 

0 

1.03e+14 

0.3 

0 

+ 

CaH4 

=  C,oHg^ 

6.15e  +  14 

0 

0 

1.03e+14 

0.3 

0 

C7H4H+ 

+ 

C2H2 

=  H7C9+ 

5.70e+14 

0 

0 

1.03e  +  14 

0.3 

0 
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Langevin  Rate  Extended  Langevin 

Parameters  Rate  Parameters 


REACTIONS 

A 

n 

^act 

A 

n 

^act 

+ 

C2H2 

=  C5H2H+ 

+ 

H2 

6.50e+14 

0 

0 

8.60e+13 

0.3 

0 

C7H4H+ 

+ 

C2H2 

= 

5.70e+14 

0 

0 

1.03e+14 

0.3 

0 

+ 

M 

=  H7C7^ 

+ 

M 

4.60e+14 

0 

0 

1.03e+14 

0.3 

0 

+ 

C2H2 

=  HyCg-" 

+ 

H2 

5.70e+14 

0 

0 

1.03e+14 

0.3 

0 

H7C/ 

+ 

C2H2 

=  H7C9^ 

+ 

H2 

5.70e+14 

0 

0 

1.03e+14 

0.3 

0 

C7H7+ 

+ 

C4H2 

=  CiiHgH+ 

6.10e+14 

0 

0 

8.69e+13 

0.3 

0 

+ 

C4H2 

=  C^HaH^ 

6.10e+14 

0 

0 

8.69e+13 

0.3 

0 

+ 

C2H2 

=  CioHg^ 

5.64e+14 

0 

0 

1.06e+14 

0.3 

0 

CsHy-^ 

+ 

C2H2 

=  HgC^o^ 

5.64e+14 

0 

0 

1.06e+14 

0.3 

0 

CsHr'" 

+ 

C3H4 

= 

+ 

H2 

6.02e+14 

0 

0 

1.06e+14 

0.3 

0 

CaHy-^ 

+ 

C4H2 

= 

5.95e+ 14 

0 

0 

8.94e  +  13 

0.3 

0 

CgHr^ 

+ 

M 

=  H7C9-^ 

+ 

M 

4.50e+14 

0 

0 

1.09e  +  14 

0.3 

0 

CgHr'^ 

+ 

C2H2 

=  C^iHaH^ 

5.60e+14 

0 

0 

1,09e+14 

0.3 

0 

+ 

C2H2 

=  C^iHaH^ 

5.60e+14 

0 

0 

1.09e  +  14 

0.3 

0 

CgHT-^ 

+ 

C4H2 

=  C^gHgHa'*’ 

5.80e+14 

0 

0 

9.16e+13 

0.3 

0 

H7V 

+ 

C4H2 

=  Ci3HgH3+ 

5.80e+14 

0 

0 

9.16e+13 

0.3 

0 

HgCio^ 

+ 

M 

=  CioHg-^ 

+ 

M 

4.50e+14 

0 

0 

1.12e+14 

0.3 

0 

CioHg^ 

+ 

C2H2 

=  C12H9+ 

+ 

Hg 

5.53e+14 

0 

0 

1.12e+14 

0.3 

0 

HgCio^ 

+ 

C2H2 

=  Ci2H9^ 

+ 

Hg 

5.53e+14 

0 

0 

1.12e+14 

0.3 

0 

HgCio^ 

+ 

C2H2 

=  H9C,2-" 

+ 

Hg 

5.53e+14 

0 

0 

1.12e+14 

0.3 

0 

CioHg^ 

+ 

C3H4 

=  C^angHa"*" 

+ 

2H2 

5.85e+14 

0 

0 

1.12e+14 

0.3 

0 

CioHg^ 

+ 

C4H2 

= 

+ 

CgHg 

5.75e+14 

0 

0 

9.40e+13 

0.3 

0 

HgCio^ 

+ 

C4H2 

= 

+ 

CgHg 

5.75e+14 

0 

0 

9.40e+13 

0.3 

0 

HgCio^ 

+ 

C4H2 

=  H9C,2-^ 

+ 

CgHg 

5.75e+14 

0 

0 

9.40e  +  13 

0.3 

0 

C„HaH^ 

+ 

C2H2 

=  C^aNgHg'^ 

+ 

Hg 

5.50e+14 

0 

0 

1.14e+14 

0.3 

0 

CnHg^ 

+ 

C2H2 

=  CisHaHa-^ 

+ 

Hg 

5.50e+14 

0 

0 

1.14e+14 

0.3 

0 

HgCn^ 

+ 

M 

=  C11H3H+ 

+ 

M 

4.50e+14 

0 

0 

1,14e+14 

0.3 

0 
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Langevin  Rate  Extended  Langevin 

Parameters  Rate  Parameters 


REACTIONS  A  n  A  n  E^^, 


+ 

^2^2 

=  C5H2H+ 

-1-  H2 

6.50e+14 

0 

0 

8.60e+13 

0.3 

0 

+ 

M 

=  C„H9^ 

+  M 

4.50e  +  14 

0 

0 

1.14e+14 

0.3 

0 

+ 

M 

= 

+  M 

4.50e+14 

0 

0 

1.14e+14 

0.3 

0 

+ 

M 

= 

+  M 

4.4764-14 

0 

0 

1.16e+14 

0.3 

0 

+ 

E 

-C2H2 

-1-  CH 

1.05e+19 

-0.5 

0 

1.05e+19 

-0.5 

0 

C4H5+ 

+ 

E 

C2H2 

^2^3 

1.2064-19 

>0.5 

0 

1.20e+19 

-0.5 

0 

C5H2H+ 

+ 

E 

-C3H3 

+  C2 

1.3364-19 

-0.5 

0 

1.336+19 

-0.5 

0 

CsHs^ 

+ 

E 

-  C3H3 

4"  ^2^2 

1.406+19 

-0.5 

0 

1.406+19 

-0.5 

0 

HsCe'" 

E 

+  C^H 

1.426+19 

-0.5 

0 

1.42e+19 

-0.5 

0 

C7H4H+ 

+ 

E 

C4H2 

+ 

1.50e+19 

-0.5 

0 

1.506+19 

-0.5 

0 

H7C/ 

+ 

E 

+  CH3 

1.526+19 

-0.5 

0 

1.52e+19 

-0.5 

0 

CsHy-^ 

+ 

E 

-  CgHg 

+  C2H 

1.61e  +  19 

-0.5 

0 

1.61e+19 

-0.5 

0 

H7C9+ 

+ 

E 

-CgHg 

+  CH 

1.706  +  19 

-0.5 

0 

1.706+19 

-0.5 

0 

+ 

E 

-(-  H 

1.796+19 

-0.5 

0 

1.796+19 

-0.5 

0 

CnHsH^ 

+ 

E 

^10*^8 

+  CH 

1.86e+19 

-0.5 

0 

1.866+19 

-0.5 

0 

CisHg'' 

+ 

E 

-*  C12H3 

+  H 

1.936  +  19 

-0.5 

0 

1.936+19 

-0.5 

0 

+ 

E 

1.99e+19 

-0.5 

0 

1.996+19 

-0.5 

0 

18 


TP-531 


TABLE  II 


COMPARISON  OF  BENZENE/OXYGEN  FLAMES  OBSERVED 
AT  AEROCHEM  AND  AT  DARMSTADT 


AeroChem 

Darmstadt 

EQUIVALENCE  RATIO 

1.8,  2.0 

1.9,  2.0 

-1-30  mol%  argon 

UNBURNED  GAS  VELOCITY,  cm 

50 

42 

PRESSURE,  kPa 

2.67 

2.67 

THRESHOLD  SOOT  FORMATION,  cp 

1.90 

1.83 

TEMPERATURE  PEAK,  K: 

1900 

2094 

-  at  mm  above  burner 

12.9 

11.6 

-  for  equivalence  ratio 

1.80 

1.95 

19 
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TABLE  III 


COMPARISON  OF  IONS  OBSERVED  IN  ACETYLENE/OXYGEN  AND 
BENZENE/OXYGEN  FLAMES  AND  IONS  OBSERVED  AT  AEROCHEM 
AND  AT  HOMANN'S  LABORATORY  IN  DARMSTADT,  GERMANY 


ACETYLENE/OXYGEN  BENZENE/OXYGEN 


ION  MASS 

ION  FORMULA 

AERO® 

DARMS^ 

AERO'" 

□ARMS' 

19 

H3O+ 

X 

24 

C 

X 

X 

31 

CHgO^ 

X 

33 

CH5O  + 

X 

39 

X 

X 

X 

X 

41 

CgHO 

X 

X 

43 

C2H3O+ 

X 

X 

X 

47 

C2H7O+ 

X 

49 

C4H  + 

X 

51 

C4H3+ 

X 

X 

X 

53 

X 

X 

X 

X 

55 

C3H3O+ 

X 

57 

C4H9'^,  C3H5O+ 

X 

X 

63 

X 

X 

X 

X 

65 

X 

X 

67 

C4H30^ 

X 

X 

X 

69 

C4H5O+ 

X 

71 

C3H3+-02 

X 

73 

CeH  +  ,  C4HgO'^ 

X 

X 

75 

CgHa"' 

X 

X 

X 

77 

CeHs'*' 

X 

X 

X 

20 
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ACETYLENE/OXYGEN  BENZENE/OXYGEN 


ION  MASS 

ION  FORMULA. 

AERO® 

DARMS*^ 

AERO‘S 

DARMS' 

79 

CeH?* 

X 

X 

X 

81 

CgHgO^ 

X 

X 

X 

83 

C4H3+-02 

X 

85 

C4H5+-02 

X 

87 

X 

X 

X 

89 

X 

X 

X 

91 

C7H7+ 

X 

X 

X 

X 

95 

C6H7O+.  C5H3'^-02 

X 

X 

X 

97 

CjH+.CsHs^-Oj 

X 

X 

99 

X 

X 

X 

103 

X 

X 

X 

X 

105 

C7H50+ 

X 

X 

X 

109 

C7H90+ 

X 

X 

X 

111 

C9H3^CeH7+.02 

X 

X 

X 

115 

X 

X 

X 

X 

117 

CgHg'' 

X 

X 

119 

CgH70+ 

X 

X 

X 

121 

C3H9O+ 

X 

X 

123 

CgH^.O-^ 

X 

X 

X 

124 

X 

X 

126 

CioHe^ 

X 

X 

129 

X 

X 

X 

X 

131 

C9H70^C^oH11^ 

X 

X 

X 

133 

C9HgO^ 

X 

X 

135 

CnHa-^ 

X 

139 

X 

X 

141 

CiiH9'^ 

X 

X 

X 

X 

21 
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ACETYLENE/OXYGEN  BENZENE/OXYGEN 


ION  MASS 

ION  FORMULj4 

AERO^ 

DARMS^ 

AERO‘S 

DARMS' 

143 

C,oH,0+ 

X 

X 

X 

145 

X 

X 

X 

147 

X 

152 

X 

X 

153 

Ci2H9^ 

X 

X 

X 

155 

X 

X 

157 

C^HgO^ 

X 

X 

159 

^13*^3^ 

X 

165 

X 

X 

X 

169 

0^2^90'*' 

X 

170 

CisHioO^ 

X 

X 

171 

X 

177 

Ci4H9^ 

X 

X 

X 

179 

X 

X 

181 

CiaHgO^ 

X 

X 

183 

X 

189 

CisHg^ 

X 

X 

X 

191 

CisHu^ 

X 

X 

X 

192 

CuHsO" 

X 

X 

193 

X 

X 

203 

X 

X 

X 

204 

C^sHgO^ 

X 

X 

205 

Ci5H90^ 

X 

X 

215 

X 

X 

X 

X 

217 

C16H9O+ 

X 

X 

219 

X 

X 

227 

X 

X 

X 

X 

22 
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ACETYLENE/OXYGEN  BENZENE/OXYGEN 


ION  MASS 

ION  FORMULA^ 

AERO® 

DARMS^’ 

AERO° 

DARMS' 

228 

Ci/HgO-^ 

X 

X 

231 

X 

X 

239 

^19^11 

X 

X 

X 

X 

240 

C,8H50+ 

X 

X 

241 

X 

243 

C18H11O  + 

X 

X 

250 

’^20*^10''' 

X 

251 

^20*^11"*^ 

X 

X 

X 

253 

X 

X 

255 

X 

263 

X 

X 

X 

X 

265 

C20H9O  .C2iH^3 

X 

X 

267 

C2oHi,0*.C2iHi5'' 

X 

276 

C22H12’'' 

X 

277 

C22H,3^ 

X 

X 

X 

279 

^21*^11*^^'  ^22*^15 

X 

X 

287 

C23H1/ 

X 

289 

*^23*^13^ 

X 

X 

X 

293 

*^22^13*^^ 

X 

X 

300 

C24H,2^ 

X 

301 

C24Hi3^ 

X 

X 

X 

X 

303 

C23H,,0+ 

X 

X 

313 

*^25*^13^ 

X 

X 

X 

317 

^24^13*^^ 

X 

324 

C26Hi2^ 

X 

325 

^26*^1 3^ 

X 

X 

X 

329 

^25*^13*^ 

X 

23 
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ACETYLENE/OXYGEN  BENZENE/OXYGEN 


ION  MASS 

ION  FORMULA 

AERO® 

DARMS^ 

AERO® 

□ARMS' 

337 

^27^13^ 

X 

X 

X 

341 

X 

348 

^28^12^ 

X 

349 

^28*^13'*' 

X 

350 

^28^14^ 

X 

351 

^28^15^ 

X 

X 

X 

353 

C2,H^30  + 

X 

361 

^29^13*^ 

X 

X 

X 

X 

363 

C28HiiO^ 

X 

X 

X 

365 

*^28^13^ 

X 

X 

374 

C30Hi4^ 

X 

375 

^30^15"^ 

X 

X 

377 

C29H,30^ 

X 

X 

385 

C3iHi3^ 

X 

386 

C3iHi4^ 

X 

X 

387 

C3iHi5^ 

X 

X 

391 

^30^150^ 

X 

398 

^32^^14'^ 

X 

399 

^32^15^ 

X 

X 

X 

403 

C3iHi50^ 

X 

411 

^33^15^ 

X 

X 

X 

X 

415 

C32H,50  + 

X 

422 

C34Hi4-^ 

X 

423 

C34Hi5^ 

X 

X 

X 

X 

427 

C33Hi50^ 

X 

435 

C35Hi5^ 

X 

X 

X 

X 

439 

C34Hi50^ 

X 

24 


TP-531 


ACETYLENE/OXYGEN  BENZENE/OXYGEN 


ION  MASS 

ION  FORMULA 

AERO^ 

DARMS^^ 

AERO° 

□ARMS' 

446 

^36^14^ 

X 

447 

X 

X 

X 

X 

459 

^37*^15^ 

X 

X 

X 

X 

470 

X 

CO 

CO 

o 

X 

471 

X 

X 

X 

472 

*^38*^16'*' 

X 

473 

^38^17'*' 

X 

X 

483 

X 

X 

X 

X 

495 

^40*^15^ 

X 

X 

497 

*^40*^17^ 

X 

506 

C4iHi4^ 

X 

X 

507 

C4iHi5'' 

X 

508 

C4iHi6^ 

X 

X 

509 

C4iHi7^ 

X 

X 

X 

520 

^42*^16^ 

X 

521 

C42H17'' 

X 

X 

X 

533 

^43^17'' 

X 

X 

X 

X 

544 

+ 

CO 

X^ 

O 

X 

545 

C44Hi7^ 

X 

X 

X 

X 

557 

^45^17^ 

X 

X 

X 

X 

568 

^46^16'' 

X 

569 

^46^17'*' 

X 

581 

^47^17^ 

X 

®AeroChem,  Ref.  38,  40. 
^Darmstadt,  Refs.  27,  28,  39. 
'^Refs.  4,  41 . 

‘^Refs.  34,  35. 
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93- 


10^  10^  10^  10'^  1 

MASS,  u 

FIGURE  1  PAH  ION  DIAMETERS 


PAH  ion  diameters  calculated  from  experimental  ion  mobility  (Langevin), 
and  spheres  of  density  1.5  cm^  (p  =  1.5)  and  2.0  cm^  (p=2.0). 
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FIGURE  2  LOGio  OF  EQUILIBRIUM  CONSTANT  FOR  SYSTEM 
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C3H3  + 
C4H5-I- 
C5H2H-I- 
H5C6  + 
C7H4H-I- 
C7H7-1- 
C8H7-H 
C9H7-h 
C10H9-f- 
dlHSH-h 
C12H9+ 
C13H6H3-H 
C14H7H4+ 
C15H10H-I- 
C16H11  -h 
H11C17-1- 
C18H10H-t- 
C19H11 + 
H11C20  + 
C21H11  + 
C22H13+ 
C23H13-I- 
C24H13  + 


-3-2-10  1  2 


LOG  10  (EXP^.x/CALC  max) 

FIGURE  3  COMPARISON  OF  EXPERIMENTAL  AND  CALCULATED 

(CALC)(naxin,uni  CATION  CONCENTRATIONS. 


90-47B 


FIGURE  4  RATE  OF  DECAY  OF  NEUTRAL  AND  IONIC  SPECIES 
CONCENTRATIONS  IN  ACETYLENE/OXYGEN  FLAMES  AT  2.67  kPa  AND 
AN  EQUIVALENCE  RATIO  OF  3.0 
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FIGURE  5  FREE  ENERGY  OF  FORMATION  PER  CARBON  ATOM  FOR  IONS, 
FREE  RADICALS  AND  MOLECULES  AT  1000  K. 


FIGURE  6  FREE  ENERGY  OF  FORMATION  PER  CARBON  ATOM  FOR  IONS, 
FREE  RADICALS  AND  MOLECULES  AT  2000  K. 
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c/o 

0.6  0.8  1.0  1.2  1.4  1.6 


1.5  2.0  2.5  3.0  3.5  4.0 

EQUIVALENCE  RATIO  • 


ION  PROFILES  THROUGH  AN 
ACETYLENE-OXYGEN  FLAME  AT  2.7  kPa 


NEUTRAL  PROFILES  THROUGH  A  PROPANE- 
OXYGEN  FLAME  at  15  kPa 


from  Bockhorn.  Petting  and  Wenz,  Ber.  Bunsenges. 
Phys.  Chem.  87:1067(1983). 


FIGURE  7  COMPARISON  OF  CHANGES  IN  CONCENTRATIONS  OF  IONS  AND 
NEUTRALS  FOR  TWO  FLAMES  NEAR  SOOT  THRESHOLD. 
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APPENDIX  A 


A  DETAILED  MODEL  OF  CATION  GROWTH  IN  FLAMES* 


H.  F.  Calcote,  R.  J.  Gill 
AeroChem  Research  Laboratories,  Inc. 
Princeton,  NJ  08542-0012,  USA 

and 


F.  N.  Egolfopoulos 

Department  of  Mechanical  Engineering 
University  of  Southern  California 
Los  Angeles,  CA  90089-1453,  USA 


ABSTRACT 

A  detailed  model  of  ion  growth  in  the  standard  acetylene-oxygen  premixed  flat  flame  (equivalence 
ratio  =  3.0,  P  =  2.67  kPa,  u  =  50  cm/s)  was  developed  to  model  cation  growth  from  CjHj"'’  (mass 
39  u)  to  CjgHg''’  (165  u).  The  model  included  all  observed  ions  and  no  ions  were  included  which 
have  not  been  measured.  The  reaction  scheme  included  ion-molecule  association  reactions  with 
neutral  reactants  (ie,  H2,  C2^2>  ^4^2  C3H4),  ion-electron  dissociative  recombination, 

isomerization  reactions  and  ambipolar  diffusion.  The  Sandia  Flame  Code  was  modified  to 
accommodate  ionic  reactions  and  non-Arrhenius  rate  coefficients  for  ion-molecule  reactions.  The 
calculated  ion  profiles  were  compared  with  the  experimental  ion  profiles  measured  for  the  standard 
flame.  In  general,  the  calculated  ion  profiles  were  in  good  agreement  with  both  the  shape  and  the 
number  density  when  compared  to  the  experimental  measurements. 


INTRODUCTION 

Very  large  ions  have  been  observed  in  flames  [1-4]  and  these  have  been  associated  with  soot 
formation  [5-6].  The  source  of  these  large  ions  has  been  debated  by  the  combustion  community  and 
many  schemes  have  been  invented  to  explain  their  origin.  These  schemes  range  from  thermal 
ionization  of  large  polycyclic  aromatic  hydrocarbons  [7]  to  thermal  ionization  of  soot  particles  [8] 
and  a  kind  of  chemiionization  of  large  PCAH  [9,10].  All  of  these  require  that  the  charge  be 
transferred  from  a  larger  moiety  to  a  small  moiety  by  charge  transfer  or  ion-molecule  reactions.  This 
is  contrary  to  thermodynamic  expectations  because  the  larger  moiety  is  more  stable  than  the  smaller 
[11-14]  so  the  reactions  would  be  expected  to  progress  from  the  smaller  to  the  larger  moiety.  We 
have  proposed  that  the  large  ions  are  derived  from  smaller  ions  by  a  series  of  ion-molecule  reactions 
and  that  this  growth  of  ions  is  an  important  mechanism  in  incipient  soot  formation  [15,16].  The 


‘submitted  to  25th  International  Symposium  on  Combustion,  University  of  California,  Irvine, 
CA,  July  31  -  August  5,  1994,  November  1993. 
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initial  ion  is  produced  by  chemiionization  which  is  recognized  as  the  source  of  highly  non-equilibrium 
concentrations  of  small  ions  eg  HCO'*’  and  C3H3"*’  [6]  and  the  C3H3'''  ion  grows  by  rapid  ion- 
molecule  reactions  to  produce  larger  ions.  Homann,  in  a  recent  paper  [10]  seems  to  have  accepted 
the  growth  of  ions  as  responsible  for  the  large  ions  in  flames. 

The  challenge  was  to  develop  a  detailed  quantitative  mechanism  which  would  account  for  the 
observed  ions.  This  was  not  a  straight  forward  problem  because  the  necessary  data  bases  were 
unavailable.  Thus,  for  large  ions  at  flame  temperatures,  the  thermochemistry,  ion-molecule  reaction 
rate  coefficients,  ion-electron  recombination  coefficients,  and  ion  and  electron  transport  properties 
all  had  to  be  developed  [11,17,18].  An  additional  complication  was  that  a  computer  code  to  handle 
ions  and  to  do  some  of  the  other  tasks  required  to  model  an  ionic  mechanism  was  unavailable. 
Eraslan  and  Brown  [19]  had  developed  a  code  to  handle  ions  but  it  lacked  some  of  the  capabilities 
which  we  needed  for  this  study.  The  Sandia  Flame  Code  was  thus  modified  to  include  ionic 
reactions  and  to  satisfy  other  requirements.  When  these  components  were  developed,  a  reaction 
scheme  was  assembled  to  simulate  ion  growth  from  small  chemiions  to  perinaphthenyl,  Cl3H6H3^ 
mass  165  u  in  the  standard  acetylene/oxygen  flame  (</)=3.0,  p=2.67  kPa,  unburnt  flow  velocity=50 
cm/s). 


Clearly,  where  so  much  had  to  be  done  to  compare  experiment  and  theory,  there  remains 
much  more  to  be  done  in  terms  of  refinement  (eg,  Ref.  20).  The  objective  of  this  work  was  to 
determine  if  it  is  feasible  to  model  ion  growth  with  a  reasonable  set  of  assumptions.  If  it  is,  then 
those  assumptions  will  be  refined  to  increase  their  scientific  precision.  If  it  is  not  possible  to  account 
for  the  large  ions  by  the  proposed  mechanism,  the  challenge  becomes  greater  because  other 
mechanisms  for  large  ion  formation  in  flames  have  been  rejected  [6]. 

In  this  paper,  the  development  of  the  necessary  data  bases,  development  of  a  detailed 
mechanism,  modifications  to  the  Sandia  Flame  Code,  and  a  comparison  of  the  mechanism  with 
experimental  data  are  briefly  summarized.  The  details  will  be  reported  elsewhere  [11,17,18]. 


THERMODYNAMICS 

The  major  source  of  uncertainty  in  the  kinetics  model  was  the  thermodynamics  of  cations. 
These  uncertainties  were  manifest  in  the  flame  simulations  where  rates  for  reverse  reactions  were 
calculated  from  thermodynamics  and  forward  reaction  rate  coefficients.  The  required 
thermodynamic  data  for  the  ions  were  either  collected  from  the  literature  or  calculated  [11]. 
Available  compilations  of  data  for  ions,  e.g.  Refs.  13,  21,  22,  contained  very  few  odd  number  carbon 
atom  ions;  experimentally  these  species  dominate  in  the  flame.  The  compilations  included  only  a  few 
isomers  of  which  there  can  be  many  for  each  ion  mass.  Isomerization  was  included  in  the  model  as 
an  elementary  reaction  step. 

There  were  several  neutral  molecules  for  which  the  thermodynamic  data  are  uncertain;  for 
example,  diacetylene.  The  available  heats  of  formation  for  diacetylene  at  298  K  range  from  439  to 
473  Id/mol.  [24-28]  We  used  440  kJ/mol  [24]  because  this  value  has  been  used  previously  in  several 
thermochemical  calculations  with  large  ions.  For  a  few  cases,  such  as  for  the  cyclopropenyl  ion, 
there  was  sufficient  experimental  information  to  employ  the  more  accurate  statistical  mechanical 
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methods  to  calculate  Cp°  and  S°.  It  is  difficult  to  assess  the  accuracy  of  the  thermodynamic  data  for 
ions  (or  neutrals!)  because  experimental  data  on  heat  capacities  and  entropies,  or  equilibrium 
measurements,  at  flame  temperatures  are  unavailable.  The  estimated  errors  in  these  quantities  can 
be  cis  large  as  10%  at  1500  K  [29], 


REACTION  RATE  COEFFICIENTS 

Experimental  rate  coefficients  for  ion-molecule  and  ion-electron  recombination  reactions 
were  available  only  for  small  ions  near  room  temperature.  It  was  thus  necessary  to  estimate  most  of 
these  coefficients. 

A.  ION-MOLECULE  REACTIONS 

In  general,  experimental  reaction  rate  coefficients  are  very  close  to  the  rate  calculated  by 
Langevin  theory  [30]  which  does  not  have  a  temperature  coefficient  for  non-dipole  molecules. 
Langevin  theory  accounts  only  for  ion-molecule  collisions  and  assumes  all  collisions  lead  to  reaction. 
Langevin  theory  is  substantiated  by  it's  agreement  with  experiment.  In  the  Langevin  theory,  the  ion 
is  treated  as  a  point  charge.  This  is  certainly  not  true  for  large  ions;  even  when  the  charge  is 
localized,  an  approaching  reactant  will  be  shielded  by  the  rest  of  the  ion.  We  have  thus  extended  the 
Langevin  theory  to  include  large  ions,  and  this  leads  to  a  small  positive  temperature  coefficient  [18]. 
Inclusion  of  the  extended  Langevin  theory  made  a  negligible  difference  in  the  calculated  ion  profiles. 

Ion-molecule  reactions  in  our  reaction  scheme  were  written  toward  increasing  molecular  size; 
thus  as  the  temperature  is  increases  the  free  energy  of  reaction  becomes  more  positive,  due  to 
entropy  effects.  To  avoid  excessively  fast  reverse  reactions,  forward  reaction  rate  coefficients  were 
adjusted  so  that  the  reverse  reaction  rate  did  not  exceed  the  extended  Langevin  rate.  One  of  the 
major  problems  in  modeling  large  ion  formation  was  structure  identification;  mass  spectrometry  gives 
mass  only.  The  number  of  carbon  and  hydrogen  atoms  was  determined  by  the  use  of  isotopes  [4] 
so  the  molecular  formulas  are  well  known.  However,  for  a  given  molecular  formula  there  can  be 
several  isomeric  structures.  Each  of  these  was  considered  as  a  separate  reactant/product  entity.  To 
avoid  an  overly  complex  mechanism,  only  the  most  stable  structures  were  considered.  To  account  for 
the  fact  that  a  Langevin  ion  complex  may  lead  to  several  reaction  products  it  was  assumed,  as  a  first 
order  correction,  that  product  branching  was  determined  by  the  reaction  free  energy  of  that  path. 
The  set  of  64  ion-molecule  reactions  used  in  this  paper  were  chosen  from  a  set  of  192  reactions 
leading  to  q3H,H3^ 

Typical  reactions  with  their  rate  coefficients  are: 

HoCin"^  +  CoHo  HqC.o'*'  +  ^7  k  =  6.7  x  10^^  x  cm^/mol/s 

CgHi^  +  C2H2  -  C11H8H+  k  =  1.1  X  10^^  X  tO-32  cm^/mol/s 

C1JH9+  -I-  M  ^  -1-  M  k  =  1.1  X  10^"^  X  T^'^^  cm^/mol/s 

The  order  of  "C"  and  "H"  in  the  ions  was  used  to  distinguish  among  the  isomers  in  the  model. 
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B.  ION-ELECTRON  RECOMBINATION 

Ions  disappear  by  either  cation-electron  or  cation-anion  recombination.  Cation-anion 
recombination  rate  coefficients  are  about  two  orders  of  magnitude  slower  than  cation-electron 
recombination  rate  coefficients.  There  are  no  good  measurements  of  anion  concentrations  in 
sooting  flames  but  there  is  evidence  for  their  presence  [32-34].  They  have  not  been  included  in  this 
mechanism;  their  inclusion  could  increase  the  cation  concentrations  but  calculations  indicated  that 
they  have  little  effect  on  the  calculated  ion  profiles.  In  choosing  product  channels  for  large  cation- 
electron  dissociative  recombination  reactions,  only  molecules  observed  by  Bockhorn  et  al.[35]  were 
considered  as  products. 


We  estimate  the  rate  of  ion  recombination,  a,  by  the  equation  for  the  rate  of  collision  of 
electrons  with  particles  [36]: 


a 


Tzd^ 

~T’ 


(  ( 


Tim 
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1  + 


(2ti  ef^d)kT 


(1) 


in  which  d  =  ion  diameter,  mg  =  electron  mass,  and  cq  =  dielectric  constant  of  free  space.  Ion 
diameters  were  calculated  from  experimental  ion  mobilities  by  use  of  the  Langevin  equation  for  ion 
mobilities.  Equation  (1)  gives  a  T“^^^  temperature  dependence  which  compares  favorably  with 
experiments  of  Ogram  et  al.[37]  for  The  values  calculated  by  Eq.  (1)  were  about  twice  the 

measured  values  of  Ogram,  so  the  calculated  values  were  divided  by  2.  From  Eq.  1,  reaction  rate 
coefficients  increase  with  the  size  of  the  ion. 


A  typical  dissociative  recombination  reaction  and  rate  coefficient  is: 

+  e  Cj2H8  -b  H  a  =  1.0  X  x  cm^/mol/s 


REACTION  MECHANISM 
The  ionic  path  starts  with  the  chemiionization  reaction: 

CH*  +  O  HCO+  +  e 

The  chemiion  EICO"''  (formyl  cation)  reacts  in  a  series  of  ion-molecule  reactions  to  produce  C3H3''’, 
an  ion  observed  in  large  concentrations  in  fuel  rich  flames.  To  minimize  the  number  of  reactions, 
only  H2,  C2H2,  C4H2  and  C3H4  were  considered  as  reactants  with  ions  in  the  ion  growth  scheme. 
As  ions  grew,  they  were  neutralized  through  ion-electron  recombination  reactions  at  rates  which 
increase  with  increasing  ion  mass,  producing  neutral  by-product  polycyclic  aromatic  hydrocarbons. 

In  the  kinetic  scheme,  experimental  profiles  for  03113'*'  (propargylium)  and  HgCg'*'  (phenyl 
cation)  were  used  because  good  experimental  profiles  for  the  HCO'*'  ion  were  unavailable  in  fuel 
rich  flames.  When  the  concentration  of  03113'*'  was  calculated,  using  the  neutral  mechanism  for  soot 
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forrnation  of  Frcnklach  and  Wang  [20],  chemiionization  reactions  to  form  HCO  and  reactions  to 
produce  C3H3‘^  from  HCO'^,  the  calculated  concentration  was  excessively  greater  than  the 
experimental  value. 

How  to  handle  C3H3''’  isomers  was  not  clear.  The  ion  has  been  identified  in  flames  by  mass 
spectrometry  which  provides  no  direct  information  on  structure.  Two  isomers,  the  linear 
propargylium  cation,  C3H3‘^,  and  the  cyclic  cyclopropenium  cation,  H3C3‘'',  have  been  considered 
as  potentially  important  [38].  Calculations  [39]  show  that  the  cyclic  cation  is  more  stable  (AjH)  than 
the  linear  cation  by  116  kJ/mol;  the  experimental  difference  is  105  kJ/mol  [40-41].  The  difference 
in  stability  decreases  with  increasing  temperature;  based  on  experimental  A^H's,  we  estimate  at  ^OOK 
that  AG(isomerization)  =  +99.4  kJ/mol,  giving  an  equilibrium  ratio  of  [linear]/[cyclic]  =  5  x  10“  , 
but  at  2000  K  the  ratio  is  0.03.  It  has  been  determined  experimentally  that  at  low  temperatures  the 
cyclic  isomer  is  not  reactive  with  acetylene  or  diacetylene,  but  the  linear  isomer  is  very  reactive  [42]. 
There  is  no  evidence  that  the  cyclopropenium  cation  isomerizes  to  the  propargylium  cation  in  the 
absence  of  encounter  complexes  with  neutral  molecules,  e.g.,  acetylene.[39]  Cameron  et  al.[39]  have 
suggested  that  in  sooting  flames  the  linear  isomer  may  be  the  most  abundant.  We  have  chosen  in 
the  mechanism  to  identify  the  linear  C3H3‘'’  ion  with  the  experimentally  observed  ion  in  the  absence 
of  evidence  to  the  contrary. 

In  developing  this  mechanism  only  ionic  species  which  have  been  observed  in  sooting  flames 
have  been  included  and  all  observed  ions  have  been  accounted  for.  This  is  a  more  stringent 
constraint  than  has  been  imposed  on  the  free  radical  mechanism  [20,43]. 


NUMERICAL  METHODOLOGY 

The  numerical  codes  developed  at  Sandia  National  Laboratories  by  Kee  and  coworkers  [44- 
46]  were  used  as  a  starting  point  for  modeling  the  ion  growth.  The  one-dimensional  flame  code  [44] 
was  chosen  because  it  has  a  very  accurate  and  efficient  integration  routine,  using  a  hybrid  time- 
integration/Newton-iteration  technique  to  solve  the  steady-state  equations  of  mass  and  energy.  For 
the  kinetics  calculations,  the  CHEMKIN-II  code  [45]  was  used,  because  it  has  capabilities  for  a 
variety  of  calculations  during  and  after  the  solution  of  the  problem,  and  it  can  be  integrated  very 
efficiently  into  the  flame  code.  For  transport  property  calculations,  the  code  of  Ref.  46  was  used, 
because  its  validity  has  been  tested  by  a  variety  of  investigators  and  can  be  integrated  efficiently  into 
the  flame  code.  Major  modifications,  however,  were  implemented  into  all  three  codes  to  allow  for 
non-Arrhenius  behavior  of  the  ion  kinetics,  ambipolar  diffusion  coefficients,  and  to  allow  the  user 
to  input  of  species  profiles  obtained  from  experiments. 

In  the  original  flame  code  [44]  there  was  provision  for  inputing  the  experimental  temperature 
profile.  This  approach  facilitates  the  calculations  greatly  because  the  need  to  model  heat  losses  in 
the  flat  flame  burner  experiment  is  eliminated.  The  code  was  further  modified  to  allow  for  the  input 
of  any  number  of  species  concentration  profiles  for  which  the  code  does  not  solve  but  simply  sets 
their  values  to  the  input  (experimental)  values.  This  modification  allows  for  the  "isolation"  of  any 
part  of  the  kinetic  scheme  so  that  uncertainties  related  to  the  chemistry  of  certain  species  do  not 
affect  conclusions  about  the  chemistry  of  the  species  of  interest.  In  the  present  study  the  purpose 
was  to  eliminate  uncertainties  related  to  the  neutrals  and  selected  ions.  For  the  electrons,  no 
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conservation  equation  was  solved;  their  concentration  was  determined  by  balancing  the  electric 
charge  between  the  electrons  and  ions.  This  is  consistent  with  the  assumption  of  ambipolar  diffusion. 
The  more  exact  method  of  solving  ion  diffusion  (Poisson's  equations)  has  been  used  by  Eraslan  and 
Brown  [19]. 

The  original  CHEMKIN-II  code  [45]  was  modified  to  accept  a  three-parameter  Arrhenius 
type  expression  for  the  elementary  reactions: 


k=AT^t\^{-EJRT)^B 


(2) 


where  B  is  a  constant.  This  rate  expression  was  used  to  fit  rate  data  for  ionic  reactions  which  were 
non-Arrhenius. 

The  molecular  transport  code  [46]  was  modified  to  use  ambipolar  diffusion  coefficients  for 
the  ions.  Specifically,  the  diffusion  coefficient  Dj  of  ion  i  in  a  gas  is  given  by  Einstein's  relationship: 

D.=2kT\iJe  (3) 


where  k  is  the  Boltzmann's  constant,  T  is  the  temperature,  /t;  is  the  ion  mobility  and  e  is 
the  elementary  charge.  The  ionic  mobility,  /tj,  at  P^  =  2.67  kPa  and  T^  =  273  K  was  calculated  [48] 
for  each  of  the  dominant  neutral  flame  gases  (j  =  neutral  gas  index)  by  first  determining  the 
coefficients,  aj  and  bj,  for  each  gas  in  the  Langevin  equation: 

by  fitting  this  to  experimental  values  of  PCAH  ion  diffusion  [49-51]  M'Wj  is  the  mass  of  the  ith  ion. 
TTie  coefficients  are: 


Gas 

a  . . 

b 

H2 

3472 

-0.499 

0.998 

O2 

743 

-0.487 

0.998 

CO2 

493 

-0.464 

0.998 

C2H2 

473 

-0.434 

0.999 

CO 

688 

-0.472 

0.998 

H20 

761 

-0.465 

0.998 

where  r^  is  the  coefficient  of  determination  of  the  fit  to  the  actual  data. 

The  ion  mobilities  calculated  above,  were  corrected  to  flame  pressures  and  temperatures,  by 
the  following  formula  (except  for  water),  where  P  is  the  flame  pressure  (kPa)  and  T  is  the  local 
temperature: 
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=  {PJP)  X  X 

For  water  vapor: 


The  ion  mobility  in  the  gas  mixture  was  calculated  by: 


(5) 

(6) 


(7) 


where  Xj  is  the  mole  fraction  of  gas  j  in  the  flame. 


RESULTS  AND  COMPARISON  WITH  EXPERIMENTS 

A  comparison  of  the  calculated  with  the  experimental  ion  concentration  profiles  [4]  for  the 
standard  acetylene/oxygen  flame  (equivalence  ratio  =  3.0,  P  =  2.67  kPa,  unburned  gas  velocity  = 
U  =  50  cm/s)  is  shown  in  Fig.  1,  and  a  comparison  between  the  maximum  calculated  and 
experimental  ion  concentrations  is  shown  in  Fig.  2.  The  experimental  profile  of  the  propargylium 
and  phenyl  cations  and  all  neutral  species  involved  in  the  ionic  mechanism  (ie,  acetylene,  diacetylene, 
allene  and  hydrogen)  were  input  profiles.  Because  the  ion  growth  paths  seemed  to  run  parallel  for 
even  and  odd  numbers  of  carbon  atoms,  and  generally  independent  of  one  another,  it  was  necessary 
to  use  both  an  odd  and  an  even  numbered  carbon  cation  as  input  ions.  Cation  reactions  with  C3H4 
couple  the  odd  and  even  ion  growth  pathways;  but  there  were  few  such  reactions  in  this  model. 

The  experimental  temperature  profile  was  input  to  the  model  to  avoid  the  necessity  of  heat 
transfer  calculations.  Temperature  was  an  unexpectedly  important  parameter  and  to  obtain 
reasonable  results,  the  experimental  temperature  profile  was  reduced  linearly  by  0%  at  the  burner 
surface  to  10%  at  the  temperature  maxima.  The  reduced  maximum  flame  temperature  was  1818  K. 
Thermodynamic  equilibrium  had  a  much  greater  effect  on  the  ion-molecule  reactions  than  expected 
(the  reverse  reaction  rates  were  calculated  from  forward  rates  and  equilibrium  constants),  this 
showed  up  as  a  very  strong  temperature  dependence  for  ion  concentrations.  At  higher  temperatures 
ion-molecule  reverse  reactions  became  more  important  because  entropy  effects  increase  with 
increasing  temperature.  This  same  effect,  of  course,  controls  neutral  mechanisms  as  well  and  may 
very  well  explain  the  high  temperature  portion  of  the  bell  shaped  curve  for  soot  formation  between 
about  1400  and  1900  K.  The  entropy  effect  is  a  more  probable  explanation  of  the  bell  shaped  curve 
than  the  explanation  based  on  the  increased  rate  of  soot  oxidation  with  temperature  [52].  A  number 
of  isomers  were  included  in  the  calculations,  but  only  the  most  abundant  calculated  cation  isomer 
for  any  ion  mass  is  plotted  in  Figs  1  and  2;  the  mass  spectrometer  is  blind  to  isomeric  structure.  It 
is,  of  course,  conceivable  that  the  most  thermodynamically  stable  isomer  is  not  formed  because  of 
reaction  kinetic  restrictions.  This  may  increase  or  decrease  the  agreement  with  experiment. 

For  most  of  the  ions  the  calculated  and  experimental  profiles  were  in  agreement  within  the 
accuracy  of  the  data.  The  ion  most  out  of  agreement  in  the  initial  calculation  was  C^ngH 
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(isomers  are  distinguished  by  the  order  of  C  and  H)  which  led  to  large  discrepancies  for  C-^2^9^ 
Cj3HgH3'^.  The  measured  value  for  exceeded  the  calculated  value  by  three  orders  of 

magnitude.  Reduction  in  the  heat  of  formation  of  by  10%  gave  the  agreement  reported 

in  Figs.  1  and  2. 

Of  concern  is  the  fact  that  some  of  the  calculated  concentration  profiles  show  an  increase 
as  the  burner  is  approached.  Double  cation  peaks  have  been  observed  experimentally  and  it  is  not 
clear  whether  these  are  related  or  not.  The  calculated  ion  concentration  increase  toward  the  burner 
is  due  to  too  a  large  temperature  effect  and  the  use  of  a  flame  temperature  profile  that  is  not 
accurate  near  the  burner  surface.  This  situation  could  be  corrected  by  modifying  the  input 
temperature  profile  even  more  than  we  have  already  done. 

To  better  evaluate  the  agreement  between  experimental  and  calculated  concentrations,  the 
maximum  concentration  ratios,  log^Q  (experimental  concentration/calculated  concentration)  are 
presented  in  Fig.  2.  The  isomer  which  had  the  highest  calculated  concentration  is  compared  with 
the  experimental  data.  In  general,  calculated  concentrations  are  within  a  decade  of  the  experimental 
data.  The  source  of  the  original  disagreement  with  the  C^j^HgH'^  cation  has  been  considered  by 
carefully  reviewing  the  thermodynamics,  including  the  possibility  of  other  more  stable  isomers, 
without  arriving  at  a  better  solution  than  reducing  the  aj-H  by  10%.  It  appears  that  there  is  little 
latitude  for  adjusting  the  rate  coefficients.  The  experimental  results  could  be  called  into  question, 
but  the  major  error  is  probably  in  the  thermodynamics.  This  will  be  pursued  further  as  we  extend 
the  model  to  larger  ions;  development  of  the  thermodynamics  and  reaction  kinetics  for  larger  ions 
has  already  been  completed. 

The  major  ion  growth  paths  (Fig.  3)  were  identified  by  integrating  ion  fluxes  over  the 
distance  grid.  This  flow  diagram  shows  the  weak  coupling  between  the  odd  and  even  ion  growth 
paths  in  this  model. 

Oxidation  of  ions  has  yet  to  be  added  to  the  mechanism.  Unlike  the  neutral  mechanism, 
wherein  oxidation  removes  the  neutral  growth  species  from  the  growth  chain,  for  ions,  oxidation 
simply  removes  carbon  atoms  from  the  ion,  making  a  smaller  ion,  or  removes  hydrogen  from  the 
growing  ion  at  a  greater  rate.  Removal  of  hydrogen  may  in  fact  contribute  to  ion  growth  toward 
soot  which  has  a  C/H  ratio  of  about  10,  while  for  the  largest  ion  treated  here  the  ratio  is  only  1.4. 
Since  the  ions  considered  in  the  model  have  a  decreasing  free  energy  per  carbon  atom  with 
increasing  size  [53],  it  is  assumed  that  the  charge  will  remain  on  the  larger  hydrocarbon  products, 
thus  e.g.: 

C13H9+  +  OH  -  C13H8+  +  H2O 
would  be  expected  to  occur,  rather  than: 

C13H9+  +  OH  -  q3H2  -1-  H3O+ 

This  of  course  requires  further  consideration  with  other  oxidizers  and  other  products. 
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SUMMARY  AND  CONCLUSIONS 

The  current  model  of  ion-growth  correctly  predicts  both  the  shape  of  the  ion  profiles  and  the 
magnitude  of  the  ion  concentrations.  Although,  adjustments  in  the  temperature  profile  (10%)  and 
the  enthalpy  of  formation  of  the  ion  (10%)  were  used  to  obtain  closer  agreement  with 

experiment,  the  results  were  very  encouraging  for  such  a  simple  model.  The  changes  in  the 
mechanism  are  physically  reasonable,  considering  the  uncertainties  in  the  thermodynamics  and 
temperature  measurements. 
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Fig.  1. 
Fig.  2. 

Fig.  3. 


FIGURE  CAPTIONS 


Calculated  (solid  lines)  and  Experimental  (dashed  lines)  Ion  Concentration  Profiles 
for  Ions  up  to  Mass  165  u. 

Horizontal  Bar  Graph  of  LogjQ  (Maximum  Experimental  Concentration  /  Maximum 
Calculated  Concentration)  for  Cation  Isomers  having  the  Highest  Calculated 
Concentration. 

Ion  Flow  Diagram  Based  on  Integrated  Flux  (over  Calculation  Distance  Grid)  for 
Even  and  Odd  Cations. 
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REVIEWER  A,  No.  362 


There  is  much  of  interest  in  this  paper  but,  ovcrail,  not  enough  new  or  definitive  material  for  me 
to  be  sure  that  the  paper  should  be  accepted  for  the  Symposium.  The  authors  have  something  of 
a  dilemma  in  trying  to  describe  their-  approach  to  the  reaciion-ldnetic  modeUing,  as  well  as 
presenting  a  meaningful  analysis  of  the  results,  within  the  constraints  of  a  Symposium 
submission  -  unfortunately,  on  neither  front  docs  the  paper  succeed  in  reaching  the  standard 
required. 

Successes  in  the  modelling  of  neutral  chemistry  in  combustion  systems  have  relied  heavily  on 
the  use  of  sensitivity  analysis,  in  order  that  the  overall  model  can  be  see  in  temts  of  its  mo.st 
crucial  parameters.  It  is  one  of  the  deficiencies  of  this  paper  that  there  is  no  coherent  approach 
to  parameter  evaluation,  with  parameters  being  adjusted  on  a  more  or  less  ad  hue  basis  to 
achieve  desired  results.  A  case  in  point  is  the  adjustment  of  the  for  CnHjH*.  The  enormous 
sensitivity  of  ion  concentrations  to  heats  of  formation  was  shown  previously  by  Eraslan  and 
Brown  (Ref.  19)  and  one  can  expect  a  more  quantitative  assessment  here. 

Thermochemistry  also  appears  in  some  other  issues.  Forward  rate  constants  have  apparendy  been 
adjusted  so  that  the  reverse  rate  calculated  from  the  equilibrium  constant  is  not  "too  fast"  to  be 
consistent  with  calculated  rates,  but  the  same  theory  has  been  used  to  estimate  the  forward  rates 
in  the  first  place.  The  alternative  of  adjusdng  the  theimocheraical  data  is  not  discussed.  A  clear 
hierarchy  needs  to  be  established  here  because  it  is  not  clear  that  the  estimation  of  the 
thermochemical  data  is  in  effect  any  more  reliable  than  the  c.stimation  of  rate  constant  in  this 
system. 

A  substantial  part  of  the  text  is  devoted  to  describing  how  ambipolar  diffusion  coefficients  have 
been  estimated.  Eraslan  and  Brown  (Ref.  19)  argued  that  it  was  not  nece.ssay  to  allow  for 
ambipolar  diffusion  and  the  authors’  earlier  work,  (Ref  4.)  shows  no  evidence  to  the  contrary.  At 
the  very  least,  the  present  Discussion  should  include  some  comments  and  results  from  the 
present  modelling  on  this  issue. 

The  input  temperature  profile  is  another  (distributed)  parameter  which  needs  to  be  assessed  more 
rigorou-sly.  The  authors  have  adjusted  the  profile  in  one  pan  of  the  flame  in  order  to  obtain  a 
better  fit  with  the  chosen  data  sec  Without  a  more  detailed  evaluation  of  the  sensitivity 
coefficients  for  important  factors  such  a.s  thermochemical  and  kinetic  parameters,  it  isj^SSsible  to 
evaluate  what  is  actually  being  achieved  by  this  manipulation  other  than  curve-fitung.  This  is 
particularly  so  when  one  considers  that  the  maximum  temperature  reported  for  the  Aerochem 
burner  was  already  lOOK  lower  than  others  had  obtained  (sec  Ref.  4)  and  further  downward 
adjustment  seems  contrary  to  the  weight  of  experimental  evidence. 

There  is  no  way  that  anyone  reading  this  paper  could  ever  duplicate  the  modelling  which  has 
been  reported  as  mechani.sm  and  thermochemistry  arc  not  reported  in  full.  Nor  is  there  any 
comparison  with  the  modelling  undertaken  by  Eraslan  and  Brown  (Ref.  19)  in  terms  of  the 
mathematical  procedures,  the  reaction  mechanism,  or  the  results.  While  I  look  forward  to  future 
work  by  the  authors,  the  present  paper  really  only  whets  the  apoeute. 


REVIEWER  B,  No.  428 


This  work  deals  with  the  question  whether  one  can  characterize  soot  formation  as  a  non- 
equilibrium  process,  with  respect  to  the  equilibrium  concentrations  of  acetylene  (a  potential 
reactant)  and  of  soot  (the  product),  as  compared  to  the  experimental  CjHj  and  soot 
concentrations.  Although  one  knows  that  soot  formation  is  kinetically  controlled  and  not 
thermodynamically,  it  is  useful  to  learn  from  the  figures  how  large  the  discrepancies  are. 

However,  the  paper  suffers  from  the  fact  that  the  results  of  complicated  quantitative  -j 
equilibrium  calculations,  which  cannot  be  reproduced  by  the  reader,  are  mixed  with  very  2 

qualitative  concepts  such  as  'driving  force'.  Deviations  from  equilibrium  for  characteristic 
reactions  could  have  been  expressed  in  terms  of  thermodynamic  affinities,  so  that  the  reader  gets 
a  more  quantitative  idea,  how  far  from  equilibrium  some  important  reactions  are.  From  the 
figures,  one  only  gets  the  impression,  that  the  real  'flame  world'  and  the  equilibrium  state  occupy 
areas  in  the  diagrams,  which  are  at  quite  a  distance  from  each  other.  The  authors  speak  of 
'driving  force'  through  supra-equilibrium  concentrations  of  H  and  ions.  In  chemical  3 
thermodynamics  'driving  force'  can  be  quantified  by  affinity  (de  Donder)  which  is  defined  for  one 
stoichiometric  equation  as  -TViji;.  In  an  ideal  system  p^=  pf+RT  Incj.  Through  this  the  'driving 
force'  is  coupled  to  concentrations.  Now,  the  concentration  of  H  is  many  orders  of  magnitude  4 
greater  than  that  of  ions.  Would  the  authors  comment  on  the  absolute  values  of  these  driving 
forces  ? 

Extended  equilibrium  calculations  have  been  performed.  However,  the  equilibrated 
reactions  considered  are  not  given,  only  a  number  of  species.  It  does  not  become  clear,  what 
reaction  is  considered,  for  example,  for  circumcoronene.  Is  it  9  +  54  C(Graphite)  or  +  54 

C(Soot)  =  C54Hig?  The  same  applies  for  other  species. 

There  are  other  shortcomings  and  non-adequate  formulations  which  the  authors  might 
consider,  in  order  not  to  convey  the  impression  that  they  are  going  to  oversimplify  things.  For 
example: 

They  say  that  the  only  flame  property  that  changes  at  the  soot  threshold  is  the  ionic  0 
character.  This  statement  is  far  too  simple  and  because  of  its  simplicity  it  is  not  correct.  There  are 
a  great  number  of  properties  concerning  flame  structure  and  composition  that  change,  although 
gradually. 

Introduction,  line  5  from  below:  The  presence  of  anions  is  not  only  indicated  by  a  change 
in  the  measured  ion  recombination  coefficient.  From  mass  spectrometry  of  flame  ions  one  knows  7 
exactly  what  kinds  of  negative  ions  are  present  for  different  fuels  when  soot  formation  starts  in 
a  premixed  flame. 

It  should  not  only  be  justified  by  quotations  (36,  37)  but  discussed  by  the  authors  why 
they  hold  a  super-equilibrium  concentration  of  H  to  be  a  'driving  force'  for  soot  formation.  There  g 
are  a  number  of  observations  of  a  presumed  influence  of  H  on  the  formation  of  benzene,  PAH 
and  soot.  But  the  results  are  not  conclusive  and  there  seems  to  be  a  strong  influence  of 
temperature. 

P.  1 1,  line  11.  What  is  the  meaning  of  the  statement  'Only  one  ion  is  required  per  soot 
particle'.  Required  for  what  ?  If  it  shall  be  expressed  that  one  positive  charge  is  required  on  the  9 
soot  particle  for  its  existence  as  a  charged  species,  it  is  trivial.  But  I  presume  that  the  authors  do 
not  have  this  in  mind. 
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Some  further  comments  and  recommendations: 

Introduction,  p.4,  line  6:  What  is  meant  by  the  statement,  that  charge  is  'transferred  from  a  larger  1  0 
moiety  to  a  small  moiety  by  charge  transfer  or  ion-molecule  reactions'  ?  What  is  the  exact 
meaning  of 'moiety'  in  this  context  ? 

Homann  is  not  correctly  quoted  by  reference  [10],  since  that  paper  deals  with  the  growth  of  -j 
charged  soot  particles  and  not  with  the  growth  of  small  molecular  ions  to  larger  ions. 

The  name  perinaphthenyl  is  not  a  lUPAC  name  for  any  C13H9",  There  are  several  possible 
structures  for  this  molecular  formular.  What  structure  is  indicated  by  the  formula  Ci3H^H3*  ?  ^ 

P.  7,  line  14:  In  which  way  were  isomers  distinguished  by  the  order  of  "C"  and  "H"  in  the  ^3 
formulas  ? 

The  experimental  profiles  of  C3H3^  and  of  the  phenyl  ion,  which  have  been  put  into  the  14 
calculations,  should  also  be  given  in  fig.  1 . 

REVIEWER  C,  IMo.  331 

3.  The  abstract  should  mention  that  ion  profiles  were 
computed  with  temperature  and  other  species  profiles 
imposed  from  the  experiment.  It  should  qualify  that  the  1 
temperature  profile  was  also  reduced  to  improve  the 
comparison,  since  by  the  authors’  admission  this  was 
unexpected . 

5.  The  model  is  new  and  its  verification  is  limited;  it  may 
be  prudent  to  delay  this  paper *s  publication.  The 
modifications  to  Sandia’s  flame  code  seem  reasonable,  yet 
critical  reviews  of  the  extensions  of  thermodynamic 
properties  and  Langevin  theory  have  not  been  completed; 
they  are  referenced  as  "submitted”  or  "in  preparation." 

The  accuracy  of  the  predictions  of  ion  profiles  in  this 
one  flame  may  be  promising,  but  many  questions  remain. 

The  proposed  mechanism  starts  with  chemiionizat ion  to 
form  HCO+,  eventually  producing  C3H3+  by  a  series  of 
reactions;  further  cation  growth  proceeds  by  the 
mechanism  in  fig.  3.  However,  the  present  results  skip 
this  initiation,  using  the  experimental  profile  of  C3H3+ 
directly,  because  HC0+  data  were  not  available,  and 
including  the  proposed  initiation  steps  produced 
"excessive"  concentrations.  This  leaves  the  initiation 
of  the  proposed  mechanism  in  doubt.  It  seems  that 
several  adjustments  were  involved  in  producing  the 
agreement  in  fig.  1.  It  is  not  obvious  why  a  factor  of  2 
is  necessary  in  the  electron  collision  rate.  In 
addition,  the  experimental  temperature  profile  and  the 
enthalpy  of  formation  of  one  ion  were  adjusted  by  10%  to 
improve  the  agreement.  The  need  to  reduce  the  flame 
temperature  was  "unexpected"  and  is  a  concern,  yet  there 
is  no  real  discussion  of  why  this  is  reasonable 
(experimental  accuracy?)  or  how  sensitive  the  predictions 
are  to  the  temperature  profile. 
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Given  these  concerns,  I  hesitate  to  judge  the  predictions 
accurate,  or  the  proposed  mechanism  ready  for  acceptance. 

9.  The  authors  should  be  able  to  address  the  questions 

regarding  some  of  their  parameter  adjustments.  However, 
justification  of  skipping  the  initiation  steps  will  take 
more  effort.  In  any  case,  this  work  is  relevant,  and 
further  effort  will  surely  make  this  a  fine  paper. 

RESPONSE  TO  REVIEWERS'  COMMENTS 


The  reviewer  is  identified  by  letter  on  his  review  copy  and  the  reviewer's  comment  is 
identified  by  a  number  on  the  right  of  the  review.  References  numbers  refer  to 
references  in  subject  paper. 


Reviewer  A.  No.  362 


1.  This  material  has  not  been  previously  published. 

2.  Sensitivity  analysis  takes  many  forms,  from  formal  mathematical  treatments  to 
inspection  of  how  specific  changes  affect  the  result.  We  used  the  latter.  That 
is  the  means  by  which  we  determined  the  desirability  of  adjusting  the 
temperature  profile  and  A,  H  for  one  ion.  Paper  length  limitations  do  not  allow 
presentation  of  the  many  versions  of  the  model  tested  in  arriving  at  the 
mechanism. 

3.  Equilibrium  is  determined  by  the  relative  forward  and  reverse  rates.  The 
exothermic  ion-molecule  reactions  are  assumed  to  have  no  activation  energies, 
consistent  with  experimental  data  on  small  ions.  The  adjustment  is  a  means  of 
determining  the  rate  for  the  forward  reaction  when  the  reverse  reaction  has  a 
negative  AG  at  flame  temperatures.  It  is  standard  procedure  to  calculate  a  rate 
coefficient  from  the  reverse  rate  coefficient  and  the  equilibrium  constant.  It  is 
treated  a  little  differently  for  ion-molecule  reactions  where  the  temperature 
effect  is  non-Arrhenius.  If  this  were  not  done  some  rates  would  exceed 
collision  rates.  The  question  of  hierarchy  does  not  arise,  this  correction  must 
be  done  -  it  is  not  arbitrary  as  assumed  by  the  reviewer. 

4.  Eraslan  and  Brown  do  not  argue  against  the  need  for  ambipolar  diffusion,  they 
treat  it  in  a  more  exact  way,  by  the  use  of  Poisson's  equations.  Ambipolar 
diffusion,  or  a  more  exact  treatment  as  used  by  Eraslan  and  Brown,  must  be 
included  for  a  neutral  plasma  made  up  of  cations  and  electrons.  It  accounts  for 
the  difference  in  diffusion  coefficients  of  electrons  and  cations  and  the 
interaction  of  the  electric  field  produced  by  the  diffusion  of  the  two  charges 
(positive  and  negative).  It  makes  a  small  difference  whether  one  uses  the 
treatment  of  Eraslan  and  Brown  or  ambipolar  diffusion,  and  neither  makes  a 
significant  difference  in  this  particular  modeling  result.  The  authors  earlier  work 
shows  the  importance  of  ambipolar  diffusion  in  flames  [Eighth  Symposium  on 
Combustion,  1962,  pg  184]. 
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5.  Adjustment  of  the  temperature  profile  was  made  only  after  consideration  of  the 
sensitivity  to  thermodynamic  and  reaction  kinetics.  The  need  for  this 
adjustment  is  of  some  concern  to  the  authors.  See  the  Introduction  for  the 
philosophy  of  our  approach.  The  maximum  temperature  for  the  data  under 
discussion  was  the  temperature  AeroChem  measured  on  the  MIT  burner  and  is 
not  already  1 00  K  lower  than  reported  by  others  [see  Ref.  4  of  the  paper].  The 
reason  for  the  differences  in  temperatures  reported  in  Ref.  4  are  well  known  - 
see  Ref.  4. 

5a.  Modeling  often  reduces  to  curve  fitting,  it  is  a  very  powerful  analytical  tool. 

6.  Symposium  restrictions  do  not  permit  reporting  thermochemistry  in  full  (this  is 
an  extensive  paper  which  stands  alone),  nor  is  it  normal  practice  to  report 
complete  thermochemistry  data  with  kinetic  mechanisms.  The  equations,  with 
a  condensed  mechanism,  is  equivalent  to  reporting  the  complete  mechanism. 

7.  Symposium  restrictions  do  not  permit  the  paper  length  which  would  be  required 
to  compare  the  modeling  and  mathematical  procedures  used  here  with  those  of 
Eraslan  and  Brown  for  non-sooting  flames,  nor  would  it  be  appropriate,  or 
useful. 


Reviewer  B.  No.  428 

1 .  The  Sandia  flame  code  is  a  well  established  code  for  this  kind  of  work.  The 
question  of  "how  the  rate  constants  describe  the  neutral  structure  is  not 
germane  because  we  start  with  two  known  ion  profiles  to  avoid  the  complexity 
of  the  neutral  chemistry.  As  pointed  out,  page  4,  the  neutral  chemistry  as  used 
by  Frenklach  over-estimates  the  neutral  radical  concentrations  and  thus  the 
initial  chemiion,  HCO'*',  concentration.  That  is  one  of  the  reasons  we  start  with 
experimental  cation  profiles.  When  the  ionic  mechanism  is  demonstrated  to  be 
adequate  for  a  number  of  different  flames,  we  will  consider  a  neutral 
mechanism  to  produce  chemiions.  Brown  and  associates  have  done  it  for  non¬ 
sooting  flames. 

2.  There  is  an  extensive  set  of  experimental  data  on  this  flame  obtained  in  a 
number  of  different  laboratories.  We  know  of  no  other  such  well  documented 
flame. 

3.  H2,  C2H2,  C4H2,  and  C3H4,  see  text,  page  4  bottom. 

4.  Some  examples  are  given  on  page  3,  space  limitation  prevents  presentation  of 
the  complete  reaction  scheme,  but  the  equations  should  permit  one  to  calculate 
them  and  compare  them  with  the  examples  given.  The  thermochemistry  is  of 
course  required,  our  program  uses  the  complete  JANNAF  tables  plus  a  very 
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large  number  of  ions  for  which  we  developed  the  data.  We  will  be  delighted  to 
furnish  a  copy  to  anyone  interested  prior  to  publication. 

5.  All  observed  ions  in  the  flame  under  study  are  in  the  mechanism  and  protonated 

polyynes  have  been  handled  the  same  as  any  other  ion,  we  have  not  checked 
to  see  if  they  are  in  equilibrium  or  not.  We  are  not  as  sure  as  the  reviewer  that 
their  equilibrium  is  "well  established".  Michaud,  et.al.  (1981)  assumed 
equilibrium  to  calculate  equilibrium  constants,  and  I  am  not  aware  that  Homann, 
et.al.  (e.g.,  1990),  who  observe  more  polyynes  than  we  do,  ever  demonstrated 
they  are  in  equilibrium  -  they  also  assumed  equilibrium.  The  calculated 
concentration  profile  for  is  about  10®  times  in  excess  of  the  equilibrium 

value.  This  is  typical  of  all  flame  ions,  and  is  explained  by  the  now  well 
established  phenomenon  of  chemiionization. 

6.  Proton  transfer  reactions  have  not  been  included  in  this  kinetic  scheme  to  date. 
Their  contribution  appears  to  be  small  because  of  the  low  concentration  of  large 
neutral  species. 

7.  The  source  of  all  ions  is  the  initial  chemiionization  step  which  produces  an 
electron  for  each  cation.  Dissociative  recombination  eliminates  one  of  each 
unless  the  electron  is  attached  to  an  entity  to  produce  an  anion  and  then 
recombination  removes  both  a  cation  and  an  anion  from  the  system.  As  a 
simplification,  the  program  assumes  an  equal  number  of  electrons  as  cations. 
This  is  a  programmers  means  of  producing  an  electron  every  time  a  cation  is 
produced  and  removing  an  electron  every  time  a  cation  is  removed.  This 
procedure  maintains  internal  consistency. 

8.  in  our  flux  analyses,  we  have  not  concentrated  on  answering  this  question. 
Obviously,  if  we  travel  with  a  gas  volume  leaving  the  burner  surface  and 
moving  through  the  flame  front,  ion  concentration  will  rise  and  fall  as  the  net 
rate  of  production  of  an  ion  increases  and  then  declines. 

9.  This  is  not  an  extended  abstract  of  a  paper  published  or  to  be  published 
elsewhere,  but  is  obviously  a  condensation  of  a  large  amount  of  work. 
Separate  papers  will  be  published  on  the  thermochemistry,  kinetics,  and 
extension  of  Langevin  theory,  as  is  appropriate. 

10.  Moiety  is  misused  in  this  text  and  can  be  replaced  by  "species".  This 
statement  describes  how  one  would  have  to  interpret  Homann's  thesis  that  the 
observed  (small)  ions  are  produced  from  the  ionization  of  particles  or  large 
molecules.  The  meaning  is  that,  for  a  large  ion  to  produce  a  small  ion,  either 
simple  charge  transfer  (an  elementary  ion-molecule  reaction)  or  transfer  of  a 
charged  fragment  (moiety)  by  an  ion-molecule  reaction  must  occur  between  the 
large  ion  or  charged  particle  and  a  small  neutral  species. 
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1 1 .  The  processes  are  related  in  terms  of  the  source  of  the  charge  on  particles. 
Homann,  in  the  past,  has  rejected  conventional  chemiionization  processes  as 
the  source  of  this  charge,  although  he  has  proposed  different  chemiionization 
processes  for  which  there  is  absolutely  no  evidence  (see  e.g..  Ref.  10,  page  72, 
last  left  column  paragraph).  The  first  reference.  Ref.  9,  to  Homann  is  for  this 
paragraph  and  the  second  reference  is  to  demonstrate  Homann's  recognition  of 
the  importance  of  chemiions,  -  -  it  is  approached  from  the  ionic  rather  than 
the  neutral  side."  The  reviewer  might  appreciate  this  point  more  by  reading 
Ref.  6  of  the  submission,  especially  page  497,  left  column.  Reference  could  be 
left  out,  it  may  be  too  subtle, 

12.  Perinaphthenyl  is  the  name  used  by  Stein,  U.S.  NIST,  and  other  U.S. 
authorities,  e.g.,  Kebarle.  It  is  a  reference  to  the  three  fused  six-  membered 
aromatic  rings.  A  number  of  isomers  of  this  ion  have  been  considered, 

13.  Arbitrary,  what  else  would  be  possible?  It  turns  out  in  practice  that  the 
simplest  relationships  e.g.,  C13H9,  H9C13  to  more  complex  relationships  C13H7H2 
are  determined  by  the  order  in  which  we  identified  the  isomer.  We  will  help  the 
reader  by  naming  them  on  the  figure.  We  use  this  relationship  to  communicate 
with  the  computer  program.  For  our  data  base  we  use  a  number  in  parenthesis 
following  the  usual  order  of  symbols;  the  computer  program  does  not  recognize 
this. 

14.  Reasonable  suggestion.  We  were  trying  to  simplify  what  is  an  already  complex 
figure. 

Reviewer  C.  No.  331 

1 .  We  agree. 

2.  It  is  the  nature  of  all  new  works  to  require  further  testing. 

3.  Correct,  but  we  believe  that  this  work  stands  sufficiently  alone  to  be  published, 
and  is  of  sufficient  importance. 

4.  A  very  large  number  of  questions  remain,  which  should  not  be  surprising  since 
this  represents  a  major  step  toward  the  solution  of  a  problem  for  which  the 
thermochemistry  and  reaction  kinetics  were  all  missing  when  we  started.  Many 
questions  still  remain  in  the  neutral  proposed  mechanisms  of  soot  formation, 
after  considerably  more  effort  by  a  far  larger  number  of  people.  As  an  example, 
it  is  still  debated  in  the  literature  how  benzene  is  formed  in  the  neutral/free 
radical  models. 

5.  The  proposed  mechanism  that  is  in  doubt  is  that  part  of  the  neutral  mechanism 
required  to  produce  the  reactants  for  chemiionization,  these  were  taken  from 
Frenklach's  neutral  mechanism  of  soot  formation.  The  indication  is  that  the 
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reactions  in  the  neutral  mechanism  are  too  fast.  Our  objective  is  to  examine  the 
ionic  growth  mechanism  unencumbered  by  idiosyncrasies  in  the  neutral 
mechanisms,  so  we  modified  the  Sandia  program  to  avoid  involvement  in  the 
same  mechanism  everyone  else  is  working  on.  Hopefully,  by  the  time  we  need 
the  neutral  mechanism  to  complete  our  model,  the  folks  working  on  it  will  have 
reached  some  kind  of  consensus. 

6.  It  is  valid  procedure  to  make  reasonable  adjustments,  see  e.g.,  Gardiner,  J. 
Phys.  Chem.  83.  37  (1979). 

7.  As  stated,  to  reconcile  theory  with  experiment.  This  was  not  an  adjustment  to 
fit  the  model  to  experiment,  as  the  reviewer  seems  to  think,  but  a  means  of 
developing  a  reliable  data  base.  Incidentally  our  sensitivity  analysis  showed 
very  little  sensitivity  to  the  recombination  rate  of  ions  with  electrons. 

8.  Modeling  is  a  process  by  which  one  adjusts  parameters  within  reason  to  fit 
experiments  in  order  to  better  understand  the  phenomena.  Fitting  only  two 
parameters  to  obtain  such  good  agreement  is  very  encouraging.  It  should  be 
rated  as  "excellent"  when  compared  to  the  agreement  reported  for  neutral 
models  of  soot  nucleation  where  it  is  common  practice  to  include  species  in  the 
mechanism  that  are  not  observed  in  the  flame  and  where  some  of  the  major 
products  which  are  observed  are  not  predicted.  We  are  concerned  about  the 
temperature  correction  and  are  trying  to  understand  why  it  was  necessary,  or 
whether  it  is  really  necessary. 

9.  Neither  would  we  propose  the  mechanism  as  "ready  for  acceptance".  It  would 
be  inappropriate  to  propose  any  current  mechanism  of  soot  formation  as  "ready 
for  acceptance"  while  there  is  so  much  disagreement  on  some  of  the  most 
elementary  aspects  of  the  problem.  Presumably  we  will  continue  to  improve 
the  model  as  we  extend  it  to  large  ions  and  to  other  flames.  How  much  of  the 
published  work  on  combustion  does  the  reviewer  rate  as  "ready  for 
acceptance"? 

10.  The  initial  steps  are  a  problem  all  to  themselves  as  we  described  above  and  will 

be  addressed  in  a  separate  paper;  they  have  no  relevance  to  the  growth 
process,  so  long  as  there  is  a  starting  point  for  the  growth  process  -  which  we 
have  taken  as  and  HgCg'^  profiles. 
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ABSTRACT 

Observed  acetylene  and  soot  mole  fractions  in  laboratory  premixed  flames  are  compared 
with  equilibrium  calculations  for  a  variety  of  fuels  under  a  variety  of  conditions.  In  all  of  the 
flames,  except  the  standard  acetylene/oxygen  flame  (0  =  3.0,  P  =  2.67  kPa,  u  =  50  cm/s),  the 
experimental  concentrations  at  C/O  <  1.0  are  in  great  excess  of  equilibrium,  sometimes 
exceeding  a  factor  of  10^^.  Equilibrium  acetylene  and  soot  are  only  formed  at  C/O  >  1.0.  In 
contrast,  in  the  standard  acetylene/oxygen  flame,  while  the  experimental  concentration  of 
acetylene  remains  greater  than  the  equilibrium  value,  the  experimental  soot  concentration  is 
more  than  an  order  of  magnitude  less  than  equilibrium.  The  observed  supra-equilibrium 
acetylene  and  soot  production  may  be  due  to  either  carbon  dioxide  production  being  greater  than 
parallel  hydrocarbon  growth  reactions,  leading  to  an  effective  C/O  ratio  greater  than  determined 
by  the  mixture  ratio;  or  the  mechanism  involves  the  production  of  supra-equilibrium 
concentrations  of  chemiions  (or  excess  hydrogen  atoms).  The  same  mechanism  would  not  hold 
for  the  acetylene/oxygen  flame  where  soot  concentrations  are  less  than  equilibrium.  The  sub¬ 
equilibrium  soot  concentration  may  be  due  to  the  lower  ion  concentration  in  these  very  fuel  rich 
flames.  The  total  ion  concentration  becomes  smaller  as  the  equivalence  ratio  increases. 


INTRODUCTION 

Our  motivation  for  studying  the  thermochemistry  of  soot  formation  is  to  understand  the 
ionic  mechanism  of  soot  formation.  At  soot  threshold,  the  only  flame  property  to  change  is  the 
ionic  character  of  the  flame.  In  a  fuel  rich  acetylene/oxygen  or  ethylene/oxygen  flame  [1-3]  the 
total  cation  concentration  decreases  with  increasing  equivalence  ratio  until  at  some  critical 
equivalence  ratio,  near  soot  threshold,  the  total  ion  concentration  increases  with  equivalence 
ratio.  The  larger  cations  increase  in  concentration  and  size  at  the  expense  of  the  smaller  cations, 
and  new  large  cations  appear.  Anions  are  observed  simultaneously  with  these  changes  in  cations, 
as  indicated  by  a  change  in  the  measured  ion  recombination  coefficients  [1].  Are  ions  a  cause  or 
effect  of  soot  formation?  Since  they  are  the  only  property  that  changes  at  the  soot  threshold,  it 
might  be  argued  that  they  are  the  cause.  Why  then  does  soot  appear  when  the  ion  concentration 
is  falling?  Are  the  large  ions  the  precursors  of  soot  or  is  soot  the  precursor  of  ions,  or  is  the 
simultaneous  appearance  of  soot  and  large  ions  an  incidental  occurrence?  The  increase  in  ion 
concentration  beyond  the  soot  threshold  might  be  explained  as  attachment  of  electrons  to  large 
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molecular  species  or  incipient  soot  particles  that  attach  electrons  to  produce  anions.  Ion-ion 
recombination  is  about  two  orders  of  magnitude  slower  than  ion-electron  recombination  and  this 
could  account  for  an  increase  in  observed  cation  concentrations. 

The  increase  in  ion  concentration  beyond  the  soot  threshold  has  been  interpreted  as 
evidence  against  the  ionic  mechanism  of  soot  formation  [3,4],  However,  the  formation  of 
observed  ions  from  soot  particles  or  larger  ions  has  been  eliminated  as  a  possible  explanation  of 
the  source  of  flames  ions  [2,5].  Ions  become  increasingly  stable  with  increase  in  size,  so  there  is 
no  driving  force  to  produce  the  observed  ions  from  charged  soot  particles  or  larger  ions.  As  an 
initial  step  to  e.xamining  these  questions  in  more  detail,  we  examine  here  the  thermochemistry  of 
soot  formation. 

When  the  ratio  of  C/0  is  greater  than  1  there  is  insufficient  oxygen  to  form  CO  and 
carbon  must  be  produced.  Thus  soot  formation  is  more  a  matter  of  stoichiometry  than 
thermochemistry.  Gaydon  and  Wolfhard  discuss  the  thermochemistry  of  soot  formation  in  detail 
in  their  book  [6,7],  in  terms  of  the  Boudouard  reaction: 

CO2  +  C  (solid)  =  2CO  (1) 

for  which  the  equilibrium  constant  is: 

Kp  =  P2(C0)/P(C02)  (2) 

This  equilibrium  becomes  effective  when  solid  carbon  is  produced  and  this  changes  the  actual 
C/0  ratio  by  a  slight  amount.  Porter,  1954,  in  a  review  paper  on  the  mechanism  of  carbon 
formation  [8]  continues  the  discussion  of  the  effect  of  equilibrium/stoichiometiy  on  soot 
formation. 

Gay  et.al.[9],  1961,  performed  complete  equilibrium  calculation  of  the  appearance  of  soot 
and  compared  their  calculations  with  experiment.  For  ethylene/oxygen  flames  they  predicted  soot 
formation  at  C/O  =  1.0  with  an  adiabatic  flame  temperature  of  2345  K.  They  observed  soot 
formation  at  this  C/O  ratio.  For  acetylene  the  C/O  ratio  was  calculated  to  be  1.08  with  an 
adiabatic  flame  temperature  of  3275  K.  They  observed  soot  at  C/O  =  1.16.  Their  observations 
appeared  to  confirm  the  equilibrium  nature  of  soot  formation.  As  will  be  discussed  later,  they 
failed  to  observe  the  non-equilibrium  nature  of  soot  formation  because  their  temperatures  were 
too  high.  Interestingly  they  choose  not  to  measure  soot  formation  in  methane,  ethane,  or  methyl 
acetylene  flames  because  their  calculations  indicated  that  soot  would  be  formed  below  1000  K, 
and  thus  beyond  the  rich  flammability  limit. 

Gaydon  and  Wolfhard  [6]  in  1953  had  interpreted  a  1945  observation  by  Whittingham 
[10]  as  evidence  that  equilibrium  does  not  control  soot  formation.  Whittingham  had  observed 
that  as  little  as  0.1%  SO3  added  to  the  air  supply  of  a  coal-gas  Bunsen  flame  caused  a  non¬ 
sooting  flame  to  become  luminous.  "Non-equilibrium  soot  formation  in  premixed  flames"  was  the 
subject  of  a  paper  by  Millikan  [11]  in  1962.  He  measured  soot,  acetylene,  and  methane 
concentrations  in  a  premixed  ethylene/air  flame  as  a  function  of  C/O  ratios  on  both  sides  of  the 
measured  critical  ratio  for  soot  formation,  C/O  =  0.62.  The  acetylene  concentration  increased 
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linearly  from  the  non-sooting  through  the  sooting  flame  and  the  methane  concentration  increased 
slightly.  The  temperatures  were  between  1800  and  1820  K.  As  discussed  below,  equilibrium  soot 
production  is  zero  at  these  equivalence  ratios.  He  also  measured  a  change  with  temperature  of 
the  critical  C/0  ratio,  0.58  to  0.63.  The  C/O  ratio  for  soot  formation  increased  linearly  with 
temperature  between  1700  and  1820  K.  He  deduced  an  activation  energy  of  142  ±  40  kJ/mol 
from  his  results,  which  he  interpreted  as  the  difference  between  the  activation  energy  for  soot 
formation  and  soot  oxidation.  This  has  been  generally  accepted  as  the  explanation  of  the  high 
temperature  side  of  the  bell  shaped  curve  for  soot  formation  [4,12]. 

Haynes  and  Wagner  [4]  discuss  why  the  critical  C/O  ratio  for  soot  formation  occurs  more 
in  the  vicinity  of  0.5  than  1.0  where  thermochemistiy  would  predict  it  to  occur.  They  observe 
that  departure  from  equilibrium  is  also  shown  by  the  significant  quantities  of  CO2  and  H2O 
which  are  observed  well  beyond  the  critical  C/O  ratio  of  1.0  where  their  concentrations  should  be 
very  low.  They  explain  this  on  the  basis  of  parallel  reaction  paths  for  production  of  CO2  and 
H2O,  and  for  production  of  hydrocarbons.  Oxygen  rapidly  forms  H2O  and  CO2  and  is  locked  in 
this  form  for  the  short  reaction  time  available  in  the  flame  because  equilibrium  with 
hydrocarbons  are  via  hydrogen  atom  reactions  which  have  high  activation  energies  and  are  thus 
slow.  Oxygen  is  thus  no  longer  available  to  react  with  the  hydrocarbons  so  the  effective  oxygen 
concentration  is  reduced  and  thus  the  effective  C/O  ratio  is  increased.  This  allows  hydrocarbon 
radicals  to  continue  to  grow  toward  soot  while  preserving  their  radical  character. 

The  non-equilibrium  nature  of  soot  formation  is  not  discussed  in  the  1965  classic  review 
by  Palmer  and  Cullis  [14]  nor  in  several  other  more  recent  reviews  on  soot  formation  in  flames 
[12,15-17].  In  this  paper  experimentally  measured  mole  fractions  of  soot  and  acetylene  are 
compared  with  equilibrium  calculations  as  a  function  of  equivalence  ratio. 


SOURCE  OF  DATA 

The  experimental  soot  and  acetylene  yields,  and  flame  temperatures  as  a  function  of 
equivalence  ratio  were  taken  from  the  literature.  This  is  not  the  most  desirable  means  of 
collecting  data  but  the  sample  is  representative  in  that  a  number  of  fuels  and  conditions  are 
covered. 

Equilibrium  soot  and  acetylene  yields  were  calculated  using  an  equilibrium  program 
developed  by  Curt  Selph  and  Robert  Hall  at  Edwards  Air  Force  Base,  California.  This  program 
handles  both  solid  and  liquid  products  and  includes  in  the  calculation  all  the  species  in  the 
JANAF  Tables.  It  is  thus  unnecessary  to  chose  the  products,  which  can  often  lead  to  errors;  the 
program  automatically  eliminates  all  products  with  mole  fractions  below  some  critical  value.  We 
have  made  corrections  to  the  thermodynamic  database  where  warranted  and  have  added  many 
additional  species,  eg  large  aromatic  hydrocarbons  and  large  cations  and  anions  that  might  be 
involved  in  soot  formation.  These  include:  coronene,  C24HJ2)  ovalene,  C52HJ4,  circumcoronene, 
C54H^g  and  buckminsterfullerene,  C^q  ;  large  cations,  eg  perinaphthenyl  cation  (C13H9+)  and 
protonated  forms  of  cyclopentadienyl(l,10)-benz(8,9)-pyrene  (C2iHi^)+,  circumcoronene 
(C54H19+)  and  circumcircumcoronene  (C9(5H25  +  );  large  anions,  eg,  phenylacetyl  (CgHg-), 
perinaphthenyl  (Ci3H9-),  and  dibenz(a,h)-anthracene  (C22Hi4-)  anions.  The  program  calculates 


B-3 


TP-524 


mole  fractions  in  which  the  mole  fraction  of  a  solid  or  liquid  substance  is  the  molar  number  of 
atomic  units  in  the  solid  or  liquid  phase.  Soot  differs  from  graphite  in  that  it  is  about  10  mol 
percent  hydrogen.  We  were  not  aware  of  any  experimental  measurements  on  the 
thermochemistry  of  soot  so  we  estimated  the  thermochemical  properties  of  a  typical  soot  and 
included  that  in  the  equilibrium  calculations. 

To  estimate  the  thermochemical  properties  of  soot,  we  examined  the  C/H  ratios  for  the 
fully  condensed  benzene  ring  compounds  of  the  circumcoronene  series  At  n  =  10, 

the  molecular  formula  is  C^^qqHj^q  with  about  the  same  C/H  ratio  as  soot.  The  gas  phase 
quantities:  AHf°(298  K),  S°(298  K)  and  Cp(T)),  were  derived  using  Benson's  group  additivity 
method  [18,19].  To  reduce  the  gas  phase  data  to  solid  phase  "soot",  the  enthalpy  and  entropy  of 
sublimation  were  obtained  from  the  boiling  point  of  C^^qqHj^q  This  was  estimated  by  Joback's 
method  [20]  modified  to  fit  the  experimental  boiling  points  for  benzene,  naphthalene, 
phenanthrene,  and  coronene.  The  atmospheric  pressure  boiling  point  for  C(3qqH(5q  was  thus 
estimated  to  be  11,000  K.  The  entropy  of  sublimation  for  CqqqH^q  was  then  estimated  from  the 
boiling  point  by  Vetere's  relationship,  and  the  enthalpy  change  was  calculated.  These  values  were 
subtracted  from  the  gas  phase  values  to  give  the  enthalpy  of  formation  and  entropy  of  solid 
^600^60-  capacities  were  approximated  as  those  of  the  gas  phase  species  allowing  the 

entropy  and  enthalpy  as  a  function  of  temperature  to  be  calculated.  The  data  for  C^ogH^Q 
(pseudo-soot)  per  carbon  atom  and  graphite  are  compared  in  Table  1. 

When  the  derived  thermochemical  data  for  "soot"  replaced  graphite  in  the  equilibrium 
calculations  there  was  no  significant  difference  in  calculated  soot  or  acetylene  yields  nor  in  the 
threshold  equivalence  ratio  for  soot  appearance.  The  threshold  equivalence  ratio  depends  more 
on  stoichiometry  than  on  the  thermochemistry  of  the  soot. 

RESULTS  AND  DISCUSSION 

Acetylene  and  soot  mole  fractions  for  a  number  of  fuel  systems  are  shown  in  Figs.  1  to  5. 
The  plotted  experimental  values  are  ma.xima.  Above  C/O  =  1,  increasing  the  fuel  mole  fraction 
leads  to  increased  equilibrium  yields  of  soot,  due  to  insufficient  oxygen  to  consume  the  carbon. 
Inclusion  of  large  PCAH,  cations  and  anions  in  the  thermochemical  database  had  negligible 
effect  on  the  calculated  equilibrium  curves.  Equilibrium  mole  fractions  of  these  species  were 
below  the  sensitivity  limit  of  the  calculation,  about  10'^^  mole  fraction.  Except  for  fuel  acetylene, 
there  is  a  dramatic  difference  between  the  experimentally  measured  and  the  equilibrium 
calculated  mole  fractions  of  acetylene  and  soot  as  well  as  the  C/O  ratio  at  which  they  appear. 
The  critical  C/O  ratio  is  much  lower  than  predicted  by  equilibrium.  The  ratio  of  observed  mole 
fractions  of  acetylene  to  maximum  equilibrium  mole  fractions  is  especially  dramatic,  varying,  from 
a  factor  of  10^  for  fuel  acetylene  to  10^  for  methane.  The  difference  between  the  adiabatic  flame 
temperatures  and  the  measured  temperatures  are  due  to  the  use  of  flat  flame,  frequently  cooled, 
burners.  Experimental  flame  temperatures  are  graphed  when  available  for  a  range  of 
equivalence  ratios,  otherwise  they  are  simply  indicated  in  the  captions. 

The  results  for  a  methane/oxygen  flame  at  100  kPa  (1  atm.)  are  presented  in  Fig.  1.  The 
maximum  equilibrium  acetylene  yield  is  on  the  order  of  10'^  mole  fraction  and  peaks  at  C/O  = 
0.93.  Equilibrium  soot  formation  is  about  0.1  mole  fraction.  By  contrast,  the  experimental 
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observations  are  in  flames  at  a  much  lower  equivalence  ratios  but  still  fuel  rich,  equivalence  ratio 

>  1.  The  maximum  equivalence  ratio  at  which  either  soot  or  acetylene  are  measured  is 
determined  by  the  fuel  rich  flame  stability  limit.  In  flames  at  C/O  about  0.5,  the  experimental 
acetylene  mole  fraction  is  of  order  10  ^  which  is  about  lO”^  times  greater  than  the  maximum 
equilibrium  mole  fraction,  and  about  10^^  times  greater  than  the  equilibrium  mole  fraction  at 
C/O  =  0.5.  The  measured  soot  mole  fractions  are  comparable  to  the  equilibrium  values  for  C/O 

>  1  but  occur  at  much  smaller  C/O  ratios.  No  equilibrium  soot  production  is  calculated  at  that 
equivalence  ratio  at  which  soot  is  observed.  Clearly  soot  mole  fractions  in  these  methane/oxygen 
flames  are  far  greater  than  equilibrium  predictions. 

A  similar  situation  is  observed  for  other  relatively  low  temperature  flames: 
propane/oxygen  flames  at  15  kPa,  Fig.  2,  ethylene/oxygen  flames  at  100  kPa,  Fig.  3,  and 
toluene/air  flames  at  100  kPa  [34]  (not  shown).  Benzene/oxygen  flame  measurements  at  5.3  kPa, 
Fig.  4,  show  the  same  non-equilibrium  formation  of  acetylene  and  soot  but  at  C/O  ratios  much 
closer  (than  for  the  other  flames)  to  the  ratio  at  which  soot  would  be  produced  in  equilibrium  . 
The  measured  flame  temperature  is  also  greater,  about  2,000  K.  At  C/O  about  0.95  the 
equilibrium  acetylene  mole  fraction  is  about  10  mole  fraction,  so  that  the  measured  mole 
fraction  is  about  10^^  greater  than  equilibrium.  This  flame  is  also  distinguished  from  the  others 
in  that  the  equilibrium  acetylene  mole  fraction  for  C/O  >  1  is  orders  of  magnitude  greater,  due 
the  greater  flame  temperature. 

For  both  neutral  and  ion  growth,  increasing  molecular  weight,  normalized  by  the  number 
of  carbon  atoms,  is  a  thermodynamically  favored  process,  i.e.  the  system  approaches  an 
equilibrium  with  increasing  molecular  size  [35].  However,  in  the  flames  discussed  above,  soot  is  a 
supra-equilibrium  product.  This  is  clearly  a  chemical  kinetically  controlled  process.  It  thus 
requires  a  driving  force  from  a  non-equilibrium  precursor.  Something  must  be  present  in  excess 
of  equilibrium  to  initiate  the  process.  Two  kinetic  explanations  have  been  offered:  (1)  Haynes 
and  Wagner's  explanation  that  fast  reactions  leading  to  CO2  and  H2O  remove  oxygen  from  the 
system  giving  an  effective  greater  C/O  ratio  [4]  producing,  e.g.,  supra-equilibrium  acetylene  mole 
fractions;  and  (2)  supra-equilibrium  ion  or  hydrogen  atom  mole  fractions  [36,37]  drive  the  growth 
reactions  to  a  pseudo-equilibrium.  We  give  evidence  below  that  a  combination  of  these  two 
mechanisms  may  be  operative. 

To  test  the  Haynes  and  Wagner  concept,  an  effective  C/O  ratio  was  estimated  for  a 
methane/oxygen  flame  in  Fig.  1  by  assuming  the  experimental  quantity  of  oxygen  in  CO2,  CO  and 
H2O  in  excess  of  the  equilibrium  mole  fractions  was  removed  from  the  system  and  a  C/O  ratio 
was  calculated  for  the  remaining  gases.  An  experimental  C/O  =  0.635  was  thus  increased  to 
about  0.9.  This  is  in  the  right  direction.  Unfortunately  this  procedure  does  not  differentiate 
from  a  mechanism  where  supra-equilibrium  atomic  hydrogen  or  ion  concentrations  with  supra- 
equilibrium  acetylene  concentration  drive  the  process  toward  soot.  On  the  other  hand  observed 
supra-equilibrium  concentrations  of  acetylene  can  not  be  accounted  for  by  supra-equilibrium 
concentrations  of  atomic  hydrogen  or  ions.  Thus  the  Haynes  and  Wagner  mechanism  must  play 
a  role  in  supra-equilibrium  soot  formation. 

For  the  free  radical  mechanism,  the  driving  force  is  through  supra-equilibrium  hydrogen 
atom  concentrations  [36].  For  the  ionic  mechanism  it  is  through  chemiions,  formed  by 
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chemiionization,  a  well  established  nonequilibrium  process  which  produces  ions  in  great  excess  of 
equilibrium  [37],  The  importance  of  supra-equilibrium  acetylene,  and  presumably  diacetylene 
and  other  hydrocarbons  is  indicated  by  the  results  reported  in  Figs.  1  to  4.  Supra-acetylene 
concentrations  must  be  considered  as  a  main  driving  force  for  soot  formation  and  not  simply  a 
building  block.  Hydrogen  atoms  or  chemiions  are  necessary  because  it  has  been  demonstrated 
that  the  mechanism  does  not  involve  polymerization  of  acetylene  alone  [38];  otherwise 
polyacetylenes  would  be  the  precursor  species.  It  is  generally  agreed  that  the  precursor  species  is 
a  polycyclic  compound,  either  ionic  or  neutral.  Much  more  acetylene  is  consumed  in  the 
nucleation  steps  as  well  as  in  the  growth  steps  than  ions  or  hydrogen  atoms.  Only  one  ion  is 
required  per  soot  particle.  Thus,  both  ions  (or  hydrogen  atoms)  and  acetylene  must  be  present, 
and  this  occurs  at  an  equivalence  ratio  where  the  ion  concentration  is  falling  and  the  acetylene 
concentration  is  rising.  This  merits  a  more  quantitative  treatment  which  we  plan  to  pursue. 

Results  for  the  well  studied  acetylene/oxygen  flame  at  2.67  kPa  and  50  cm/s  unburned  gas 
flow  on  a  flat  flame  burner  are  shown  in  Fig.  5.  The  equilibrium  mole  fractions  are  calculated 
for  the  experimental  temperatures.  These  results  are  completely  different  from  the  flames 
already  discussed.  Like  the  benzene  flames,  the  measured  flame  temperatures  are  close  to  2000 
K,  yet  soot  is  observed  only  where  equilibrium  indicates  it  should  be  observed.  Further,  soot  is 
observed  at  a  smaller  mole  fraction  than  predicted  by  equilibrium.  On  the  other  hand,  the 
acetylene  mole  fraction  is  considerably  in  e.xcess  of  that  observed  in  the  other  systems.  The 
flame  temperatures  above  C/O  =  1  remain  high  because  acetylene  is  a  highly  endothermic 
compound. 

Sub-equilibrium  concentrations  of  soot  in  the  acetylene  flame  can  be  explained  by  two 
effects  of  the  high  flame  temperature;  (1)  the  high  temperature  increases  the  general  rate  of 
reactions  toward  equilibrium  and  (2)  entropy  becomes  important  in  the  growth  of  large  molecular 
species  at  high  temperature.  Any  system  with  a  decrease  in  the  number  of  species  has  a  negative 
entropy.  In  soot  formation  many  small  molecules  combine  to  make  large  molecules  and  then 
particles.  Entropy  effects  become  increasingly  important  as  the  temperature  increases: 

AG(r)  =  AH(r)  -  T  AS(r)  (3) 

For  reactions  producing  fewer  species,  at  some  temperature,  TAS(r)  becomes  more  important 
than  AH(r)  and  the  reaction  becomes  less  favorable,  and  may  reverse  direction.  We  suggest  that 
in  the  acetylene  flame  the  elementary  kinetic  steps  in  the  reactions  leading  to  equilibrium  soot 
formation  are  slowed  down  in  the  forward  direction.  In  the  other  flames  above,  except  the 
benzene  flame,  the  temperatures  are  lower  so  that  AG(r)  for  the  elementary  steps  is  dominated 
by  AH(r)  and  the  reactions  are  driven  toward  soot. 

It  is  unfortunate  that  the  standard  low  pressure  acetylene/oxygen  flame  which  has 
received  so  much  attention  behaves  so  differently  from  other  hydrocarbons  flames.  The  question 
is  raised,  do  arguments  for  the  mechanism  of  soot  formation,  eg  the  ionic  mechanism,  derived 
from  studies  of  this  flame  apply  to  soot  formation  in  other  flames.  This  question  deserves 
further  attention. 
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IV.  SUMMARY 

Experimental  mole  fractions  of  acetylene  in  all  flames  studied  are  in  great  excess  of 
equilibrium.  The  observed  soot  mole  fractions  also  exceed  the  equilibrium  values  in  all  the 
flames  except  for  the  acetylene/oxygen  flames.  The  differences  sometimes  exceed  10  .  There 
are  two  explanations  which  are  simultaneously  operative.  In  one,  due  to  Haynes  and  Wagner, 
oxygen  is  tied  up  in  water  and  carbon  dioxide  which  are  slow  to  attain  equilibrium  with  the  rest 
of  the  system  so  that  the  growing  hydrocarbons  are  in  an  effective  C/O  mixture  >1.0  in  which 
equilibrium  favors  soot  formation.  In  the  other  mechanism,  the  initial  shuffle  of  reactions  in  the 
flame  front,  as  the  system  goes  from  a  highly  non-equilibrium  reactant  mixture  to  equilibrium, 
acetylene  and  some  other  hydrocarbons  are  produced  in  supra-equilibrium  concentrations. 
Chemiionization  and  formation  of  excess  hydrogen  atoms  are  a  part  of  this  process.  These  drive 
the  reactions  to  produce  soot  in  pseudo-equilibrium  concentrations.  The  acetylene/oxygen  flame 
has  a  greater  flame  temperature  so  that  general  equilibrium  is  more  nearly  reached.  At  the  same 
time  entropy  effects  slow  down  the  rate  of  individual  reactions  leading  to  soot  so  soot  never 
reaches  the  equilibrium  concentration.  The  formation  of  acetylene  in  supra-equilibrium 
concentrations  plays  a  major  role  in  the  supra-equilibrium  formation  of  soot  beyond  that  of 
simply  acting  as  a  building  block. 
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TABLE  1 


THERMOCHEMISTRY  QUANTITIES  FOR 
GRAPHITE  AND  SOOT  (AS  C^^ooH^q) 


GRAPHITE 

SOOT 

AH°£,  kJ/mol 

0.00 

5.95 

S°,  J/mol  K 

Cp,  J/mol  K 

5.80 

8.86 

298 

8.54 

9.22 

500 

14.5 

17.1 

1000 

21.6 

23.8 

1500 

23.9 

25.5 

2000 

25.1 

26.3 

2500 

25.9 

26.7 
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FIGURE  1  EXPERIMENT  VS.  EQUILIBRIUM 
METHANE/OXYGEN  FLAMES,  100  kPa 

Equilibrium  values;  broken  lines,  calculated  for  adiabatic  flames. 
Experimental  results:  C2H2:  A  -  D'Allesio,  et  al  (22);  Soot:  ♦  -  AeroChem; 

•  -  D'AIessio,  et  al  (22);  ■  -  Flower  (24);  Temperature:  O,  D'AIessio,  et  al  (22,25). 


FIGURE  2  EXPERIMENT  VS.  EQUILIBRIUM 
ETHYLENE/AIR  FLAMES,  100  kPa. 

Equilibrium  values:  broken  lines,  calculated  for  adiabatic  flames. 
Experimental  results:  C2H2:  D  -  Millikan  (11);  O  -  Harris  &  Weiner  (26); 

Soot:  •  -  Harris  &  Weiner  (26);  A  -  Bonig,  et  al  (27); 
Experimental  flame  temperatures:  Scattered.  1650  -  1820  K. 


FIGURE  3  EXPERIMENT  VS.  EQUILIBRIUM 
PROPANE/OXYGEN  FLAMES,  15  kPa 

Equilibrium  values:  broken  lines,  calculated  for  adiabatic  flames. 

Experimental  results,  Bockhorn,  et  al  (28):  A  -  C2H2;  ■  -  Soot.  Temperature:  1800  K. 


FIGURE  4  EXPERIMENT  VS.  EQUILIBRIUM 
BENZENE/10%  ARGON/OXYGEN  FLAMES,  5.3  kPa. 

Equilibrium  values:  broken  lines,  calculated  for  adiabatic  flames. 
Experimental  results,  McKinnon  (29):  A  -  C2H2'>  ®  '  Soot;  O  -  Temperature. 


FIGURE  5  EXPERIMENT  VS.  EQUILIBRIUM 
ACETYLENE/OXYGEN  FLAMES,  2.7  kPa. 

Equilibrium  values:  broken  lines,  calculated  for  adiabatic  flames. 
Experimental  results;  C2H2:  O  -  Homann,  et  al  (30),  Bonne,  et  al  (31); 
Soot:  •  -  Homann,  et  al  (30),  Bonne,  et  al  (31);  a  -  Wersborg,  et  al  (32); 

Temperature:  □  -  AeroChem. 
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REVIEWER  A,  No.  081 

I  cannot  recommend  acceptance  of  this  paper  for  the  Combustion  Symposium  in  any  form. 

The  authors  resurrect  themodynamic  analysis  for  soot  formation,  whereas  it  is  well  established  1 
and  well-accepted  that  soot  formation  is  a  kinetically  controlled  process.  Thus,  this  paper's 
"findings"  add  little  to  the  current  knowledge. 

It  is  not  clear  what  authors'  motivation  was  for  this  work.  In  the  opening  sentence  of  the  2 

Introduction  they  state  that  their  motivation  was  "to  understand  the  ionic  mechanism  of  soot 
formation."  Forgetting  for  a  second  that  the  ionic  mechanism  has  been  rejected  by  the  research  3 
community  working  on  soot  formation,  I  find  nothing  in  this  paper  that  addresses  the  ionic 
mechanism  of  soot  formation.  Look  for  example  at  the  Summary:  there  is  nothing  there  about 
the  ionic  mechanism.  There  is  also  nothing  on  this  in  the  Discussion. 

The  arguments  presented  in  the  manuscript  are  not  clear,  if  not  incorrect.  How  can  one  deduce  ^ 
from  a  thermodynamic  analysis  an  underlying  mechanism?  Thermodynamics  describes  state 
properties,  i.e.,  those  independent  of  reaction  path. 

The  analysis  itself  is  not  performed  consistently.  The  final  thermodynamic  state  is  strongly  j_ 

dependent  on  temperature.  The  authors  assumed  an  adiabatic  flame  temperature  for  some  of  the  ^ 

conditions  and  the  experimental  one  for  others.  Using  adiabatic  flame  temperature  is  certainly 
not  a  correct  assumption,  which  may  explain  some  of  the  "disagreements"  in  the  results. 

So,  what  are  the  accomplishments  of  this  work?  A  thermodynamic  analysis  (performed 
inconsistently)  does  not  explain  the  experiment.  The  authors  suggest  that  something  must  be  in  6 
disequilibrium,  the  conclusion  that  no  one  doubts.  And  some  of  the  suggestions  do  not  make  any 
physical  sense;  e.g.,  the  last  sentence;  "The  formation  of  acetylene  in  supra-equilibrium 
concentrations  plays  a  major  role  in  the  supra-equilibrium  formation  of  soot  beyond  that  of 
simply  acting  as  a  building  block."  Then,  what  is  this  role? 

Some  minor  points: 

The  authors  made  an  attempt  to  provide  a  historical  review  of  thermodynamic  analysis  of  soot 
formation.  Yet,  they  left  out  two  important  studies:  that  of  Farmer  who  suggested  to  use  C96H24  7 

for  a  soot  model,  very  much  what  the  present  authors  did,  and  the  work  of  AJberty,  who  provided  0 
the  most  elegant  thermodynamic  analysis  of  PAH  formation. 

Many  questions  that  the  authors  pose  in  this  manuscript  have  been  answered  by  others  in  terms 
of  non-ionic  mechanisms,  to  which  the  present  authors  do  not  even  make  a  reference.  ^ 
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This  work  is  one  of  the  first  attempts  to  model  ion  concentrations  in  a  sooting  flame  for 
which  the  20  Torr  acetylene/oxygen  flat  flame  with  C/0  =  1 .2  has  been  chosen.  Of  course,  there 
are  still  a  great  number  of  shortcomings,  but  this  is  a  beginning  and  one  cannot  learn  how  to 
improve  things  if  one  does  not  try  them  out.  What  disturbs  me  most  with  this  paper  is  that  it 
seems  to  be  an  'extended  abstract'  of  (a)  far  more  detailed  paper(s)  in  preparation  [11,  17,  18]. 

Those  seem  to  contain  some  of  the  really  interesting  questions,  assumptions  and  conditions  used 
in  the  calculations,  for  example 

How  well  does  the  Sandia  flame  code  and  the  applied  total  mechanism  and  rate  constants  ^ 
describe  the  neutral  structure  of  this  flame  which  must  be  counted  to  the  more  heavily 
sooting  flames  (threshold  C/0  =  0.85). 

Why  have  the  calculations  been  applied  to  such  a  rich  flame  and  not  performed  for  a  2 
sooting  flame  nearer  at  the  soot  threshold  ? 

Which  of  the  experimental  profiles  of  the  neutrals  have  been  put  into  the  calculation  ?  3 

What  are  the  parameters  of  equation  (2)  used  for  the  different  ion-molecule  reactions  ?  4 

How  have  the  well-established  equilibria  between  the  protonated  polyynes  been  handled,  ^ 
or,  if  they  were  disregarded,  how  does  the  calculated  profile  for  C5H3  compare  with  the 
equilibrium  concentrations  ? 

How  have  the  proton  transfer  reactions  been  handled  which  belong  to  the  fastest  ion-  6 
molecule  reactions  ? 

It  does  not  become  clear  how  electro-neutrality  has  been  accounted  for.  I  understand  that 
the  experimental  profiles  of  C3H3*  and  CgH5*  have  been  put  into  the  calculation.  On  the  7 
other  hand,  electron-ion  combinations  are  included  which  cause  a  loss  of  these  species. 
However,  this  does  not  seem  to  be  accounted  for,  since  the  experimental  profiles  of  these 
two  kinds  of  ions  are  used  throughout  the  flame.  The  electrons  are  consumed  by  the 
reactions  and  put  again  into  the  calculation  to  maintain  electro-neutrality,  into  which  also 
enter  the  two  preset  profiles  ?  How  then  is  internal  consistency  maintained  ? 

Does  the  calculation  indicate  steady-state  conditions  for  any  kind  of  ions,  and  if  so,  for  8 
which  ? 

This  list  of  open  questions  is  not  complete.  It  is  very  difficult  for  the  reviewer  to  come  to  a 
conclusion,  whether  this  manuscript  should  be  recommended  for  the  Symposium.  On  the  one 
hand,  this  is  an  interesting  piece  of  work  which  continues  with  the  efforts  to  model  flames  to  the  0 
best  of  our  knowledge.  On  the  other  hand,  symposium  papers  should  not  be  'devaluated'  to 
extended  abstracts  of  papers  published  somewhere  else.  I  must  leave  this  decision  to  the  Program 
Committee.  If  this  kind  of  manuscript  is  generally  accepted,  it  ranks  under  "good",  if  the  authors 
can  address  at  least  some  of  the  problems. 
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RESPONSE  TO  REVIEWERS'  COMMENTS 


The  reviewer  is  identified  by  letter  on  his  review  copy  and  the  reviewer's  comment  is 
identified  by  a  number  on  the  right  of  the  review.  References  numbers  refer  to 
references  in  subject  paper. 


Reviewer  A,  No.  081 

1.  It  would  be  useful  if  the  reviewer  would  furnish  at  least  one  reference  which 
e.g.  shows  that  is  it  well  known  that  in  most  flames  kinetics  controls  the 
production  of  soot  in  excess  of  equilibrium,  and  that  in  the  standard  flame  soot 
production  is  less  than  equilibrium,  or  the  extent  to  which  soot  production 
differs  from  equilibrium. 

2.  See  paragraph  top  pg.  2  and  continual  comparison  of  H  atom  or  ions  as  driving 
force  for  soot  formation. 

3.  We  are  not  aware  that  the  Board  of  Directors  of  The  Combustion  Institute  or 
the  National  Academy  of  Sciences  or  any  such  authoritative  body  has  made 
such  a  decision. 

4.  Done  all  the  time,  a  kinetic  mechanism  must  be  consistent  with 
thermochemistry.  In  this  case  for  most  of  the  flames,  our  consideration  of 
thermochemistry  has  shown  that  non-equilibrium,  e.g.,  supra-equilibrium, 
concentration  driven  reactions  lead  to  soot  formation,  but  in  the 
acetylene/oxygen  flame  soot  yields  are  sub-equilibrium,  indicating  kinetic 
limitations. 

5.  Adiabatic  flame  temperatures  were  only  used  where  experimental  temperatures 
were  unavailable.  This  does  not  account  for  any  of  the  observations,  or 
"disagreements". 

6.  Acetylene  is  frequently  considered  a  "building  block"  for  soot  formation,  ie  the 
addition  of  acetylene  to  the  growing  PCAH,  neutral  or  ion,  is  the  main  source 
of  carbon.  This  could  occur  whether  the  acetylene  is  present  in  equilibrium  or 
in  excess,  so  long  as  there  is  sufficient  concentration.  Its  presence  in  supra- 
concentration  means  it  could  contribute  to  the  production  of  supra-equilibrium 
concentrations  of  soot. 

7.  Farmer  [Wang,  Matula  and  Farmer,  18th  Combustion  Symposium,  p  1 149]  do 
not  propose  "C96H24  for  a  soot  model"  but  as  a  precursor.  Their  choice  was 
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based  on  thermodynamics  of  formation  and  that  this  particular  molecule, 
circumcircumcoronene  is  stable  in  their  equilibrium  calculation  at  about  the 
temperature  soot  is  formed.  They  in  fact  do  not  subscribe  to  the  reviewer's 
statement  (1 ).  See  e.g.,  pg.  1154"---  and  that  if  C96H24  does  not  tend  to  be 
formed  (as  indicated  by  equilibrium  calculations  )  soot  will  not  be  formed." 

8.  Agreed  that  Alberty  "provided  the  most  elegant  thermodynamic  analysis  of  PAH 
formation".  We  have  used  his  work  and  many  others  in  developing  the 
thermodynamics  used  in  this  work,  but  he  did  not  estimate  the  thermodynamics 
of  soot. 

9.  Unfortunately,  we  are  limited  to  5,500  words  and  found  it  impossible  to 
reference  all  prior  works;  there  are  already  38  references. 


Reviewer  B.  No.  428 

1 .  The  reason  Reviewer  B  cannot  reproduce  the  equilibrium  calculations  should  be 
clear  from  the  following  responses. 

2.  "Driving  force"  is  not  a  qualitative  concept;  and  is  a  commonly  used  term.  The 
quantitative  aspect  in  this  paper  is  the  ratio  of  measured  to  equilibrium 
concentrations. 

3.  It  is  suggested  that  the  differences  we  express  in  concentration  could  be  better 
expressed  as  "thermodynamic  affinities".  One  of  the  authors  (H.F.  Calcote)  has 
been  involved  in  studying,  and  of  course  reading,  non-equilibrium  chemistry  for 
47  years  and  does  not  recall  one  author  discussing  combustion  non-equilibrium 
in  terms  of  "thermodynamic  affinities".  This  term  serves  a  completely  different 
need  than  describing  non-equilibrium. 

4.  The  absolute  H  atom  concentration  is  greater  than  the  ion  concentration  but  the 
ratio  of  measured  to  equilibrium  concentration  is  greater  for  ions.  The  driving 
force  is  not  the  concentration  but  the  ratio  of  measured  to  equilibrium 
concentrations. 

5.  Complex  equilibria  are  not  calculated  as  a  series  of  individual  equilibrium 
reactions  as  assumed  by  the  Reviewer.  That  would  be  far  too  complicated 
when  such  a  large  number  of  species  are  involved.  Computer  programs  have 
been  developed  which  do  this  on  the  basis  of  heats  of  formation  of  reactants 
and  free  energies  of  formation  of  all  possible  products  as  a  function  of 
temperature.  Any  combination  of  reactants  and  products  which  can  be  written 
will  of  course  be  in  equilibrium  and  can  be  calculated  based  on  the  final 
temperature  and  appropriate  equilibrium  calculated  species. 
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6. 


Perhaps  we  should  have  said  the  only  property  shown  to  dramatically  change 
at  soot  threshold  is  the  ionic  structure  of  the  flame.  It  would  be  helpful  if  the 
reviewer  would  give  references  to  the  "great  number  of  properties"  he  refers 
to,  we  are  unaware  of  them.  Some  of  the  few  such  data  we  are  aware  of  are 
shown  in  Figure  7  of  this  final  report.  The  dramatic  change  in  ionic  character 
are  compared  with  changes  in  concentrations  of  the  PCAH,  etc  (which  are 
important  in  the  neutral  models  of  soot  formation).  This  supports  our 
statement,  and  should  be  recognized  as  additional  evidence  for  the  ionic 
mechanism, 

7.  We  know  of  no  such  quantitative  concentration  measurements  of  anions  in 
sooting  flames. 

8.  Supraequilibrium  H  atom  concentrations  is  the  driving  force  for  soot  formation 
in  proposed  radical  mechanisms.  We  did  not  mean  to  imply  that  this  concept 
has  been  proven,  just  that  it  is  a  part  of  the  (unproven)  mechanism,  much  as 
the  role  of  ions  in  the  (unproven)  ionic  mechanism  of  soot  formation.  However, 
the  length  of  the  paper  was  limited  by  Symposium  rules,  and  we  felt  that  these 
concepts  are  generally  understood  in  the  soot  mechanism  community.  They  are 
also  described  in  the  given  references.  If  we  had  more  space  we  would 
dedicate  it  to  more  relevant  issues  involving  the  ionic  mechanism. 

9.  In  the  ionic  mechanism,  one  ion  grows  to  produce  one  soot  particle,  where 
many  acetylenes  are  required  to  produce  one  soot  particle. 
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Simulation  of  Electric  Field  Effects  in  Premixed  Methane  Flames 

TIMOTHY  PEDERSEN  and  ROBERT  C.  BROWN 

Depanment  of  Mechanical  Engineering.  Iowa  State  University,  Ames,  lA  50011 


The  objective  of  this  study  is  to  predia  the  effect  of  electric  fields  on  the  ionic  strucnire  of  one-dimensional 
methane  flames.  The  chemical  kinetic  mechanism  devised  for  this  study  combines  existing  methane  oxidation 
with  a  scries  of  chemiionization,  ion-molecule,  and  dissociative-recombination  reactions  thought 
to  be  important  to  the  ionic  structure  of  methane  flames.  The  governing  equations  include  Poisson’s  equation 
to  describe  the  electric  field  intensity  in  the  flames.  The  model  correctly  prediaed  several  features  in  the 
ionic  structure  of  methane  flames.  The  model  was  used  to  predict  saturation  currents  as  a  function  of  peak 
flame  temperature  from  which  the  apparent  activation  energy  for  the  saturation  process  was  determined.  The 
predicted  apparent  aaivation  energy  compared  weU  with  the  experimental  range  of  values.  Thw  analysis 
determined  that  the  rate-limiting  reaction  in  the  overall  ionic  reaction  mechanism  is:  C;H30*-i-  H,  -* 
HjO'-i-  CjH,. 


NOMENCLATURE 

A  burner  cross-sectional  area 

Cp  specific  heat  at  constant  pressure 

Z),  ambipolar  diffusivity 

Dij  binary  diffusivity  of  species  i  through  j 

Di_  „  mixture  diffusivity  of  species  i 

dI  t  mixture  thermal  diffusivity  of  species  i 

e  electron  charge 

E  electric  field  intensity 

Ep  activation  energy 

Eppp  apparent  activation  energy 

G  generalized  parameter  (see  Table  3) 

A/fy  heat  of  formation 

J,  mass  flux  of  species  i 

j  electric  current  flux 

saturation  current  flux 

Kt  mixture  thermal  diffusion  ratio  of  spe- 
cies  i 

k  rate  coefficient 

kf,  Boltzman’s  constant 

M  molecular  weight 

n  J.  positive  ion  molar  concentration 
rt  _  negative  charge  molar  concentration 

R  universal  gas  constant 

T  temperature 

Tf  peak  flame  temperature 
t  time 

u  velocity 

V  voltage 

X  mole  fraction 
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X  spatial  variable 

Y  mass  fraction 

permittivity  of  free  space 
r,  F'  generalized  parameters  (see  Table  3) 

O  equivalence  ratio 

6  generic  variable 

p  gas  density 

binary  mobility  of  species  i  through  spe¬ 
cies  j 

lij  „  multicomponent  mobility  of  species  i 
'V  transformed  coordinate  system  variable 

(o  species  production  rate  per  unit  volume. 

Subscripts 

e  electron 

i  species  designation 

j  species  designation 

nspc  number  of  species 

o  initial  value 


INTRODUCTION 

The  objective  of  this  study  is  to  predict  the 
effect  of  electric  fields  on  the  ionic  structure  of 
methane  flames.  An  understanding  of  this  phe¬ 
nomenon  could  be  imponant  in  attempts  to 
control  blowoff  limits,  flame  speed,  and  soot 
formation  by  the  application  of  electric  fields. 
The  source  of  flame  ions  is  generally  ac- 
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cepted  to  be  the  chemiionization  reaction 

CH  +  O  ^  CHO*+  e-, 

as  first  proposed  by  Calcote  [1],  Cool  and 
Tjossem  [2]  determined  that  the  reaction  of 
electronically  excited  CH  with  oxygen  atoms  to 
be  2000  times  faster  than  the  same  reaction 
with  ground-state  CH.  Accordingly,  the  reac¬ 
tion 

CH*(A-A)  +  O  CHO^+  e- 

may  also  be  a  strong  source  of  flame  ions,  a 
hypothesis  that  is  supported  by  the  computa¬ 
tional  study  of  Eraslan  and  Brown  [3]. 

The  CHC  ion.  is  not  the  dominant  ion  in 
hydrocarbon  flames.  Calcote  [4]  noted  that 
CHO"  would  be  quickly  consumed  by  HiO  to 
produce  via 

CHO^a  H,0  ^  CO  ^  H3O". 

The  H3  0~  ion  can  form  other  flame  ions  by 
a  series  of  ion-molecule  reactions  or  it  can 
undergo  dissociative  recombination  [4] 

H3O-A  e"- H3O  ^  H. 

A  large  ratio  of  electron  concentration  to 
negative-ion  concentration  has  been  found  in 
hydrocarbon  flames  [5].  The  electron,  with  its 
small  mass  and  high  mobility,  has  a  high  diffu¬ 
sion  rate  out  of  the  flame,  resulting  in  a  charge 
distribution  in  the  flame.  The  result  is  a  self-in¬ 
duced  electric  field  through  the  flame  that 
retards  the  charge  separation.  This  process  of 
ambipolar  diffusion  results  in  electrons  and 
positive  ions  diffusing  out  of  the  flame  front  at 
equal  rates. 

When  an  external  electric  field  is  imposed 
across  the  flame  front,  electrons  and  ions  are 
removed  at  a  rate  proponional  to  the  electric 
field  strength.  Eventually  a  field  strength  is 
reached  at  which  electrons  and  ions  are  re¬ 
moved  from  the  flame  front  as  fast  as  they  are 
created  by  chemiionization.  The  electric  cur¬ 
rent  flux  drawn  from  the  flame  under  this  limit 
condition  is  termed  the  saturation  current  flux, 

j  S' 

Lawton  and  Weinberg  [6]  appear  to  be  the 
first  to  report  a  comprehensive  study  of  satura¬ 
tion  currents.  They  measured  saturation  cur¬ 
rents  with  a  porous  disc  burner,  which  served 


as  the  cathode,  and  a  plane  electrode  mounted 
parallel  to  and  slightly  above  the  flame  to  serve 
as  the  anode.  The  peak  flame  temperature  7^ 
of  this  flat  flame  could  be  controlled  by  varying 
the  flow  rate.  The  researchers  carried  out  ex¬ 
periments  on  methane,  ethylene,  and  pro¬ 
pane-air  mixtures.  They  also  investigated  hy¬ 
drogen/  air  mixtures  to  which  trace  amounts 
of  hydrocarbon  were  added.  Saturation  cur¬ 
rents  as  a  function  of  final  flame  temperature 
were  not  documented  by  Lawton  and  Wein¬ 
berg  [6]  for  hydrocarbon  flames.  However,  the 
saturation  current  as  a  function  of  equivalence 
ratio  for  each  hydrocarbon  flame  was  docu¬ 
mented.  For  the  hydrocarbon  seeded  hydro¬ 
gen-air  flame,  the  variation  of  saturation 
currents  with  final  flame  temperature  was 
recorded.  The  hydrogen— air  flames  were 
seeded  with  up  to  1%  hydrocarbon. 

Lawton  and  Weinberg  [6]  plotted  ]n(/j)  ver- 
sus  1/7^  to  find  the  apparent  activation  energy 
^app  for  the  ionization  process.  Apparent  acti¬ 
vation  energies  reported  for  the  investigated 
flames  ranged  in  value  from  37.5  ±  2  to  67  ±  4 
kcal/ mol  depending  on  equivalence  ratios  and 
hydrocarbon  fuel.  These  researchers  recog¬ 
nized  that  the  apparent  activation  energies 
were  much  larger  than  the  activation  energies 
for  proposed  chemiionization  processes  in 
flames.  They  did  not  draw  any  firm  conclusions 
about  the  resulting  values  of  apparent  activa¬ 
tion  energies.  It  was  emphasized  by  the  re¬ 
searchers  that  the  work  presented  was  purely 
to  establish  a  practical  method  for  the  study  of 
ion  generation  rates  in  flames. 

In  1968  Peeters  and  Van  Tiggelen  [7]  mea¬ 
sured  the  rate  of  chemiionization  in  flames 
with  equivalence  ratios  dose  to  unity.  Their 
experimental  setup  included  a  Powling-Eger- 
ton  type  burner  open  to  the  atmosphere.  The 
burner  served  as  the  cathode  and  a  conducting 
plate  placed  above  the  flame  served  as  the 
anode.  They  measured  the  saturation  current 
for  various  flame  temperatures.  The  apparent 
activation  energy  was  found  to  be  73  kcal/mol. 
However,  no  attempt  was  made  to  relate  this 
activation  energy  to  a  detailed  reaction  mecha¬ 
nism. 

About  a  decade  later,  Bertrand  et  al.  [8] 
measured  the  burning  velocity  of  a  stoichio¬ 
metric  methane  flame  stabilized  on  a  flat  flame 
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burner.  The  researchers  measured  the  varia¬ 
tion  of  saturation  current  with  final  flame  tem¬ 
perature.  The  apparent  activation  energy  for 
this  flame  was  equal  to  53  kcal/mol.  Again,  no 
attempt  was  made  to  relate  the  apparent  acti¬ 
vation  energy  to  a  reaction  mechanism. 

Numerical  methods  used  in  flame  simula¬ 
tions  have  progressed  to  the  point  where  it  is 
now  possible  to  predict  saturation  currents  in 
flames.  The  focus  of  this  research  is  to  relate 
the  apparent  activation  energies  obtained  from 
the  predicted  saturation  currents  to  a  chemical 
kinetic  reaction  mechanism.  We  choose  to  in¬ 
vestigate  methane  flames  with  equivalence  ra¬ 
tios  ranging  from  0.2  to  2.13  since  considerable 
experimental  data  exist  for  these  flames  and 
the  (nonionic)  chemical  kinetic  mechanism  is 
relatively  well  understood. 

CHEMICAL  KINETIC  MECHANISM 

The  computational  model  requires  a  chemical 
kinetic  mechanism  that  includes  both  neutral 
species  and  ionic  species  reactions.  Our  study 
e.xamines  saturation  currents  for  both  methane 
flames  and  methane-seeded  hydrogen  flames; 
accordingly,  we  have  attempted  to  devise  a 
neutral-species  kinetic  mechanism  that  can 
handle  both  types  of  flames.  The  ionic  mecha¬ 
nism  must  allow  for  the  production  and  de¬ 
pletion  of  ions  through  chemiionization,  ion- 
molecule  reactions,  and  dissociative  recombi¬ 
nation.  Only  major  ions  are  considered  in  the 
ionic  mechanism  since  they  constitute  more 
than  85%  of  the  total  ion  current.  The  com¬ 
plete  chemical  kinetic  mechanism  is  listed  in 
Table  1. 

Our  neutral  species  mechanism  draws  heav¬ 
ily  upon  the  methane  oxidation  mechanism  of 
Coffee  [9]  who  assembled  data  from  several 
sources.  The  first  21  reactions  were  taken  from 
Dixon-Lewis  [10].  These  were  verified  in  stud¬ 
ies  performed  on  Hj/O^/N,  and  CO/H,/ 
O^/ N2  flames  [11].  The  inclusion  of  these 
reactions  make  it  possible  to  study  the  conse¬ 
quences  of  seeding  hydrogen  flames  with 
methane. 

Dean  et  al.  [12,  13]  performed  a  series  of 
shock  tube  experiments  that  are  the  basis  for 
including  reactions  22-29.  The  first  four  reac¬ 
tions  of  this  set  allow  for  the  formation  of  CHj 


from  CH4.  Formaldehyde  (CHiO)  is  then 
formed  from  CHj.  The  reactions  that  include 
CHO  (reactions  30-33)  have  reaction  rates 
that  are  not  well  known.  The  values  from 
Dixon-Lewis  [10]  are  used  for  the  lack  of  bet¬ 
ter  data.  Reactions  34-37  include  the  oxida¬ 
tion  of  CH4.  This  set  of  reactions  were  shown 
by  Gelinas  [14]  to  be  imponant  in  modeling 
CH4  flames. 

The  primary  source  of  C,  hydrocarbons  is 
reaction  38.  The  rate  for  this  reaction  was 
taken  from  Wamatz  [15]  and  Dixon-Lewis  [10], 
which  was  then  modified  for  atmospheric  con¬ 
ditions  by  Luther  and  Troe  [16].  Reactions 
38-56  provide  the  sequence  by  which  is 

converted  to  CHj  and  CO: 

C;Hj  ^  C,H5  ^  C;H4  -»  C2H3 

-  C^H,  -  CHj.CO. 

Also,  reactions  46  and  47  form  a  pathway  for 
the  production  of  CHiO. 

The  chemistry  for  describing  CH,  is  pro¬ 
vided  by  reactions  57-62  [14].  These  reactions 
were  first  used  in  the  modeling  of  CH  4/air 
flames  by  Coffee  [9].  Reactions  for  the  forma¬ 
tion  of  ground  state  CH  and  electronically 
excited  CH*  (A*A)  come  from  the  acetylene 
flame  mechanism  assembled  by  Eraslan  and 
Brown  [3]  from  several  sources. 

The  ionic  mechanism  consists  of  reactions 
74-86.  The  precursor  ion  CHO^  is  assumed  to 
be  produced  by  reaction  between  ground-state 
or  electronically  excited  CH  with  oxygen  atoms. 
In  addition  to  the  precursor  ion  CHO^,  the 
model  simulates  the  major  ions  HjO^,  CH3*, 
CtHjO"^,  and  C3H3’’.  Experiments  have 
shown  other  ions  to  be  present  at  concentra¬ 
tions  that  are  relatively  small  [17]  and  are 
neglected  in  the  present  study.  The  H30'^  ion 
is  assumed  to  be  produced  by  [4]: 

CHO^-l-  H3O  CO  -1-  H30^. 

The  CH3‘^  ion  is  assumed  to  be  produced  by 
ion-molecule  reaction  of  CH3  with  either 
CHO"^  or  H30‘'.  The  C2H30'^  ion  is  formed 
by  three  ion-molecule  reactions  (reactions  77, 
81,  82).  Brown  and  Eraslan  [3]  proposed  that 
C3H3“  is  formed  by 

CH3"+  C3H3  ^  C3H3^+  H3. 
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TABLE  1 


Reaction  Mechanism  for  Ions  in  Methane  Flames^ 

Reaction 

A 

B 

E. 

1 

OH  +  H,  H,0  +  H 

LITE  +  09 

UOOO 

1825.0 

2 

H  +  O;  i  OH '+  0 

L42E  +  14 

0.0000 

8250.0 

3 

0  +  Hj  s*  OH  +  H 

1.80E  +  10 

1.0000 

4480-0 

4 

H  +  0-,  +  M'  ^  HO,  +  M' 

L03E  +  18 

-0.7200 

0.0 

5 

H  +  HO,  OH  +  OH 

L40E  +  14 

0,0000 

540.0 

6 

H  +  HO,  0  +  H,0 

l.OOE  -h  13 

0.0000 

540-0 

7 

HO,  +  H  0,  +  H, 

L25E  +  13 

0.0000 

0.0 

8 

HO,  +  OH  =-  6,  +  H,0 

7oOE  +  12 

0.0000 

0.0 

9 

HO,  -r  0  0,  +  OH 

L40E  -r  13 

0.0000 

540.0 

10 

H  4-’  H  +  H,  «  H,  +  H, 

9.20E  -r  16 

-0.6000 

0.0 

11 

H  +  H  +  0,  ^  H,  +  O,' 

l.OOD  -f  18 

- 1-0000 

0.0 

12 

H  -r  H  +  h',0  ^  H,  +  H,0 

6.00E  +  19 

-1.2500 

0.0 

13 

H  -r  H  +  CO  H,  '-r  CO’ 

l.OOE  +  18 

-1.0000 

0.0 

14 

H  4-  H  +  CO,  H,  ^  CO, 

5.49E  +  20 

-2.0000 

0.0 

15 

H  4-  H  4-  CH4  H,  4-  CH^ 

5.49E  +  20 

-2.0000 

0.0 

16 

H  4-  OH  4-  M"  =«  h',0  4-  M’ 

L60E  +  22 

-2.0000 

0.0 

17 

H  4-  0  4-  M'  OH ’4-  M" 

6.20E  +  16 

-0.6000 

0.0 

18 

OH  4-  OH  =-«  H,0  4-  0 

5.72E  +  12 

0.0000 

390.0 

19 

OH  4-  CO  -  CO,  +  H 

L57E  -r  07 

uooo 

-385.0 

20 

0  4-  CO  4-  M'  CO,  4-  M' 

5.40E  +  15 

0.0000 

2300.0 

21 

H  4-  CO  4-  M'  5-  CHO  ^  M' 

5.00E  ^  14 

0.0000 

755.0 

CH4  4-  0  CH,  -r  OH 

4.07E  14 

0.0000 

7040.0 

23 

CH,  ^  H  -  CHj  +  H, 

7.24E  ^  14 

0.0000 

7590.0 

24 

CH,  4.  OH  -  CHj  -  H,0 

L.55E  ^  06 

2.1300 

1230.0 

25 

CH,  -4  M  =-  CHj  -  H  M 

4.68E  ^  17 

0.0000 

46910.0 

26 

CH,  -  0  CH,0  4-  H 

6.02E  -  13 

0-0000 

0-0 

27 

CH,0  -4  0^  CHO  -  OH 

L82E  -r  13 

0.0000 

1550.0 

28 

CH,0  -4  H  -  CHO  -  H, 

3.31E  ^  14 

0.0000 

5290.0 

29 

CH,0  -4  OH  *  CHO  4-  H,0 

7.58E  ^  12 

0.0000 

72.0 

30 

CHO  4-  0,  *>«  CO  4-  HO, 

3.00E  ^  12 

0.0000 

0.0 

31 

CHO  4-  h'=-  CO  t  H, 

4.00E  ^  13 

0.0000 

0.0 

32 

CHO  4-  OH  ^  CO  H,0 

5:00E  ^  12 

0.0000 

0.0 

33 

CHO  4-  0  *•  CO  4-  OH* 

l.OOE  ^  13 

0-0000 

0.0 

34 

CH,0  -4  CHj  ;=-  CHO  4-  CH, 

2.23E  ^  13 

0.0000 

2590.0 

35 

CHj  +  OH  5-  CH,0  4-  H, 

3-98E  ^  12 

0,0000 

0.0 

36 

CHj  4-  HO,  5-  CH4  4-  0, ' 

1.02E  +  12 

0.0000 

200.0 

37 

CO  +  HO,  5*  CO,  4-  OH 

L50E  -f-  14 

0.0000 

11900.0 

38 

CHj  4-  CHj  ^  C,'H4 

4.56E  -r  17 

-7.6500 

4250.0 

39 

CjH^  4-  0  =-  C,Hj  -4  OH 

2J1E  +  13 

0.0000 

3200.0 

40 

C,Hj  4-  H  CjH,  +  H, 

5.00E  +  02 

3J000 

2620.0 

41 

CjHj  4-  OH  ^  CjHj  4-  H,0 

6.63E  4-  13 

0.0000 

675.0 

42 

CjHj  4-  H  s*  C,Hs 

7.23E  +  13 

0.0000 

0.0 

43 

CjH,  4-  H  =•  CHj  4-  CHj 

3.73E  +  13 

0.0000 

0.0 

44 

CjHj  ^  C,H,  4-  H 

2.29E  -r  11 

0.0000 

19120.0 

45 

C,Hj  4-0,5*  C2H4  -4  HO, 

1J3E  +  12 

0.0000 

2446.0 

46 

CjHj  4-  0  5*  CH,  4-  CH,o' 

153E  4-  13 

0.0000 

2516.0 

47 

C,H4  4-  oh  CH,0  4.  CHj 

5.00E  +  13 

0.0000 

3020.0 

48 

C,H,  4-  0  C,Hj'-4  oh 

2.53E  4-  13 

0.0000 

2516.0 

49 

C,H4  4-  0,  5-  C,Hj  -4  HO, 

1J3E  +  15 

0.0000 

27680.0 

50 

CtH^  4-  H  5*  C,Hj  4“  H, 

lOOE  4.  15 

0.0000 

10000.0 

51 

CjH*  4-  OH  5*  CjHj  4-  H,0 

4.40E  +  14 

0.0000 

3270.0 

52 

CjHj  4-  M  5=5  C,H,  -4  H  4-  M 

3.01E  +  16 

0.0000 

20380.0 

53 

C,Hj  4-  0,  »  C,H,  -4  HO, 

U7E  4-  13 

0.0000 

5030.0 

54 

C,Hj  4-  H  ^  "4  H, 

7.i3E  4-  13 

0.0000 

0.0 

55 

C,Hj  -4  OH  5-  C,H,  -4  H,0 

l.OOE  -r  13 

0.0000 

0.0 

56 

C,H,  -r  OH  -  CHj  *4-  CO* 

5.48E  ^  13 

0.0000 

6890.0 
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TABLE  1 — Continued 


Reaction 

A 

B 

E. 

57 

CHj  +  H  -  CHj  -1-  Hj 

2.00E  +  11 

0.7000 

- 1500.0 

58 

CHj  OH  --  CHj  -1-  HjO 

6.00E  10 

0.7000 

1010.0 

59 

CHj  +  0,  CHO  -K  OH 

l.OOE  +  14 

0.0000 

1860.0 

60 

CHj  Oj  CH,0  +  O 

l.OOE  +  14 

0.0000 

1860.0 

61 

CHj  Oj  ••  COj  -1-  Hj 

l.OOE  +  14 

0.0000 

1860.0 

62 

CHj  +  H  CH  -I-  Hj 

4.00E  +  13 

0.0000 

0.0 

63 

CH  +  0  ^  CO  +  H 

4.00E  +  13 

0.0000 

0.0 

64 

CH  +  0, CO  +  OH 

lOOE  +  13 

0.0000 

0.0 

65 

CjH  +  6  CO  +  CH 

l.OOE  +  13 

0.0000 

0.0 

66 

CH*  +  M  ^  CH  +  M 

4.00E  +  10 

0.5000 

0.0 

67 

CH*  +  Oj  -  CH  -1-  0, 

2.40E  +  12 

OJOOO 

0.0 

68 

CH»  ^  CH 

1.90E  +  06 

0.0000 

0.0 

69 

C,H  +  Oj  CH*  +  CO, 

4.50E  +  15 

0.0000 

25000.0 

70 

C,H  +  0  CH*  -i-  CO 

7.10E  +  11 

0.0000 

0.0 

71 

C,H,  +  H  -  C,H  -h  H, 

6-OOE  +  13 

0.0000 

23651.0 

72 

C^H^  +  OH  C,H  -  H,0 

l.OE  13 

0.0000 

7000.0 

73 

cIh '+  0,  -  CO  +  CHO 

5.00E  +  13 

0.0000 

1505.0 

74 

CH  +  0  ^  CHO*-!-  e* 

2.52E  +  11 

0.0000 

1700.0 

75 

CH*  +  0  »■  CHO*+  e- 

5.04E  14 

0.0000 

1700.0 

76 

CHO*+  HjO  HjO*-!-  CO 

l.OOE  ^  16 

-0.0897 

0.0 

77 

HjO^-F  C,H,  C,HjO’-!-  H, 

8.39E  15 

0.0000 

0.0 

78 

CHO*-^  CH,'-  CHj^t  CO 

5.62E  ^  14 

-0.0060 

0.0 

79 

HjO*-r  CH,  -  CHj*-)-  H,0 

6.17E  -  14 

-0.0060 

0.0 

80 

CH3*-rC,H, -CjHj'-i-'H, 

7.24E  -  14 

0.0000 

0.0 

81 

CjHj'-r  H,6  -  C,H;,0"-'CH, 

7.24E  ^  14 

0,0000 

0.0 

82 

CHj'-r  CO;  -  C,H,6^^  O 

7.24E  14 

0.0000 

0.0 

83 

HjO'-r  H,6  H 

2.29E  18 

-0.5000 

0.0 

84 

C3H3'-!-e"—  products 

1,50E  -  19 

-0.5000 

0.0 

85 

CH3'-r<--  CH,  H 

2.29E  ^  18 

-OJOOO 

0.0 

86 

C,H30*+e--  M  -r  H 

2.29E  -r  IS 

-OoOOO 

0.0 

Reaction  rates  are  in  the  form  k  »  AT^zxpi-E^/RT)  and  have  units  of  cm^/mol-s. 


There  are  four  possible  mechanisms  for  the 
recombination  of  ions  in  flames.  Three  of  these 
mechanisms,  three-body,  radiative,  and  mutual 
neutralization  recombination,  have  been  ruled 
out  by  previous  research  [4]  as  unimportant 
pathways  for  ionic  recombination.  The  fourth 
and  most  widely  accepted  mechanism  is  disso¬ 
ciative  recombination  represented  by  reactions 
83-86  in  Table  1. 

Reaction  rate  coefficients  have  been  mea¬ 
sured  for  only  a  limited  number  of  ion-mole¬ 
cule  reactions.  Accordingly,  coefficients  have 
been  estimated  from  the  average  dipole  orien¬ 
tation  (ADO)  theory  [18-20]  for  reactions  for 
which  experimental  data  do  not  exist.  The  ADO 
theory  predicts  a  temperature-dependent  rate 
coefficient  that  is  correlated  to  the  parametric 
rate  expression 

k  =AT''. 

The  values  of  n  range  from  0.006  to  0.09. 


The  heats  of  formation  for  ions  have  proved 
very  imponant  in  predicting  ionic  structure  of 
flames  [3].  The  values  used  in  the  present 
simulation  are  given  in  Table  2. 

MODEL  FORMULATION 

The  flame  simulated  in  this  work  is  a  one-di¬ 
mensional,  laminar,  premixed  flame.  The  equa¬ 
tions  describing  the  structure  of  this  flame  in 
the  presence  of  an  electric  field  include  the 

TABLE 2 


Heats  of  Fonnation  for  Ions 


Species 

(kcai/mol) 

Ref. 

H3O* 

139.0 

21 

CHj^ 

256.0 

22 

CHO" 

199.1 

21 

C.Hj* 

255.0 

23 

C.H3O" 

151.0 

23 
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mass  and  species  conservation  equations  and 
Poisson's  equation  to  describe  the  electric  field 
intensity. 

Mass  conservation: 


dp  d{  pu) 

—  +  -  =  0. 

dt  dx 

Species  conservation: 


dt 


pu 


d_^ 

dx 


1  d(AJ,) 


dx 


+  O):. 


(1) 


(2) 


The  neutral  and  ionic  species  are  assumed  to 
obey  the  boundary  conditions  [3]: 

y  =  Yo  at  x  =  0  (3) 

and 


dY 

—  =  constant  at  x  =  3=. 


The  unsteady  terms  in  the  conservation  equa¬ 
tions  are  retained  as  a  computational  device 
only.  Their  inclusion  allows  steady-state  solu¬ 
tions  to  be  solved  by  time  marching,  an  ap¬ 
proach  proposed  by  Spalding  and' Stephenson 
[25]. 

Poisson’s  Equation: 
d^V  1 


Substituting  the  current  flux 


j  Een^p.^  (9) 

into  Eq.  8  gives  an  expression  for  the  force  per 
unit  volume  in  terms  of  mobility  and  current: 


(10) 


Multiplying  Eq.  10  by  the  distance  from  the 
flame  to  an  electrode  will  give  an  expression 
for  the  pressure  drop  across  the  flame: 


Ap  =jx\ - 

U- 


(11) 


Lawton  and  Weinberg  [17]  predict  the  maxi¬ 
mum  AR ,  from  Eq.  11,  to  be  on  the  order  of 
0.0004  atm.  Accordingly,  ionic  wind  effects  can 
be  neglected  in  the  flame. 

The  temperature  field,  in  principle,  can  be 
solved  from  an  energy  equation  applied  to  the 
flame.  In  practice,  however,  unaccountable  heat 
loss  from  the  flame  can  greatly  influence  the 
temperature  field,  making  prediction  difficult 
in  some  circumstances.  We  have  chosen  to  use 
e.xperimental  temperature  profiles  in  place  of 
the  energy  equation  in  our  simulations. 

Mass  flux,  Jj,  is  obtained  from  a  simplified 
expression  proposed  by  Wamatz  [26]  that  we 
have  modified  to  include  electric  mobility  ef¬ 
fects: 


The  boundary  conditions  for  Poisson’s  Equa¬ 
tion  applied  to  the  flame  are 


v  = 

and 

0 

at  X  =  0 

(6) 

r-t 

II 

(7) 

where  x  =  0  and  L  are  the  positions  of  the 
two  electrodes  in  the  one-dimensional  flame. 

The  pressure  across  a  flame  is  usually  taken 
to  be  constant.  Accordingly,  the  momentum 
equation  is  usually  not  solved  in  laminar  flame 
simulations.  However,  the  imposed  electric  field 
across  the  flame  gives  rise  to  a  force  on  charge 
earners.  The  neglect  of  the  momentum  equa¬ 
tion  must  be  justified  in  our  simulations. 

The  net  force  per  unit  volume,  F,  produced 
by  the  electric  field  is  (17) 

F  =  Ee(n^-n_).  (8) 


=  -pA., 


dx 


-  D 


i.T  ' 


dlnT 

17 


~  Fi. 


pY.E, 

(12) 


where  the  mixture  diffusivity  of  species  /,  D,- 
is  calculated  from  binary  diffusivities,  Z?,^: 


D,  „  = 


1 

nspc  V 

E  — 

y-i  A.. 


(13) 


and  the  thermal  diffusivity  of  species  i,  D,  j-,  is 
calculated  from 


A.r  = 


1  -  y 


j-i 

j*i 


(14) 


Wamatz  [26]  has  shown  that  the  first  two  terms 
in  Eq.  12  are  accurate  to  within  a  few  percent 
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of  the  exact  solution  of  the  Stefan-Maxwell 
equations.  In  this  study,  thermal  diffusion  is 
calculated  only  for  H  and  H2  because  these 
small  species  are  the  only  two  with  significant 
thermal  diffusion  rates. 

The  binary  mobilities,  can  be  related  to 
binary  diffusivities  by  the  Einstein  relationship 
[27]: 

A/  KT 

_ii  =  _i_  =  7  X  8.617  X  10-*.  (15) 


Little  experimental  data  are  available  on  the 
diffusion  of  ions.  It  is  assumed  that  the  diffu- 
sivity  of  ions  is  equal  to  that  of  the  correspond¬ 
ing  neutral  species.  From  these  assumed  ion 
diffusion  coefficients,  the  mixture  mobility, 
/i,,  can  be  determined  [27]: 

cm-‘V-'s-'.  (16) 


The  mobility  of  electrons  in  nitrogen  at  1000 
K  is  approximately  [27] 


Mr.  m 


nspc  Y 

r  — 

1 


2.0  X  10- 


cm  *  V“'  s“'. 


n_ 


(17) 


where  /i  _  is  the  electron  density.  The  diffusion 
coefficient  for  electrons  can  then  be  deter¬ 
mined  by  rearranging  the  Einstein  relation¬ 
ship: 

D,  =  8.93  X  10^  cmVs.  (18) 

The  diffusivity  is  assumed  to  scale  with  tem¬ 
perature  to  the  3 /2  power  [17].  Thus,  electron 
diffusivity  as  a  function  of  temperature  is 

D,  =  0.2824  X  T'-5  cmVs.  (19) 


In  the  absence  of  more  complete  data,  this 
relationship  is  used  to  estimate  electron  diffu¬ 
sion  in  this  work. 

Calculation  of  ion  diffusion  currents  is  com¬ 
plicated  by  the  self-induced  electric  field  that 
develops  among  the  charge  carriers  due  to  the 
diffusion  of  electrons  and  positive  ions  away 
from  each  other  [17].  This  electric  field  im¬ 
pedes  the  further  diffusion  of  charge  carriers, 
with  the  result  that  all  charge  carriers  can  be 


described  by  a  common  diffusion  coefficient 
called  the  ambipolar  diffusion  coefficient, 
[17]: 


A 


D-m- 


(20) 


where  the  -i-  and  -  subscripts  denote  proper¬ 
ties  of  positive  and  negative  charge  carriers, 
respectively.  However,  this  description  of  ion 
diffusion  is  approximate  and  inconvenient  when 
treating  a  large  number  of  different  charge 
carriers.  We  prefer,  instead,  to  explicitly  calcu¬ 
late  the  self-induced  electric  field  arising  from 
charge  separation.  In  this  approach,  ion  diffu¬ 
sion  and  ion  mobility  are  distinct  processes. 

This  calculation  procedure  consists  of  solv¬ 
ing  Poisson’s  equation  (Eq.  5)  to  find  the  elec¬ 
tric  field  intensity  generated  by  the  distribution 
of  charge  carriers  at  a  given  instance  in  time. 
The  resulting  electric  field  intensity  E  is  used 
to  calculate  the  contributions  of  ion  currents 
to  the  overall  mass  flux  of  individual  ions  as 
given  by  Eq.  12.  The  mass  fluxes  J,  of  individ¬ 
ual  ions  is  substituted  into  Eq.  2  and  solved  for 
the  distribution  of  charge  carriers  during  the 
next  time  step.  Poisson’s  equation  is  again 
solved  to  find  the  electric  field  intensity  corre¬ 
sponding  to  the  new  charge  distribution.  Like 
chemical  species,  the  electric  field  intensity 
settles  down  to  a  steady-state  profile  after  a 
sufficient  number  of  time  steps. 

Solution  of  the  continuity  and  species  con¬ 
servation  equations  is  greatly  simplified  by  the 
coordinate  transformation  [28]: 

_.p  and  ---pn.  (21) 

Overall  continuity  is  automatically  satisfied  and 
the  neutral  species  conservation  equation  be¬ 
comes 


d<b 

d 

f  ^S4>] 

P — 

=  n— — 

pr— 

dt 

dip 

(22) 


where  is  a  generic  variable  representing 
mass  fractions  YJ.  F  and  G  for  neutral  species 
are  defined  in  Table  3  [28].  The  transformed 
equation  for  ionic  species  conservation  is 

d  I  d4>  \  ^ 

p—  =  P—  ^ 

<^f  d{If\  dip  I  . 

(23) 
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where  the  expressions  corresponding  to  F,  F', 
and  G  for  ions  are  also  given  in  Table  3.  The 
appropriate  boundary  conditions  for  the  trans¬ 
formed  coordinates  are 

<f>  —  <pQ  at  i//  =  0  (24) 

and 

d<b 

—  =  constant  at  »  (25) 

dill 

Solving  Eqs.  22  and  23  using  the  above 
boundary  conditions  gives  the  solution  in  the 

coordinate  system.  A  simple  transformation 
back  to  the  x  coordinate  system  is  obtained 
from 

rii  diij 

x  =  f  — .  (26) 

The  resulting  system  of  stiff  partial  differen¬ 
tial  equations  is  solved  by  split  operator  tech¬ 
niques  [29,  30]  and  a  specialized  method  of 
lines  algorithm  described  in  detail  in  Ref.  28. 

Poissons  equation  (Eq.  5)  must  be  solved 
periodically  as  the  electron  and  ionic  species 
profiles  evolve  in  time.  After  one  complete 
iteration  on  all  species  conservation  equations, 
Poisson’s  equation  is  solved.  The  numerical 
solution  of  Poisson’s  equation  is  simplified  by 
transforming  it  into  a  pair  of  first  order  ordi¬ 
nary  differential  equations.  This  is  done  by 
replacing  voltage  with  electric  field  intensity: 


dE  1 

—  =  — (rt.-  n_)e, 
dx  €0 

(27) 

dV 

—  = 

dx 

(28) 

with  V{x  =  0)  and  V{x  ■=  D  as  prescribed 
boundary  conditions.  These  boundary  values 
are  zero  when  no  external  electric  field  is  ap¬ 
plied  to  the  flame.  The  resulting  boundary 


TABLES 

Values  of  r,  r\  and  G  for  the  Species 
Conservation  Equations 


Equation 

r 

r 

G 

Neutral  model 

PD,.„ 

— 

Ionic  model 

P^i.m 

value  problem  is  solved  using  a  shooting  meth¬ 
od. 

Flow  simulations  are  performed  on  a  HDS 
AS/9160.  The  number  of  grid  points  used  for 
the  simulations  range  from  50  to  70.  The  aver¬ 
age  CPU  time  required  to  reach  steady  state 
was  20  min  for  a  single  run. 

RESULTS  AND  DISCUSSIONS 

The  focus  of  this  study  is  simulation  of  ionic 
species  profiles  and  saturation  currents  in 
methane  flames.  However,  the  importance  of 
accurate  neutral  species  profiles  for  the  predic¬ 
tion  of  ion  profiles  in  flames  cannot  be  over¬ 
stated.  Unfortunately,  most  flame  studies  have 
concentrated  on  measurement  of  either  ion  or 
neutral  profiles  to  the  exclusion  of  the  other. 
We  have  been  careful  to  select  methane  flames 
for  which  experimental  data  exists  on  both 
neutral  and  ionic  species  concentration  pro¬ 
files.  We  found  experimental  data  at  three 
different  equivalence  ratios,  6,  representing 
lean,  stoichiometric,  and  rich  methane  flames 
that  are  complete  enough  to  compare  with  our 
model  predictions.  We  also  found  limited  data 
on  a  hydrogen  flame  seeded  with  1.8%  meth¬ 
ane.  For  each  flame,  we  compare  model  pre¬ 
dictions  to  experimental  data  for  both  neutral 
and  ionic  species  concentration  profiles  to  vali¬ 
date  our  kinetic  mechanisms.  These  validation 
comparisons  are  made  in  the  absence  of  exter¬ 
nally  applied  electric  fields.  After  validating 
the  ionic  mechanism  of  our  model,  we  per¬ 
formed  simulations  with  externally  applied 
electric  fields  to  investigate  the  phenomenon 
of  saturation  currents. 

Lean  Flame  ( p  =  40  torr,  cj)  =  0,2) 

The  experimental  data  for  the  neutral  species 
profile  comes  from  the  work  of  Peeters  and 
Mahnen  [31].  The  CH4,  O2,  HiO,  CO,  and 
CO2  shown  in  Fig.  la  and  lb  are  in  good 
agreement  with  the  experimental  data.  The 
concentration  profiles  of  the  minor  neutral 
species  H,  H,,  O,  HO,,  and  CHtO,  shown  in 
Figs.  Ic  and  Id,  are  within  the  uncertainty  of 
experimental  data.  The  model  tends  to  over¬ 
predict  the  experimental  data  in  the  post  flame 
region,  but  ionic  structure  is  most  strongly 
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Fig.  1.  Predicted  and  experimental  neutral  species  profiles  for  a  40  torr,  d  —  0.2  methane  flame  (lines  are  model 
predictions,  symbols  are  experimental  data  from  Peeters  and  Mahnen  [31]).  (a)  CH^,  CO,  and  COj.  (b)  H,0  and  O;.  (c)  H, 
H;,  0,  and  OH.  (d)  HO,  and  CH,0. 


influenced  by  conditions  in  the  flame  front 
where  model  predictions  are  good. 

The  results  of  Olsson  and  Andersson’s  [32] 
sensitivity  study  on  a  lean  (4>  =  0.2),  low  pres¬ 
sure  (40  torr)  methane -oxygen  flame  is  rele¬ 
vant  to  our  attempts  to  simulate  neutral  species 
profiles.  Their  smdy  indicated  that  small  per¬ 
turbations  in  the  reaction  rates  for  reactions 
that  produce  CO  from  CHO  could  result  in 
variations  of  up  to  50%  in  some  species  con¬ 
centrations.  TTie  concentration  profiles  for 
CHj,  H,  Hj,  O,  and  HOj  were  found  to  be 
sensitive  to  changes  in  the  reaction  rate  for 

CHO  -4  M  -*>  H  +  CO  -1-  M. 

The  reaction  rate  for  this  particular  reaction 
has  a  high  degree  of  uncertainty  associated 
with  it  [33].  Olsson  and  Andersson  [32]  showed 
that  increasing  the  reaction  rate  by  a  factor  of 
five,  well  within  its  range  of  uncertainty,  pro¬ 
duced  changes  in  the  concentrations  of  the 
aforementioned  species  in  excess  of  50%.  Con¬ 


sidering  the  sensitivity  of  the  species  profiles  to 
the  rate  of  this  reaction,  it  can  be  concluded 
that  the  predicted  concentration  profiles  are 
within  the  uncertainty  of  experimental  data. 

Goodings  et  al.  [34]  used  a  mass  spectrome¬ 
ter  to  determine  ion  currents  of  several  ionic 
species  as  a  function  of  position  in  the  flame 
front  of  methane/oxygen  flames  at  40  torr  and 
an  equivalence  ratio  of  0.2.  The  experimental 
profiles  are  shown  in  Fig.  2a.  The  ion  concen¬ 
tration  profiles  predicted  by  the  model  for  this 
flame  are  shown  in  Fig.  2b.  Although  direct 
comparison  between  ion  currents  and  ion  con¬ 
centrations  cannot  be  made,  several  qualitative 
observations  are  possible  from  the  figures. 

The  model  predicts  to  be  the  domi¬ 

nant  ion  with  an  ion  density  of  5  x  10‘°/cm^. 
The  peak  value  of  this  ion  occurs  slightly 
downstream  of  the  peak  concentrations  for  the 
other  ions.  The  slower  decay  of  HjO"*"  is  at¬ 
tributed  to  the  reaction 

CHO^-4  H,0  ^  H30*+  CO, 
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Fig.  2.  Predicted  and  e.xperimental  ionic  species  profiles 
for  a  40  torr.  <i>  =  0.2  methane  flame,  (a)  Expenmentai 
data  from  Goodings  et  al.  [34].  (b)  Model  predictions. 

which  produces  the  majority  of  the  HjO"  ions 
in  a  region  where  there  is  a  rapid  increase  in 
the  production  of  H,0.  The  other  major  ions 
decay  more  rapidly  at  rates  that  appear  to  be 
comparable  for  all  ions.  This  rapid  decay  oc¬ 
curs  because  the  reaction  intermediates  that 
form  these  ions  by  ion-molecule  reactions  are 
depleted  very  rapidly  in  the  flame.  This  general 
behavior  is  observed  in  the  experimental  data 
of  Goodings  et  al.  [34]  shown  in  Fig.  2a.  Cal- 
cote  [35]  has  noted  that  this  slow  decay  of 
and  rapid  decay  of  the  other  major  ions 
characterizes  all  hydrocarbon  flames. 

Since  ion  currents  measured  by  mass  spec¬ 
trometry  cannot  be  easily  converted  into  ion 
concentrations,  we  tabulate  relative  ion  con¬ 
centrations  in  Table  4  for  a  more  quantitative 
comparison  between  model  predictions  and  ex¬ 
periment.  The  peak  fractional  concentrations 
predicted  for  HjO^,  C2H30~,  and  CHj' 
compare  reasonably  well  with  the  correspond¬ 
ing  experimental  values.  The  ion.  which 

was  not  detected  experimentally,  was  predicted 


TABLE  4 


Relative  Peak  Ion  Concentrations  for  a  CH^/O,  Flame 
( p  =  760  torr,  <(>  —  1.0) 


Species 

ExperimenLal 

Ratio'' 

Predicted 

Ratio"* 

0.804 

0.770 

C.H3O" 

0.158 

0.182 

CH3" 

0.023 

0.011 

CHO* 

0.004 

0.036 

C3H3" 

— 

0.0001 

“  Referenced  to  total  ion  concentration. 


to  appear  at  concentrations  only  1/1000  of 
that  of  the  predominant  ion,  C3H3O".  The 
major  discrepancy  in  the  model  was  overpre¬ 
diction  of  CHO"  by  a  factor  of  ten.  Very  likely 
this  overprediction  is  the  result  of  our  exclu¬ 
sion  of  minor  ionic  species  in  the  model,  which, 
taken  together,  could  significantly  deplete 
CHO",  the  precursor  of  larger  ions. 

Stoichiometric  Flame  ip  =  760  torr,  <b  =  1.0) 

The  experimental  data  for  neutral  species  are 
taken  from  Bechtel  et  al.  [36]  for  a  flame  at 
atmospheric  pressure  and  an  unbumed  gas 
temperature  of  298  K.  The  predicted  and 
experimental  profiles  for  O,,  H^O,  and  tem¬ 
perature  are  shown  in  Fig.  3a.  The  model 
accurately  predicts  the  measured  O,  and  H.,0 
profiles.  Figure  3b  shows  that  the  H,  and  CO 
concentrations  are  correctly  predicted  to  peak 
in  the  flame  front.  The  general  trend  of  the 
predicted  profiles  agree  with  the  experimental 
data.  However,  the  predicted  peak  concentra¬ 
tion  of  CO  falls  short  of  the  experimental 
value.  Bechtel  et  al.  [36]  indicate  that  the  error 
associated  with  the  peak  value  of  CO  is  ±  10%. 
Therefore,  the  discrepancy  between  predic¬ 
tions  and  experiments  is  within  experimental 
uncenainty. 

The  experimental  ionic  profiles  for  the  stoi¬ 
chiometric,  atmospheric  methane-oxygen  flame 
come  from  the  work  of  Peeters  et  al.  [37]. 
Figure  4a  contains  the  experimental  ion  con¬ 
centration  profiles  while  Figure  4b  contains  the 
predicted  ion  concentration  profiles.  Similar  to 
the  lean  flame  results,  the  HjO^  ion,  in  both 
model  predictions  and  experimental  data,  peaks 
late  and  decays  slowly  compared  with  other 
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fb) 

Fig.  3.  Predicted  and  experimental  neutral  species  profiles 
for  a  760  torr,  (b  *  1,0  methane  flame  (lines  are  model 
predictions,  symbols  arc  experimental  data  from  Bechtel 
et  ai.  [36]),  (a)  H2O,  O^,  and  flame  temperature,  (b)  CO 
and  Hn. 


(b) 


Fig.  4.  Predicted  and  experimental  ionic  species  profiles 
for  a  760  torr.  <i>  »  I.O  methane  flame,  (a)  Experimental 
data  from  Peeters  et  ai,  [37],  (b)  Model  predictions. 


flame  ions.  The  experimental  ion  profiles  are 
notably  broader  than  the  predicted  profiles, 
but  in  the  narrow  flame  front  of  an  atmo¬ 
spheric  flame  this  discrepancy  may  well  repre¬ 
sent  experimental  difficulties  in  spatial  resolu¬ 
tion. 

More  revealing  are  comparisons  of  the  peak 
concentrations  for  the  various  ions,  shown  in 
Table  5.  The  model  correctly  predicts  the  rela¬ 
tive  ordering  of  ions  by  abundance.  Further¬ 
more,  predicted  peak  concentrations  for  the 
four  major  ions  HjO^,  CHO"^,  and 

are  within  a  factor  of  four  of  their 
experimental  counterparts.  The  final  flame 
temperature  used  in  the  simulation  is  2240  K, 
which  comes  from  the  measured  temperature 
profile  of  Bechtel  et  al.  [36].  Peeters  et  al.  [37] 
document  a  range  of  final  flame  temperatures 
varying  from  1830  to  2084  K.  In  the  same 
paper,  Peeters  et  al.  [37]  show  the  ion  yield  to 
increase  approximately  20%  per  75  K  increase 
in  temperature.  Therefore,  differences  in  flame 


temperatures  between  simulation  and  experi¬ 
ment  can  account  for  much  of  the  discrepan¬ 
cies  in  the  ion  concentrations. 


Rich  Flame  ip  =  760  torr,  <j>  =  2,13) 

The  rich  flame  experiments  of  Goodings  et  al. 
[23]  are  notable  in  that  data  for  both  ionic  and 
neutral  species  were  collected.  The  predicted 
and  experimental  neutral  species  profiles  for 
CH^,  O.,  H.,  and  CO  are  illustrated  in  Fig. 
5a.  The  predicted  O.  and  CH^  profiles  show 

TABLES 


Peak  Ion  Concentrations  for  a  Stoichiometric  CH^/02 
Flame  (;?  »  760  torr,  o6  *  1,0) 


Species 

Experimental 

Concentration 

(mol/cm^) 

Predicted 

Concentration 

(mol/cm^) 

HjO* 

1.5  X  10“ 

3.1  X  10“ 

C-HjO- 

4.3  X  10’ 

2.0  X  10’ 

C3H,* 

1.2  X  10’ 

1.5  X  10’ 

CHO' 

1.9  X  10’ 

6.1  X  10’ 
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POSmON  (cmi 


(C) 


Fig,  5,  Predicted  and  experimental  neutral  species  profiles 
for  a  760  torr,  6  ^  2.13  methane  flame  (lines  are  model 
predictions,  symbols  are  experimental  data  from  Goodings 
et  al.  [23D.  (a)  CH^,  CO,  H,,  and  O.,  (b)  CHj,  (c)  CO., 
C^Hj,  andH.O, 

good  agreement  with  the  experimental  data. 
The  post  flame  concentrations  of  major  stable 
products,  Ht,  H^O,  and  CO  are  well  pre¬ 
dicted,  although  CO,  in  panicular,  plateaus 
earlier  than  predicted.  We  attribute  some  dis¬ 
crepancies  in  the  preheat  zone  of  the  flame  to 
apparent  errors  in  the  temperature  profile  re¬ 
ported  by  Goodings  et  al.  [25]. 

The  model  predicts  CH3  concentration, 
shown  in  Fig.  5b,  that  is  a  factor  of  two  lower 
than  the  data  from  Goodings  et  al.  [23].  This 


underprediction  of  CHj  for  rich  flames  arises 
from  its  rapid  conversion  to  CiHj.  This  same 
trend  was  observed  in  the  rich  flame  simula¬ 
tion  performed  by  Coffee  et  al.  [9].  They  re¬ 
ported  an  underprediction  in  the  CH3  profile 
by  a  factor  of  three  compared  with  the  experi¬ 
mental  data  for  a  rich  flame.  However,  the 
actual  equivalence  ratio  used  by  Coffee  et  al. 
[9]  was  not  documented. 

The  predicted  value  of  CjHj,  shown  in  Fig. 
5c,  is  also  a  factor  of  two  lower  than  the 
experimental  data.  CHj  is  directly  responsible 
for  the  production  of  C^Hj,  thereby  causing 
the  factor  of  two  discrepancy  in  the  data.  These 
differences  can  be  attributed  to  an  incomplete 
understanding  of  the  intermediate  reactions 
responsible  for  the  oxidation  of  CH^. 

The  predicted  and  experimental  ionic 
profiles  for  this  rich  flame  are  compared  in 
Fig.  6.  The  relative  concentrations  of  H3O", 
CjHjO'',  and  C3H3’  are  correctly  predicted, 
with  H30“  and  C3H3~  being  the  dominant 
ions  in  the  flames.  Figure  6a  illustrates  that  the 
location  of  peak  concentrations  are  correctly 
predicted  for  and  C3H3''  although 

the  experimental  profiles  are  broader  than  pre¬ 
dicted.  The  H3  0^  ion,  illustrated  in  Fig.  6b, 
does  not  peak  downstream  of  the  other  posi¬ 
tive  ions,  as  is  expected  from  the  experimental 
data. 

The  CHO"  and  CH3'^  ions,  illustrated  in 
Fig.  6c,  are  poorly  predicted  in  terms  of  both 
peak  concentrations  and  profile  shapes.  Van 
Tiggelen  [38]  reports  that  for  rich  methane 
flames,  concentrations  of  CHj"  and  CHO'^ 
are  about  one  tenth  that  of  H3O".  Peeters  et 
al.  [65]  report  that  concentrations  of  CH3'^ 
and  CHO"^  are  a  factor  of  100  smaller  than 
that  of  H3  0‘^.  Goodings  et  al.  [32]  show  the 
ion  concentration  of  CH3'^  to  be  a  factor  of 
400  smaller  and  the  concentration  of  CHO  *  to 
be  a  factor  of  600  smaller  than  H30‘^.  The 
prediction  by  this  model  show  CH3‘^  and 
CHO"^  concentrations  to  be  three  orders  of 
magnitude  smaller  than  the  concentrations  of 
H30’'.  The  underprediction  of  these  two  ions 
can  be  traced  to  the  neutral  species  model, 
which  underpredicts  CH3  and  CH,  by  a  factor 
of  two.  These  two  species  participate  in  reac¬ 
tions  to  produce  CH  and  CH,,  the  precursors 
of  CHO~  and  CH3’’.  Although  Goodings  et  al. 
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Fig.  6.  Prediaed  and  experimental  ionic  species  profiles 
for  a  760  torr,  -  2.13  methane  flame  (lines  without 
symbols  are  model  prediaions,  lines  with  symbols  are 
experimental  data  from  Goodings  et  al.  [23].  (a)  C1H3O* 
and  CjHj*.  (b)  H3O".  (c)  CHO*  and  CH,*. 

[34]  did  not  report  concentrations  for  CH  and 
CH,,  their  origin  suggests  that  our  model  un¬ 
derpredicts  them,  giving  rise  to  underpredicted 
CHj^andCHOT 

The  experimental  profiles  for  CH3’’  and 
CHO'*’  reported  by  Goodings  et  al.  [34]  show 
trends  that  are  not  in  accord  with  observations 
of  most  other  researchers.  The  experimental 
data  shows  the  CHj"  ion  to  peak  before  the 
flame  front  with  a  slow  decay  that  extends 


throughout  the  post-flame  region.  The  experi¬ 
mental  CHO^  profile  exhibits  two  peaks,  one 
before  and  one  after  the  flame  front.  These 
discrepancies  may  be  the  result  of  expansion 
cooling  in  the  sampling  nozzle  during  the  ex¬ 
periment. 

Methane-Seeded  Hydrogen  Flame 

The  model  was  used  to  simulate  a  stoichiomet¬ 
ric,  atmospheric  hydrogen  flame  seeded  with 
1.8%  methane  reported  by  Peeters  et  al.  [37]. 
The  major  ion,  HjO"*",  is  predicted  to  peak  just 
after  the  flame  front  and  exhibits  a  slow  decay 
in  the  post-flame  region.  The  predicted  ion 
density  of  1.0  x  10"/cm^  for  agrees 

well  with  the  value  of  1.8  x  reported 

by  Peeters  et  al.  [37].  The  experimental  profile 
also  shows  a  slow  decay  in  the  postflame  re¬ 
gion. 

The  model  predicts  an  ion  density  on  the 
order  of  lOVcm^  for  th®  C^HjO"  ion,  which 
agrees  well  with  the  value  of  10®/ cm^  reported 
by  Peeters  et  al.  [37].  The  model  predicts  all 
other  ions  to  have  ion  densities  less  than 
10®/cm^.  The  CHO”,  CH3“,  and  C3H3~  ions 
did  not  appear  in  measurable  quantities  in  the 
flame  of  Peeters  et  al.  [37]. 

In  their  experiments,  Peeters  et  al.  [37]  mea¬ 
sure  a  large  concentration  of  CH5O”  in  the 
flame.  It  is  reported  to  have  an  ion  density  of 
2.5  X  10‘°/cm^.  The  reaaion  mechanism  used 
in  the  present  research  does  not  take  into 
account  the  formation  of  this  ion,  which  has 
not  been  reported  as  a  major  flame  ion  by 
other  researchers. 

Saturation  Currents 

Saturation  currents  were  determined  from  the 
simulations  in  a  manner  similar  to  that  which 
would  be  employed  in  actual  experiments.  A 
series  of  flame  simulations  with  specified  initial 
conditions  and  temperature  profile  were  run  in 
which  the  electric  field  voltage  Vq  applied  to 
the  electrodes  was  increased  from  one  run  to 
the  next  until  the  electric  current  flux,  ;, 
through  the  flame  became  constant.  This  cur¬ 
rent  was  recorded  as  the  saturation  current 
flux,  y,,  for  the  specified  flame  conditions.  The 
apparent  activation  energy,  for  the  satu- 
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ration  current  was  determined  from  the  slope 
of  the  line  produced  by  plotting  log( /,)  versus 
1/7}.  In  actual  experiments,  the  peak  flame 
temperature  7}  would  be  controlled  by  varying 
the  flow  rate  of  gas  mixture  to  the  burner. 
However,  as  previously  noted,  predicting  heat 
transfer  from  a  laboratory  burner  flame  is 
problematic  unless  more  information  is  pro¬ 
vided  about  burner  operation  than  is  typical  in 
the  literature.  In  the  absence  of  detailed  tem¬ 
perature  profiles  for  the  experimental  flames 
in  which  saturation  currents  were  measured, 
we  simply  applied  a  scaling  factor  to  known 
temperature  profiles  of  comparable  flames  to 
achieve  desired  peak  flame  temperatures  in 
our  simulations. 

Simulations  were  performed  for  the  lean, 
stoichiometric,  and  rich  flames  described  previ¬ 
ously  as  well  as  for  a  stoichiometric  hydrogen 
flame  seeded  with  1%  methane.  Apparent  acti¬ 
vation  energies  are  given  in  Table  6. 

The  saturation  current  is  predicted  to  peak 
at  an  equivalence  ratio  close  to  unity.  Al¬ 
though  we  do  not  have  experimental  data  for 
methane  flames  with  which  to  compare  this 
prediction,  it  is  supported  by  experimental  ob¬ 
servations  by  Peelers  et  al.  [37]  for  other  hy¬ 
drocarbon  flames.  This  behavior  reflects  the 
existence  of  a  maximum  in  ion  generation  rate 
at  stoichiometric  conditions. 

The  apparent  aaivation  energy  for  the  satu¬ 
ration  current  in  methane  flames  is  predicted 
to  range  between  41.4  and  46.8  kcal/mol.  Al¬ 
though  the  variance  in  the  experimental  data  is 
rather  large,  it  agrees  reasonably  well  with 
model  predictions. 


TABLE  6 


Activation  Energies  Obtained  from  Saturation  Currents 
for  Atmospheric,  Methane  Flames  and  a  1% 
Methane-Seeded  Hydrogen  Flame 


Flame 

Equivalence 

Ratio 

Experimental 

(kcai/raol) 

Predicted 

(kcai/mol) 

CH4/0, 

02 

— 

4 1. 36 

CH,/0,' 

l.O 

46.0 

46.8 

CH,/0- 

2.13 

— 

42o 

l.O 

2S.1 

42.6 

Lawton  and  Weinberg  [6,  17],  whose  value  of 
activation  energy  agrees  within  a  few  percent 
of  our  values,  were  unable  to  identify  the  rate- 
limiting  process  responsible  for  this  apparent 
activation  energy.  They  hypothesized  that  free- 
radical  reactions  might  play  a  role.  However, 
we  think  ion— molecule  reactions  are  a  better 
explanation. 

We  employed  sensitivity  analysis  to  trace  the 
origin  of  the  apparent  activation  energy  in  the 
simulations.  A  simple  but  effective  sensitivity 
analysis  systematically  changes  rate  constants 
to  evaluate  the  imponance  of  individual  chem¬ 
ical  reactions  on  some  overall  system  parame¬ 
ter.  Preexponential  factors  are  customarily  al¬ 
tered  to  affect  changes  in  rate  constants;  we 
chose  to  vary  rate  constants  through  changes 
in  activation  energies  of  individual  ion-mole¬ 
cule  reactions.  This  modified  approach  to  sen¬ 
sitivity  analysis  gave  us  a  direct  test  as  to 
whether  a  single  reaction  was  responsible  for 
the  apparent  activation  energy  of  ion  satura¬ 
tion  current.  We  determined  that  the  orisin  of 
the  apparent  activation  energy  is 

CtHjO  Hi  — »  HjO'-r  CiHi, 

which  has  an  activation  energy  of  42.19  kcal/ 
mol.  To  understand  the  surprising  rate-limiting 
action  of  this  reaction,  the  major  pathways  for 
charge  transfer  from  chemiionizaton  to  disso¬ 
ciative  recombination  were  traced.  Figure  7 
summarizes  this  analysis  for  stoichiometric 
methane  flames,  showing  major  reaction  path¬ 
ways  between  ions  by  bold  arrows  and  minor 
pathways  by  thin  arrows.  The  rate  limiting  ac¬ 
tion  of  the  reaction  appears  to  arise  from  the 
channeling  of  a  large  fraction  of  positive 


Fig.  7.  Positive  ion  pathways  in  a  i  -  t.O  methane  flame 
(heavy  arrowy  denote  pnmary  pathwavs). 
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charges  through  on  their  way  to  dis¬ 

sociative  recombination  by  the  reaction 

H30~+  «■-»  HjO  H. 

This  pathway  was  prominent  for  all  equivalent 
ratios  examined  in  this  study. 

Simulations  of  the  methane-seeded  hydro¬ 
gen  flame  produced  an  apparent  activation  en¬ 
ergy  comparable  to  the  methane  flames.  This 
prediction  is  in  accordance  with  the  fact  that 
chemiionization  arises  from  the  hydrocarbons 
in  the  flame  and,  hence,  should  closely  follow 
the  methane  reaction  mechanism.  However, 
the  activation  energy  obtained  e.’cperimentally 
by  Lawton  and  Weinberg  was  only  half  that 
predicted  by  the  simulation,  although  the  satu¬ 
ration  current  was  of  the  right  order  of  magni¬ 
tude.  We  see  no  plausible  e.xplanation  for  such 
a  small  activation  energy  in  methane-seeded 
hydrogen  flames;  funher  experimental  verifi¬ 
cation  of  this  activation  energy  would  appear 
to  be  in  order. 

CONCLUSIONS 

A  simple  ionic  reaction  mechanism  is  proposed 
that  predicts  the  general  features  of  ions  in 
lean  to  stoichiometric  methane  flames.  The 
proposed  reaction  mechanism  provides  a  means 
for  the  production  and  recombination  of  major 
ions  thought  to  exist  in  methane  flames.  The 
largest  discrepancy  between  predictions  and 
experiments  occurred  for  rich  flame  simula¬ 
tions  in  which  the  model  underpredicted  con¬ 
centrations  for  CHO"^  and  CH3  ions.  This 
discrepancy  is  attributed  to  neutral  species  re¬ 
action  mechanisms  associated  with  reaction  in¬ 
termediates. 

Using  this  ionic  reaction  mechanism,  the 
saturation  currents  for  lean,  stoichiometric,  and 
rich  methane  flames  were  calculated.  The  lim¬ 
ited  literature  on  saturation  currents  rarely 
reports  saturation  currents  as  functions  of  tem¬ 
perature.  This  deficiency  does  not  allow  for  a 
direct  comparison  between  the  experimental 
and  the  predicted  saturation  currents.  How¬ 
ever,  the  literature  does  report  apparent  acti¬ 
vation  energies  for  saturation  currents,  which 
we  compared  to  values  obtained  from  simula¬ 


tions.  The  predicted  values  range  from  41.36 
kcal/mol  for  the  lean  flame  to  46.81  kcal/mol 
for  the  stoichiometric  flame.  These  values  fall 
within  the  range  of  values  reponed  in  the 
literature  for  methane  flames. 

The  predicted  apparent  activation  energies 
were  used  to  determine  the  rate  limiting  reac¬ 
tion  in  the  ionic  reaction  mechanism.  We  have 
traced  the  activation  process  for  all  equiva¬ 
lence  ratios  to  the  ion-molecule  reaction: 

C2H30^+  H,  H30"^  C.H,. 

This  reaction  is  responsible  for  channeling  a 
large  fraction  of  positive  charge  to  dissociative- 
recombination  reaction. 

Some  of  this  work  was  performed  with  support 
from  the  Air  Force  Office  of  Scientific  Research 
(Contract  Mo.  F49620-91-C-0021)  while  under 
subcontract  to  Aerochem  Research  Laboratory. 
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